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WAVE FORCES ON SUBMERGED STRUCTURES 


Ernest F. Brater,! M. ASCE, John S. McNown,2 M. ASCE 
and Leslie D. Stair3 
(Proc. Paper 1833) 


SYNOPSIS 


The magnitude and characteristics of forces resulting from oscillatory 
javes were determined for models of submerged barge-like structures. The 
ave profiles and the variation with time of the resulting horizontal and 
ertical forces were obtained for various wave heights, wave periods, and 
ocations of the barge with respect to the water surface. The basic model 
jas a right parallelepiped having the proportions of a typical barge. Modi- 
ications of this model were made to determine the effects of rounding the 
ides and of the inclusion of a slot in one end of the barge. Horizontal force 
aeasurements were also made for a flat plate having the same dimensions 
s the long vertical face of the barge. These studies provided design data 
eeded for the application of analytical methods in the determination of forces 
n some types of submerged structures. 


INTRODUCTION 


The design of off-shore structures such as oil drilling platforms and 
adar stations requires knowledge of the magnitude of the wave forces. 
fethods of calculating such forces have not been fully determined. Labora- 
wry investigations have been made of wave forces on vertical piles, (1,2,3,4) 
orizontal struts(5,6) and on larger and more complicated structures. (7 
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o laboratory investigations of forces on simple barge-like shapes have been 
sported. 

The forces produced by oscillatory waves result from the orbital motion 
-the water. Acceptably accurate estimates can be made of orbital velocities 
id accelerations, and the associated pressure variations can be determined 
ialytically. Also required for the computation of the drag and inertial forces 
re values of the coefficient of drag (Cg) and the coefficient of inertial re- 
istance (Cy). The drag coefficients for various bodies have been determined 
<perimentally for steady flow, but little is known about their magnitudes for 
isteady flow. Values of C,, have been determined theoretically as well as 
<perimentally for certain shapes, but here again the applicability of such 
aulues to the conditions existing during wave motion has not been tested. 
hese studies were planned to determine the necessary coefficients and a 
litable method for estimating wave forces on structures placed in the sea. 


Testing Facilities 


The experiments were conducted in the Lake Hydraulics Laboratory(8) of 
ie University of Michigan. The wave tank used for these tests is 86 feet 
mg and 40 feet wide and has walls 3‘ 10"' high. The testing frame was 
cated in a channel 13' 10'’ wide separated from the remainder of the tank 
7 a temporary wall as shown in Fig. 1. Waves were created by a wave 
achine of the plunger type located at one end of the channel. The wave 
achine could be regulated to produce waves of various heights and periods. 
he model was centered in the channel approximately 45 feet from the 
unger. After passing the model, the waves were reflected out of the channel 
to the larger portion of the tank where their energy was dissipated by means 
an artifial beach (Fig. 1). Wave heights were measured by means of re- 
rding electric resistance gages. 

The testing frame was mounted on two heavy I-beams supported on the 
alls of the wave channel. A photograph showing the method of supporting 
e model and some of the instrumentation is shown in Fig. 2. Details of the 
mamometer frame are shown in Fig. 3. Two inverted U-shaped frames (U) 
ved as the supporting structure for the models. A horizontal plate (P) 
4S suspended from the U-frames by four slender aluminum rods (r) 4 feet 
ng. The dynamometer plate was restrained from horizontal movement by 
dynamometers designated in Fig. 3 as hg, hh, and N. The first two pre- 
mted the plate from moving in the direction of wave travel, and the third 
pplied horizontal restraint in the direction normal to wave motion. The 
odels were suspended from the aluminum plate by two struts (s) stream- 
1ed in the direction of wave motion. The length of the struts was adjusta- 
e, so that the vertical position of the model could be varied. 

Dynamometers hg and hp consisted of small aluminum cantilevers fas- 
ned to one of the I-beams, and connected to the dynamometer plate by 
eans of small bars in such a manner as to be unaffected by vertical forces. 
mamometer N was fastened to a cross member as shown in Fig. 3. Elec- 
ical resistance gages were cemented on opposite sides of the cantilevers 
d connected to opposite arms of a Wheatstone bridge in a Brush strain am- 
ifier. The dynamometers were calibrated by subjecting them to known 
rizontal forces and recording the corresponding change in resistance by 
eans of an oscillograph. The relationship between the force on the dyna- 
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TESTING ARRANGEMENT 
Fig. 2 


mometers and the force on the model was determined by applying loads at 
various model positions and noting the resulting forces on the dynamometers. 

Vertical forces were measured by means of two dynamometers, vg and 
Vp, Which consited of slender horizontal beams through which the four sup- 
porting rods were connected to the testing frame as shown in Figs. 2 and 3. 
Strain gages were placed on top and bottom of the center of these beams. 
Vertical loads on the model produced bending moments in the beam and 
changes in the resistances of the strain gages which were recorded by means 
of an oscillograph. These dynamometers were so arranged that they were 
not affected by horizontal forces. 


Test Procedure 


An extensive series of preliminary investigations was made before be- 
ginning the principal tests in order to determine the behavior of the waves 
in the testing area. One phase of the preliminary investigation was planned ~ 
to determine the effect of the presence of the model on wave heights in the 
vicinity of the model. The placing of the model in the water caused vari- 
ations in the wave height for some distance upstream. However, a point 11.0) 
feet upstream from the model was found to be free from this effect for all 
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conditions, and was, therefore, selected as the wave gage location. 

The variation in wave height in the direction perpendicular to the wave 
motion was also studied. For certain wave frequencies, standing transverse 
waves were developed which produced a variation in wave height from one 
side of the model to the other of as much as 40 per cent. For such conditions, 
a single measurement of wave height at the center of the model provided a r 
value which differed as much as five per cent from the average effective 
wave height. This phenomenon was caused by the establishment of a resonant 
motion in the channel and tank, and no method was found to eliminate it. 
However, these standing transverse waves did not appear until approximately 
seventeen waves had passed the model following the beginning of a test. Also, 
the first 4 or 5 waves of a group changed in size and length as they progressed 
along the channel from the wave rod to the model. Both difficulties were 
avoided by using only those waves between the 5th and 17th in a series. 

The dynamometers were calibrated by loading the model with known | 
weights. The periodic loading resulting from wave motion was simulated by : 
manual manipulation of the weights. The calibration coefficient for the sum 
of the two principal horizontal dynamometers differed by 2.5 per cent from i 
the beginning until the end of the testing program. For the vertical dyna- 
mometers there was no measurable difference in the calibration coefficient 
from the beginning until the end of the testing program. 

The wave rod was calibrated by raising and lowering the rod known 
amounts and noting the corresponding oscillograph reading. It was necessary 
to conduct the calibration in the wave tank itself because the calibration was : 
affected by the chemical composition of the water or the grounding conditions. 
The calibration curve was nearly linear, and, although it remained virtually 
constant, a new calibration was made each day. 

The desired water depth was maintained by adding water to compensate } 
for evaporation. The water depth was measured at the front of the model. } 
The concrete floor of the tank was relatively smooth but there were undu- 
lations as great as 1/4 inch in the region between the model and the wave 
machine. i 

Prior to starting a test run, the wave machine was adjusted to generate } 
a wave of the desired height and period, and the model was placed at the i] 
desired elevation. During the first run simultaneous measurements were 
made of the wave height at a point 11 feet upstream from the model and of 
either the horizontal or the vertical forces. Such a test required three os- 
cillograph channels because both the horizontal and vertical forces were 
obtained by adding the forces in two dynamometers. Upon completing such a 
run, one dynamometer was switched off and a record was obtained from a 
wave rod located at the side of and opposite the center of the model. The 
wave heights measured by this gage were distorted by the presence of the 
model and were used only to obtain a phase relationship between the wave 
and forces. In the early tests, records were also taken from dynamometer 
N to determine forces normal to the direction of wave motion. This was 
discontinued because these forces were found to be very small. 


The Experimental Program 


The testing program was planned to determine horizontal and vertical 
forces on barges of several shapes for a variety of wave characteristics and 
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MODELS 
FIGURE 4 


or various barge locations. Representative results were sought throughout 
o as to establish trends rather than to provide complete results for all 
ossible combinations of variables. The depth of water (d) was kept constant 
t one foot throughout the tests. This depth, as well as the model dimensions 
nd the wave characteristics, were selected as being typical of representative 
rototype values reduced in a linear scale ration of 1 to 100. 

The basic model was a rectangular parallelepiped 2' 6" long, 10 inches 
ide, and 1-7/8 inches high as shown in Fig. 4. Two modifications of the 
asic model were studied. One modification consisted of the inclusion of a 
lot in one end of the barge to simulate a typical drilling structure, and the 
econd was the rounding of the upstream and downstream sides as shown in 
ig. 4. Tests were also made on a flat plate having the same area as one 
ide of the barge and a thickness of 1/4 inch. All tests were conducted with 
le top edges of the models horizontal and with the long side (2' 6"' x 1-7/8'') 
Irmal to the direction of wave motion. 


For each barge position and wave period, waves of three sizes were pro- 
jected against the model. The smallest wave height (H) in such a series was 
usually about 0.1 foot, the intermediate one approximately 0.2 foot and the 
largest greater than .25 foot. Three wave periods (T) were used in the tests 
on the basic model. These were approximately 1.4 seconds, 1.1 seconds, and 
0.85 second, the corresponding wave lengths (L) being 7.1 feet, 5.2 feet, and 
3.5 feet, respectively. Values of d/L thus varied from 0.14 to 0.29. Actual 
wave periods for individual tests differed slightly, the values given above 
being approximately the averages for the three groups of tests. For the 
period of 1.1 seconds, the model was tested at 10 locations varying from a 
position with the center of the barge 0.17 inches below the water surface to 
a location as near the bottom as possible without touching the bottom. For 
periods of 1.4 seconds and 0.9 second, all tests were conducted with the barge ~ 
at mid-depth. The slotted and rounded barges were tested only for a wave i 
having a period of 1.1 seconds. 


1833-8 HY 6 November, 1958 


Oscillatory Wave Motion | 
Oscillatory waves are so called because the water particles oscillate in j 
nearly closed orbits, there being substantially no forward motion of the fluid. ? 
In relatively deep water (d/L>0.5), the orbits are circular at the water 
surface, becoming slightly elliptical and much smaller as the bottomis ap- 
proached. In shallow water (d/L< 0.5), the particle orbits are elliptical 
throughout the depth. The analysis of the test data requires the use of theo- 
retical relationships which define the form and orbital motion of the waves. . 
Several rather complete summaries of wave theory are available. (9,10,11,12) © 
Equations for the shape of and motion in oscillatory waves were developed 
by Airy.(13) Although his derivation was based on the assumption that wave 
heights are small, his results are sufficiently accurate for many purposes 
even though used to describe waves with appreciable amplitudes. Stokes(14) 
developed equations which include the effect of wave height, but they are ex- 
ceedingly cumbersome. The small refinement provided by them was not: 
justified in the analysis of the test results. Consequently, the relationships 
utilized and presented in the following paragraphs are those credited to Airy. 
The wave length (L), the distance from crest to crest of adjacent waves, 
is related to the other characteristics of a wave by the following expression: © 


2 
24 
L = tanh =. (1 


in which T is the wave period in seconds and d is the depth of water. The 
relationship between wave period, wave velocity (C), and wave length is. 
L = CT 


The approximate equation for the wave profile is 


y - 5 cos oie} 


in which x and y are the horizontal and vertical coordinates of the water 
surface and H is the wave height. Values of x are positive in the direction of 
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wave motion and y is measured positively upward from the still water level 
(s.w.1.) as illustrated in Fig. 5. The expression 27x/L is the phase angle 
in radians which is designated as 6. In degrees, § varies from 0 at a crest 
to 360° at the preceding crest, as shown in Fig. 5. The phase angle can also 
be expressed in terms of time as 27 t/T where t varies from 0 to T during 
a wave cycle. 

The equations for the horizontal and vertical components of the orbital 
velocity are 


= 7H cosh 27 (d + z)/L 


oo sin Ta cos 8 (4) 


and 


~ 7H sinh 27 (d + z)/L 


a i sinh 27rd 


sin 6 (5) 


In these equations, z is the vertical distance from the still water surface to 
the point where the velocity is being determined, upward being taken as posi- 
tive. The corresponding expression for horizontal and vertical accelerations 
can be determined by differentiating Eqs. (4) and (5) partially with respect 

to time. The resulting expressions are 


ou, 2074 cosh 27° (d + z)/L 


-_—_- 1—————— sin © 
ot T sinh 277 d/L (8) 
and 
gu, -27° 74 sinh 277 (da + z)/L (1) 
= ———_—_—————- cos 8 
ot 7? sinh 277° d/L 


Oscillatory Wave Forces 


The forces produced on submerged structures by oscillatory waves are of 
two types: the drag force, resulting from the orbital velocity; and the in- 
ertial force, resulting from the orbital acceleration. The drag force can be 
expressed as follows, 3 

ee eS (8) 
in which Cq is a drag coefficient, P is the density, A is the cross-sectional 
area of the structure in a plane normal to the force, and u is the component 
of the velocity in the direction of the force. The flow pattern referred to is — 
that which would occur if flow were not disturbed by the presence of the 
model. Because of the nature of the velocity variation, the horizontal drag 
force varies from a maximum downstream value under the crest of a wave 
(x/L = 0) toa maximum upstream value under a trough (x/L = 0.50) as 
illustrated in Fig. 5. 

For steady flow, the drag coefficient has been examined for many differen 
shapes and has been found to be a function of the form and of the Reynolds 
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umber .(15,16) However, it was not known whether such values could be 
ised for the periodic unsteady flow produced by oscillatory waves. An ad- 
ditional complication arises for the case of large structures because of the 
variation of the motion within the dimensions of the body. 

The inertial force is expressed by the equation, 


re) 
Fi - Cv SY (9) 


in which C,, is the coefficient of inertial resistance, V is the displaced 
volume and 0u/@t is the component of the acceleration in the direction of the 
force for undisturbed flow. This type of force occurs for two reasons. The 
acceleration of the fluid is caused by a pressure gradient so that the pressure 
on the upstream face of the submerged structure is not the same as that on 
he downstream face. The resulting force can be shown to be equal to the 
oroduct of the acceleration and the mass of the displaced volume. Also, 
yecause of the presence of the body, the accelerations in the flow around the 
20dy are increased, thus requiring an additional force. This force divided 
dy the acceleration of the ambient flow is called the virtual mass, and is 
sonveniently expressed by its ratio to the mass of the fluid displaced by the 
90dy. Thus, C,, is the sum of a term, approximately unity, for the pressure 
zradient force, and the coefficient of virtual mass. 

The variation of the horizontal acceleration, Eq. (6), is such that the 
10rizontal inertial force varies from a maximum downstream value at 
(/L = 1/4 to a maximum upstream value at x/L = 3/4 as illustrated in Fig. 5. 

Values of the virtual mass coefficient have been determined theoretically 
ind experimentally for bodies of various shapes. For irrotational flow with- 
jut separation, values are available for elliptical cylinders and ellipsoids(17) 
ind for infinitely long bodies of rectangular cross-section(18) as shown in 
Fig. 6. The corresponding values obtained experimentally, (19,20) by 
scillations at high frequency and small amplitude, correspond fairly well 
vith the theoretical values. Values of C,, obtained in this manner may not 
ye applicable to flows in which separation plays a predominant role. In fact, 
he development of a wake has been found to cause marked changes in both 
Sq and Cy. 9) Variations in flow over the length of large bodies cause 
urther uncertainties. Consequently, the applicability of theoretical results 
o flows produced by oscillatory waves must be tested experimentally. 

The total force is the sum of the two forces given by Eqs. (8) and (9). 
F=C,A -3~ +, us (10) 
3ecause the two forces are 90° out of phase, the maximum combined down- 
tream force may occur at any value of x from x/L = 0 to x/L = 1/4, depend- 
ng on the relative magnitude of the components. The magnitude and location 
f the maximum force for the case of cylindrical piles has been discussed 
y Morison and others. (2,3 

For large structures, it is convenient to express the horizontal inertial 
orce in terms of the pressures exerted on the body. If the rectangular out- 
ine shown in Fig. 5 is assumed to contain only the undisturbed fluid, and if 
onvective accelerations are neglected, then an application of the equations 
f motion to this volume of fluid in the x-direction yields the following 
xpression 


(Py - Pg) Ay =v 23% (11) 


1833-12 HY 6 November, 19 


ll 2 4 6 8 1.0 2 a& 6 8 
a/b 


Cm FOR TWO DIMENSIONAL FLOW 
PAST RECTANGULAR BODIES 
Fig. 6 


in which p, and pg are the pressures on the two sides of the volume of fluid, — 
Ax is the cross-sectional area of the volume taken normal to the force, and 
Oux/Odt is the average horizontal acceleration of the fluid. Substitution in 
Eq. (9) results in the following expression for the inertial force 


Fy = "Gy, (py o Po) Ax (12) 


From the equation of motion written for the vertical direction, tiie follow- 
ing relationship can be obtained, 


Og ou 4) 
Sz =W + P ee (23 
in which w is the unit weight of the fluid. 
The substitution of the value of Ou,/ Ot from Eq. (7) and integration 
yields 


p = -wz + et? Gosh At SG 227 ae nmin (14 


re sinh 277d/L 


The constant C can be evaluated because of the fact that p = 0 for points on 
the water surface, i.e., for z = y. The expression for the pressure then 
becomes 
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Le w(-z+y)+@ LH cosh 27(d+z)/L - cosh 27(d+y)/L cos 0 (15) 
T sinh 27 d/L 
‘he value of (p; - Pg) to be used in Eq. (12) is found by inserting in Eq. (15) 
ppropriate values of x and y for the two sides of a submerged body. 
For the particular case of a rectangular barge centered at the quarter 
int of a wave, substitution in Eq. (15) yields 


2177 X, 9 (7H cos 277 x,/L 
(Py-Po) mr ME 208 top R= SEO a oe, (16) 


n which K is defined as follows, 


cosh 27(d + z)/L 
cosh 277° d/L (17) 


K = 

nd x; is the value of x at the upstream side of the model as shown in Fig. 5. 
q. (16) indicates that the pressure difference on the two sides of the body 
esults from the difference in the water surface elevation, 4 H, 
AH = H cos 2 7 x1/ L) modified by the term in the brackets which corrects 
yr vertical accelerations. 

The difference in the pressures occurring under a crest and a trough can 
lso be obtained from Eq. (16) by letting x; = 0 andz=-d. Then 


©. 2 PP 2 

yp = K - 2 sinn® TB (18) 
1 which p, is the pressure under a crest and p; the pressure under an 
djacent trough. Folsom(21,22) measured values of the ratio shown on the 
ft side of Eq. (18) as part of a study of a pressure wave gage and compared 
em with K. He found that the measured values of the ratio were between 10 
ad 20 per cent less than K. The correction provided by the additional term 
| the right side of Eq. (18) improves somewhat the correspondence between 
leory and the results of Folsom’s experiments. 


Wave Forces ona Flat Plate 


The results of the tests for horizontal forces on the flat plate are pre- 
snted in Table 1. The wave heights (column 4), measured forces (columns 
and 7), and location of model at the time of maximum force (columns 6 and 
| were determined by averaging values from three or four successive waves. 
eights of individual waves were taken as the average height of a crest above 
jjacent troughs. Separation of the upstream force from the downstream 
ulue was based on the assumption that position of the oscillograph pen cor- 
sponding to zero force did not shift during a test. Although this change 
as small, it may have caused minor inconsistencies in individual values. 

The wave lengths shown in column 3 were computed by means of Eq. (1). 
he measurement of wave lengths was not part of the routine testing program. 
evertheless, a number of wave lengths were measured and the observed 
lues corresponded with the computed values within 5 per cent. 
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In the analysis of the wave forces on a flat plate, both the drag force and 
the inertial force were found to be significant. The total horizontal force is 
determined from Eq. (10) by using the horizontal components of the velocities 
and accelerations. Values of F were computed by inserting into Eq. (10) the 
quantities which apply to this case. Values of Cq and C,, were selected to 
give the best correspondence with the measured phase angles and magnitudes 
of the maximum forces. A value of 3.5 was found for Cg and the value of 
CmV was found to be 1.75 times the volume of the circumscribing cylinder, 
which is 1.75 times the theoretical value(20) for an infinitely long thin flat 
plate. Because these values are large, an effect of the unsteadiness of the 
flow on the shape of the wake is indicated. Values of uy and Ou,/0O t were 
computed from Eqs. (4) and (6), respectively. 

For any particular value of wave height and wave period and for a par- 
ticular model location, F is a function only of the phase angle and Eq. (10) 
can be expressed as follows: 


2 

F = C,' cos" + C' sin @ (19) 

in which 

2-2 2 
C7 “H“ |cosh 277 (d+z)/L 
Ce herm CA — 20 
d a” er? sinh 27°d/L (au) 
and 

27H cosh 27 (d+z)/L (21) 


oe % ce 2 sinh 277d/L 


The phase angle corresponding to the maximum force is obtained by setting 
dF/d@ = 0, with the result shown in Eq. (22). 
(8 1 
m 


2 
2c (22) 


Computed values of 6,4, expressed in terms of x/L are shown in column 9 
of Table 1 and the corresponding values of Fg, Fj, and F are shown in 
columns 10, 11, and 12, respectively. Computed phase angles and forces are 
presented together with measured values in Table 1. These values are also 
compared graphically in Fig. 7, in which computed values are plotted against 
both upstream and downstream measured values. Although there is consider- 
able scatter in the data for the phase angle (Fig. 7a), the trends indicate good 
correspondence between measured and computed values. The correspondence 
between measured and computed forces, shown in Fig. 7b, is good for 

z/d = -0.25 and -0.50. For z/d = -0.75 the measured values exceed the 
computed values by about 30 per cent. 


sin 6 = 
max 


Horizontal Forces on Barge Models 


Results of tests for barge models are presented in Tables 2 and 3. In 
Table 2 are shown results of a number of tests for which the data were 
analyzed in considerable detail. The values of wave height, maximum force, 
and phase angle at the time of the maximum force were determined by 
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x/L Computed 


o2 


x/L Measured Measured Force in Pounds 


(a) (b) 


FLAT PLATE DATA 
Fig. 7 


averaging the results for three or four successive waves. The downstream 
forces were separated from the upstream values and are shown separately 
in Columns 8 and11. The upstream forces were, on the average, 6.6 per 
cent smaller than the downstream forces. These values are subject to the 
same instrumental limitations as are those for the flat plate. 

The difference in water surface elevation on the two sides of the barge, 
aH, at the time of the occurrence of the maximum force were scaled from 
the oscillograph charts. The values are shown for the maximum downstream 
forces in column 9 and for the maximum upstream forces in column 12. 

In Table 3 are shown the results of a number of tests which illustrate the 
effect on the force of the position of the barge in the vertical. These tests 
were conducted on the rectangular model for three different wave heights 
with the wave period held constant at approximately 1.1 seconds. The locatio 
of the center of the barge (-z) was varied from 0.08 feet to 0.92 feet below 
the water surface. In the latter position, the barge was very close to, but 
not touching, the bottom. Three different wave heights were used as shown 


WAVE FORCES 1833-17 


\SCE 


6hI°O zZ9°T | 9tt°O ee | 2L-Ot Lor 0 6r'% 6t’o | £60" | 9LGO |e See oiet S270 it 
080°0 zt | S20°0 zs°t | 8z°0 | $90°0 Zp't gz°o | 6TE0” | 891°0'| 92'S | TIT Psz70- peiiots 
PET'O LUT 6Z°ci lr P26 o9so° | s6z°o | 9z°s | It°T | S2°0- % 
¢so°O SPT ZU°t | FL°0 O20" || LZT"0 |-9e, Gh > ELoran Gece F 
Zito .Se"s 16°T | 92°0 ozv0' | 2¥Z°0 | 9%°S | TI*T | sZ°0- FA 
zOT°O eS°T zs°z | 18°0 9Zv0° | p2z°0 | 92'S | TI°T | $z°0- pepunoy 


910°0 9F°T €90°0 96°T 9L°0 
0v0°O 6S°T 0r0'°O To°T SL°O 
LET O GS TL OTTt’O Se°s cL°O 
T¥O0°O cS I €r0 0 LL’°0O LL‘O 
9€0°0 <éL’T se0°0O v8°0 LL’O 
G00 L9°T vEO’O 89°0 8L°0 
ceé0°0 ast Z£0°0O 28 °0 9L°0 
of0°O 9L°T 0£0°0 c3°0 08°0 
O€T°O 6S°T €It’oO v9°G LL’O 
6eT°O PET 60T’°O LE°S 9L°0 
621°O Ter 611°0 ZE°S LL‘O 
621°O G6E°T €cT°O 8Z°E LL’O 
2OL:O OL LT 160°0 €2°S €L°0 
COT, 0 0G" 160°0 v0°? LL’O 
S60°0 6%°T ¥60°0 Il’? 6L°0 
TOTO 6S'T ¥60°0 89°S LL‘0 


96z0° | 21z°0 | ZT°Z | Th°T | 0S’O- a 
ESTO” PCIE OLLE L | 214. 0S:0— a 
GLSO° | O1z°0 | 99° Jezs° | OSs*o- re 
OLZTO’ | 190°0 | ss°e |z98° | OS‘O- Pr 
€vto’ | yL0'0 | 61'S | OTT | Z6°0- ; 
6ET0° | €20°0 | 9z°S | TIT | Sz°0- a 
6EI0’ | €20°0 | 92's | TI°T | oS°o- ns 
PZ10° | ¥90°0 | 9z°S | TI°T | SZ‘o- is 
SzGO""| LS" O41 SES. 1 Sto ce 0= ‘“ 
sIso* | Z2Z°0']'92°S |. TT°T} Sz°0- ‘i 
61S0° | 692°0 | 61T°S | OT°T | OS*0- uo 
LIso" | 222°.) 92°S | 1°11] S$3°0- Ps 
FOrO’ | 20Z°0 | Z's | 6O°T | 26°0- ¥ 
ZIv0° | 60Z°0 | 90°S | 80°T | Sz°0- 
s6eo° | z0z°0 | Z2t°s | 60°T | OS‘o0- “ 
ZOvO’ | 902°0 | I'S | 6O'T | SZ°O- | reTNSuezoOYy 


400g ye0q 44 998 TSPpon 
peynd 8% HY 1/x SB 4 p/z jo odfy 

-ul09 

HY eoi10g weer3sdyn e@o010g weat}zsuMog ‘ 

ot. eteel) etl [° tt. fb ot 8 eee La ee a ee ee ee T 


STAGON ADYVA NO SHOUOA TVLNOZIYOH 


2 ATTEAVL 


1833-18 HY 6 November, 1958 
TABLE 3 


HORIZONTAL FORCES ON RECTANGULAR BARGE 


For T = 1.1 Seconds 
Average Maximum Force in Pounds 
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in Table 3. Individual wave heights varied by less than five per cent from 
the average values for each group. The forces shown in Table 3 are the 
average of the upstream and downstream forces. From the results shown 
in Table 2, it can be estimated that the downstream forces were about 3 per 
cent larger than the average values shown. 

The maximum horizontal forces were created when the quarter points of 
a wave (x/L = 1/4 or 3/4) were approximately at the center of the barge as" 
shown in Table 2, columns 7 and 10. Hence, the resistance caused by the 
velocity was negligible, and the maximum horizontal force was of the inertial - 
type expressed by Eq. (9) or (12). Computed values of AH are shown in 
Table 2, column 14, together with the measured values for the rectangular 
and slotted barges (columns 9 and 12). The average of the computed values 
is 2 per cent higher than the average of the measured values for the down- 
stream forces and 11 per cent higher than the measured upstream values. 

Values of C,, were computed for each test from Eqs. (12), (16), and (17) 
and are shown in column 13. The value of x; is (L/4 - L'/2) in which L' is 
the width of the barge. For the case of the rounded barge an average width 
of 0.800 feet was used instead of the maximum width of 0.833 feet. The 
method of determining the values of C,, is illustrated by the following com- ~ 
putation for the test recorded at the top of Table 2. The value of K, from 
Eq. (17), is 0.784, and 4H = 0.101 (from column 14). The value of (1 - pe), 
from Eq. (16) is 4. 94 psf. Finally, Cy, is found from Eq. (12). | 
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Ms Fy 3.07 
(P] - Pg) Ax 4.94 x% 0.39 ~ 1-59 


Cy 


For the slotted barge the displaced mass is used as a reference. In the 
omputation of Cy,, the forces on the inner faces of the slots (p,' and pg) 
nust be considered. The expression for the force is as follows: 


F = Ch [@ - Pg) Ax - (py' - pg") x'| (23) 


n which A," is the area of the inner faces of the slots. 

The test data are summarized graphically in Fig. 8. The coefficient of 
nertial resistance for the rectangular barge is plotted against the ratio of 
vave height to wave length for three values of z/d with T constant at 1.1 
econds in Fig. 8a. The values of C;, decrease significantly as the waves 
wecOme steeper. The curves indicate the approximate trends. 

In Fig. 8b, values of C), for the rectangular barge are plotted against the 
elative wave height for three wave periods, the barge position (z/d) being 
0.50. In Fig. 8c, Cy, is plotted against relative wave height for three differ- 
nt barge shapes, for constant values of both z/d and T. 

The wave force is plotted against z for three different wave heights in 
"ig. 9. These are the data from Table 3 with all force values increased by 
hree per cent to transform them from average to downstream forces in ac- 
ordance with the results presented in Table 2. The theoretical curves for 
he three wave heights were determined by means of Eqs. (12) and (16). The 
everal constant values of C,, used in determining these curves were selected 
0 give a good average fit. 

The tabulated data and the graphs indicate a general conformity between 
heory and test results. However, forces produced by the largest waves were 
ess than would be expected from the results for the smaller ones. This is 
llustrated in Fig. 8 by the decrease in C,, with increasing wave steepness 
nd also in Fig. 9 in which the value of C,,, which produces the best fit for 
he largest waves, is much less than that for smaller waves. Some deviation 
rom the theoretical curves is also indicated in Fig. 9 for tests with the 
20del very near the bottom or very near the water surface. 

The forces on the rounded barge were approximately 14 per cent less than 
10Sse on the rectangular barge. The reasons for the lower force on the 
ounded barge can be explained with reference to Eq. (9), which shows the 
3rce to be proportional to the product of Cy, and the displaced mass. Four 
er cent of the reduction is the result of the smaller displaced mass of the 
ounded barge. The remaining reduction in force was caused by the in- 
reased streamlining with the resulting lower values of Cy, as illustrated 
1 Fig. 8c. 

Although the displaced mass of the slotted barge is 26 per cent less than 
iat of the rectangular barge, the horizontal forces were only approximately 
9 per cent less. Thus, the reduction in mass is partially offset by the ad- 
itional disturbance attributable to the odd configuration of the slotted barge 
nd the resulting increase in Cj, illustrated in Fig. 8c. 

Values of C,, for the rectangular model varied from 1.31 to 1.76. The 
alue obtained from theoretical considerations of two-dimensional flow for 
lis shape is 1.33 as shown in Fig. 6. Differences are expected because of 
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ap z/d =-25 

z/d*~50 + T#Q.87 Sec. 
2/d =-.75 1.8 ° T= 1.11 See. 
x Ts 1.41 Sec. 


1.6 1.6 x 
Cm Cai e+ ‘ 
x 
1.4 1.4 
Rect. Barge Rect.Barge = 
T=L.ISec. z/d=.50 


+ Slotted Barge 
e Rectangular Borge 
* Rounded Barge 


COEFFICIENTS OF 


1.6 INERTIAL RESISTANCE 
For Horizontal Forces 
C 
. on Barge Models 
1.4 Fig. 8 
T=1.1Sec. 
z/d= .75 
12 
Oo 02 04 .06 
H/L 
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effects of periodic wake formation and because the effective dimension of | 
the barge is an appreciable fraction of the wave length. 


Vertical Forces on Barge Models 


Test results for the determination of vertical forces are shown in Table 
In this case, the dynamometers were not sensitive enough to warrant the 
separation of upward forces from downward forces. However, the records 
do indicate that upward and downward forces were nearly equal for inter- 
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X He 0.07 Ft. 
+ He 0.21 Ft. 
e H= 0.27 Ft. 


Cm=!.47 


Cm = t.52 


Rect.Barge 
T=1.1 Sec. 


Force in Pounds 


HORIZONTAL FORCES 


on Barge Models 
Fig. 9 


nediate values of z, and that the upward force was greater than the downward 
orce if the barge was located near the water surface or near the bottom. 

‘he values of vertical forces, shown in Table 4, are the averages of the up- 
7ard and downward forces. The location of the barge with respect to a wave 
rest at the time of occurrence of the maximum force is shown for a number 
f tests on each of the barges in columns 6 and 7. Wave lengths, shown in 
olumn 4, were computed from Eq. (1). 

The maximum downward vertical forces on the barge occurred if the crest 
as directly above the center of the barge (x/L = 0) as shown in column 6 of 
‘able 4. The maximum upward force occurred in most cases if the trough 
as nearly over the center of the barge (x/L = 1/2), as shown in column 7. 
hus, the vertical forces are also primarily inertial. Vertical forces cn be 
nalyzed by means of Eqs. (7) and (9). 

Values of C,,, obtained from Eqs. (7) and (9) and enumerated in Table 4, 
an be compared with the classical values obtained from theory shown in 
‘ig. 6.(18) Although theoretical results for rectangular bodies are available 
nly if the length is infinite (two dimensional flow), experiments which have 
een conducted (19,20) indicate that the end effects for the lengths used are 
eductions in the virtual mass coefficient of five to ten per cent. For the 
roportions of the barge used in these tests, the theoretical value for two 
imensional flow is 6.1, which is nearly the same as the average of the ob- 
erved values. The average observed value for the rounded barge (referred 
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Type of 
Model 


Rectangular 


Rounded 


"W 


-0.25 
-0.50 
-0.75 
-0.92 
-0.50 
-0.50 
-0.50 
-0.50 
-0.50 
-0.17 
-0.33 
-0.42 
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TABLE 4 


H 
tt] Feet 


0.215 
0.198 
0.198 
0.196 
0.306 
0.252 
0.066 
0.210 
0.110 
0.207 
0.187 
0.193 
0.199 
0.194 
0.205 
0.229 
0.230 
0.228 
0.147 
0.230 
0.135 
0.141 
0.148 
0.147 
0.152 
0.321 
0.322 
0.145 
0.151 
0.289 
0.316 
0.317 
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Downward| Upward 
Force Force 


November, 1958 


Av .Max 
Force 
Pounds Cm 


0.26 | 7.31 
0.42 | 4.54 
0.50 | 2.15 
0.49 | 1.33 
0.27 | 5.63 
0.39 | 6.48 
0.46 | 2.06 
0.41 | 3.45 
0.45 | 1.89 
5.85 
6.67 
5.69 
3.96 
3.14 
1.72 
0.281 97.44 
0.37%) 3.72 
0.42 | 2.38 
0.53 | 1.80 
6.00 
5.96 
5.26 
2.60 
2.22 
0.45 | 1.50 
0,327) 7a 40 
0.34 | 1.49 
3.51 
2.29 
6.25 
4y 14 
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) the somewhat smaller displaced mass of this shape) was 5.9, which is not 
ignificantly different from that of the rectangular barge. In comparison, 
r an elliptical section with the same proportions, an extreme in rounding, 
theoretical value of 6.33 is obtained for two dimensional flow. Finally, an 
verage value of Cy, of 3.8 was obtained from the tests on the slotted barge 
sing its displaced mass as a reference. Although no theoretical value can 
2 obtained for this shape, the weighted average of the theoretical values for 
vo dimensional flow, for the unslotted portion (6.1), and for the legs (2.5), 
; 4.9. In view of geometrical differences and possible effects of separation, 
ie theoretical and observed values of Cm correspond surprisingly well. 
The data are also shown graphically in Fig. 10. The forces on the rec- 
ngular and rounded models were plotted for various values of z in Fig. 10a. 
he group for the rectangular barge for a wave period of 1.1 seconds com- 
rises the values shown in Table 4 for which the wave heights were approxi- 
ately 0.20 feet. The group for the rounded barge is made up of the values 
om Table 4 having wave heights of approximately 0.23 feet. For this group, 
e forces were reduced by the ratio 0.20/0.23 before plotting, because, for 
small range in wave height, the forces vary linearly with the wave height 
_accordance with the theory. The curves shown in Fig. 10a were deter- 
ined from Eq. (7) and (9), using a C,, of 5.5. This value was selected to 
ve the best fit for this group of data. In Fig. 10b are shown groups of 
ints for the rounded and slotted models. These values are the ones shown 
Table 4 for wave heights of approximately 0.15 feet. The curves in Fig. 
ib were determined from Eqs. (7) and (9) using the values of C,, shown. 
Good agreement was found between analytical and test values of the verti- 
il forces on the three models used in this study. Values of C,, correspond 
lite well with the theoretical values for two dimensional flow except for the 
ises in which the models were very near the surface or the bottom. In 
ese cases, there was some indication that the upward forces were larger 
an the downward ones so that the values reported here (which are the 
rerage of the upward and downward forces) should be increased somewhat 
r design purposes. Further studies are required to determine the exact 
ture of the forces in these locations. 


SUMMARY 


Forces exerted on submerged bodies by waves can be reliably predicted 
;a combination of the inertial and drag forces if the appropriate coeffi- 
ents can be determined. The tests on the models of barge-like structures 
owed that the inertial forces were predominant so that only the coeffi- 
ent of inertial resistance was needed to predict the maximum force. For 
e horizontal forces, the measured coefficient was approximately 1.5 
pared to 1.33 as estimated from theory for flow without separation. For 
e vertical forces, the measured and estimated values of the coefficient 
sre approximately six for the unslotted barge. Some systematic variation 
the coefficient was found with relative wave height and with the vertical 
sition of the model. Measured values for the slotted and rounded barges 
sre also in good correspondence with theoretical values. Because the forces 
sre primarily inertial, the rounding of the edges of the barge did not cause 
rge changes in the forces. 

The maximum horizontal forces on a vertical plate placed normal to the 
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—+— Rectangular Model 
—~+*-- Rounded Model 
—_— mised Model 


=8 T= 1.1Sec. 
H= 0.2 Ff. 
Cm = 5.5 


0 2 4 6 8 ° 2 4 6 8 
Force in Pounds Force in Pounds 
(a) (b) 


VERTICAL FORCES 


on Barge Models 
Fig. 10 


direction of wave propagation were predictable if a drag coefficient of 3.5 
and a virtual mass equal to 1.75 times that of the circumscribing cylinder 
were used. Both of these values are large because the growth of the wake in 
unsteady flow increases the drag coefficient and because the presence of an 
wake increases the virtual mass. The phases of the wave which produced 
maximum forces corresponded satisfactorily with the predicted values. 
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SNOWMELT RUNOFF 


J. Harold Zoller,! M. ASCE and Arno T. Lenz,2 M. ASCE 
(Proc. Paper 1834) 


SYNOPSIS 


The primary factors relating to the melting of snow were evaluated for the 
nowmelt periods of 1938 through 1952 for the Big Eau Pleine River basin in 
fisconsin. Convection and condensation melt potentials were computed from 
quations developed by previous investigators. A nomograph was devised for 
omputing radiation snowmelt potential from insolation, as related to sun po- 
ition determined by date, modified by cloud cover and albedo. The as- 
umption was made that the runoff for any day would be a function of the sum- 
ation of melt potentials for that day. A snowmelt unit hydrograph was de- 
eloped and relative daily hydrographs were computed for each day. These 
elative hydrographs were then combined by a trial and error procedure to 
eproduce the recorded snowmelt hydrograph for each year. From this 
rocess the daily runoff values which contributed to the total hydrograph were 
Imputed. These daily runoff values were compared with the total melt po- 
sntials for the same days, and the correlation noted was expressed by an 
juation. 


INTRODUCTION 


Spring runoff caused by the melting of the winter’s accumulation of snow 
iten provides an opportunity for the annual storage of a major part of the 
sgional water resources. Reservoirs may be drawn down to produce water 
ywer during the high demand winter months with certain knowledge that melt- 
i of the existing snowpack will replenish the reservoir supply in the spring. 


ote: Discussion open until April 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1834 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. HY 6, November, 1958. 

Prof. of Civ. Eng., Univ. of New Hampshire, Durham, N. H. 
Prof. of Civ. Eng., Univ. of Wisconsin, Madison, Wisc. 
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Research has revealed techniques by means of which the quantity of water 
which will become available in the snowmelt runoff of any basin may be esti- 
mated with considerable accuracy if certain minimum data are acquired. 

The pattern of key research in the past has been to concentrate in small 
control areas the instrumentation necessary to measure every conceivable 
quantity suspected of being an agent in the melt phenomenon. This approach 
has yielded, and is continuing to yield, information of a most vital nature. 
However, to be of any practical value, this information must also be applicabl 
to large basins where average instrumentation exists for the recording of ‘ 
pertinent data. If such applicability can be universally demonstrated, then 
ultimate goal of the researchers will be attained. 

Most of the concentrated research in the field of snowmelt has been con- — 
fined to the mountainous areas of western United States, where much of the 
annual precipitation occurs in the form of snow. Here, snow may accumulate 
to great depths in the winter months and its melting season may extend from 
March to July or even later at the higher elevations. In these areas melting 
is almost always confined to the daylight hours and occurs chiefly as a result 
of insolation, air temperature and air turbulence produced by wind and to- 
pography. 

In contrast, the snowpack in central Wisconsin is seldom over two feet in 
thickness at the time of initial melting and is commonly nearer one foot in 
thickness. It lies largely in open fields on gently rolling topography. Its me 
ing season which also begins in March, never extends beyond April and has a 
continuous duration of from four to twenty days. During periods of high melt — 
potential, melting may occur throughout the night as well as during the day- — 
light hours, and the melting achieved by condensation of moisture on the cold 
snow surface may also be a major factor. Thus the problem of forecasting 
the rate of runoff from melting snow may be recognized as more difficult in 
many ways under the conditions imposed by Wisconsin’s climate than would t 
the case for many regions farther west. 


Melt Potential Theory 


The primary factors relating to the melting of snow are: Temperature of 
the air, temperature of the soil, wind velocity and air turbulence, humidity, 
rainfall and cloud cover. These factors are measurable, within certain limit 
and from their recorded values the potential melting capacities attributable t¢ 
convection, condensation, radiation, conduction and rainfall may be deduced. 
In addition, certain basic characteristics such as elevation, aspect, topograph 
cover, latitude and soil conditions are influential in the final melting capaci ty 
of a given melt potential. ; 

A study was made of the melt potential—runoff characteristics of the Big 
Eau Pleine River Basin of Central Wisconsin. In this study the assumption 
was made that runoff in the snowmelt season was a direct function of the hea’ 
exchange relations between the snowpack and its environment. Accordingly, 
the total melt potential was determined for every six-hour interval in the 
snowmelt periods of the years 1938 to 1952. This total melt potential for any 
interval was considered to be the algebraic sum of the radiation melt po- 
tential, the convection melt potential and the condensation melt potential for 
that interval. It was further assumed that the runoff which occurred for any 
day would occur in a distribution indicated by the variation in melt potential 
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¢ that day. For example, if the total melt potential on a given day for the 
<-hour interval beginning at 9 a.m. were 1.5 inches and the total melt po- 
itial for the interval beginning at 3 p.m. were 0.5 inches, all other melt po- 
itials for the day being negative, then it was assumed that the runoff hydro- 
aph for that day’s melt was a combination of the 9 a.m. and the 3 p.m. hydro- 
aphs in the ratio of 3 to 1 respectively. This procedure involved the de- 
rmination of a snowmelt unit hydrograph and the determination of average 
ily discharges due to periods of melting beginning at different hours in the 
y. Negative melt potential summations were ignored because their effect is 
inimized by the low conductivity of the snowpack and because the time du- 
tion of negative melt potentials was short. Finally, it was reasoned that, if 
= theory and assumptions listed above were correct, a correlation would be 
ind between the daily melt potential and the daily runoff for the various 

ars being investigated. 


mnvection Melt Potential 


The convection melt potential is the quantity of water in inches which theo- 
tically would be melted from snow of 100 per cent quality in a given time 
-erval by the convective transfer of heat. The factors which are of im- 
rtance in determining convective heat transfer are air temperature and 
rbulence, as indicated by wind velocity. Because preliminary studies with 
ervals of time equal to one day proved ineffective, the unit of time finally 
lected for investigation was six hours. Air temperatures and wind velocities 
re recorded at the Wausau CAA airport at 12:30 a.m., 6:30 a.m., 12:30 
m. and 6:30 p.m. for the period of study. The 9 a.m. to 3 p.m. six-hour 
riod is generally the most important snowmelt period of the day and the 
:30 p.m. temperature and wind velocity data may be considered fairly repre- 
ntative of that period. Thus, the four quarter-day periods were selected to 
gin at 3 a.m., 9a.m., 3 p.m., and 9 p.m. It was assumed in all cases that 
2 wind velocity and temperature readings were average values for the six- 
ur periods in which they occurred. 

Although the Wausau CAA Airport is located completely outside the upper 
¢ Eau Pleine Basin, approximately 30 miles east of the center of the basin, 
> Wausau data are considered applicable. Temperatures at Wausau and 
sdford, 20 miles northwest of the center of the basin, exhibit a very good 
rrelation for the snowmelt period. 

The convection melt potential was determined from Wilson’s (1941) 
uation: 


M = KV(T-32° 


here M is the depth of water in inches melted from snow in six hours, V is 
> wind velocity in mph, T is the dry-bulb temperature in degrees F and K 

a constant involving the latent heat of ice, exposure of instruments, air 
nsity, conversion of units, and certain considerations involved in the theory 
turbulence.” A value of K = 0.002 was used. This value of K corresponds 
ysely to Light’s (1941) value of K = 0.00176 as computed for the physical 
nditions imposed by the Big Eau Pleine Basin. The convection melt po- 
itial data for the snowmelt period of the year 1943 are tabulated in Table I. 


The contensation melt potential may be defined as the amount of water in 


ndensation Melt Potential 
shes — from snow, plus condensate, as condensation occurs on the 
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CONVECTION AND CONDENSATION MELT POTENTIALS - 19))3 
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irface of the snow. The factors which determine the condensation melt po- 
ntial are the dew point and the wind velocity, as related by Wilson’s (1941) 
juation: 


M = 2 V (e-6.11) 


here M is the depth of water melted plus condensate per six hours in inches, 
is a constant involving the latent heat of ice, exposure of instruments and 
nversion units, and is equal to 0.0032 for basins with elevations less than 
100 feet, V is wind velocity in miles per hour, e is the vapor pressure of the 
3 in millibars, and 6.11 is the saturation vapor pressure in millibars at 

°F. Records of dew point and wind velocity at Wausau were again con- 
eo applicable to the Big Eau Pleine Basin. The condensation melt po- 
mtial data for the snowmelt period of 1943 are tabulated in Table I. 


adiation Melt Potential 


Adequate data for the direct determination of incident solar radiation for 
.e Big Eau Pleine Basin were not available. Thus it became necessary to 
construct the insolation characteristics from the six-hourly cloud cover 
ita, recorded at the Wausau CAA Airport. Since records of both cloud cover 
id incident solar radiation were available for Madison, Wisconsin, a corre- 
tion between these factors was established for use with the Wausau data. 

Cloud cover is an estimated quantity, reported in tenths by a ground ob- 
srver. A cloud cover of zero represents a cloudless sky, 10 represents a 
ymplete overcast and intermediate values represent intermediate conditions. 
must be emphasized however, that although two days may be reported 
yrrectly as having 10-tenths cloud cover, the insolation reaching the earth 
r those two days may be vastly different since cloud canopies differ in their 
ipacity to filter out solar radiation. A further indication of the inaccuracy 
volved in the use of cloud cover data in determining incident radiant solar 
1ergy is illustrated by several periods in which the noon cloud cover report- 
1 was 10-tenths and the duration of sunshine reported for the hours preceding 
id following the noon cloud cover observation were each 60 minutes. In other 
ords, although the sky was completely overcast, the clouds were of such a 
ture that the sunshine was visibly able to penetrate them continuously. 
bviously, the insolation corresponding to 10-tenths cloud cover on such a day 
ill not correlate properly with the insolation corresponding to 10-tenths 
oud cover on a day when heavy clouds are extant. 

The Madison insolation records for several 15-day periods were checked 
r the 19-year period 1931-1949 inclusive, and the highest insolation value 
r each daylight hour within each interval was determined. From these 
lues were determined the maximum composite insolations which are equal 
the daily insolation in langleys which would have occurred at Madison if the 
ghest insolation values for each hour within a given period had all occurred 
1 the same day. These values were plotted against time to give the Maximum 
ymposite curve shown on Fig. 1. This curve represents the daily insolation 
immation which is approximately the highest insolation that could possibly 
cur at Madison for each day represented. The average insolation experi- 
iced in Madison is seen from Fig. 1 to be a function of this maximum com- 
site insolation. 

The latitude correction to adapt the Madison data to the Big Eau Pleine 
isin is 0.97. This correction was considered insignificant in comparison 
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th uncertainties involved in other assumptions and was not used. Its neglect 
sults in a constant error for which compensation may be made in the se- 
‘tion of the proper albedo value. 

Insolation values from Madison were converted to percentages of the maxi- 
im composite insolation by dividing the insolation value recorded by the 
iximum composite insolation for that day as determined from Fig. 1. These 
ta were then plotted against various combinations of cloud cover data. It 
S found that a significant correlation existed between the per cent of maxi- 
im composite insolation per half day and the sum of either the 6:30 a.m. 

d 12:30 p.m. cloud cover values or the 12:30 p.m. and the 6:30 p.m. cloud 
ver values, expressed in tenths. Table II summarizes the results of this 
rrelation study. It is significant to note that errors resulting from approxi- 
iting the insolation from cloud cover data are later minimized materially 
corrections for albedo and back radiation. It is also emphasized that the 
eatest accuracy of insolation determination occurs during periods of mini- 
im cloud cover, or during periods when maximum melting from insolation 
uld occur. 

Thus, from the cloud cover data at Wausau for any day the per cent of the 
iximum composite insolation for that day may be determined from Table IL. 
r example, on March 30, 1943, the 6:30 a.m., 12:30 p.m. and 6:30 p.m. 
gud covers were 0.4, 0.5 and 0.8, respectively. The sum of the 6:30 a.m. 

d the 12:30 p.m. cloud cover values is 0.9 and the sum of the 12:30 p.m. and 


TABLE II 


INSOLATION-CLOUD COVER RELATIONS 


Sits, RG Sei, ieee 
tothe seo beeee inet adie ian 
Lptheme tad per half day 

P dh 10 3h 

1 L3 rl 33 

; 43 12 31 

; Me a 29 

: Ls 27 

5 ka 1s fs 

é ho 16 23 

u 38 17 1 

P 37 18 19 

3 36 19 a 


5 
w 
Fag 
Oo 
oO 
et 


1834-8 HY 6 November, 195 


the 6:30 p.m. cloud cover values is 1.3. From Table II the morning insolatiol 

was 36 per cent of maximum for the half day and the afternoon insolation wa } 
29 per cent of maximum, so the daily total was 65 per cent of the maximum 
composite insolation for that day. From Fig. 1 the maximum composite inso= 
lation (March 30) was 658 langleys. Thus, the computed insolation for March 
30, 1943 was 428 langleys. 

The radiation components effective in producing snowmelt potential may be 
considered to follow the pattern shown in Fig. 2. The shortwave incident radi 
ation (428 langleys on March 30) consists of the direct and diffuse solar radi= 
ation. The net longwave emissive radiation consists of back-radiation from 
snow, partially balanced by incident diffuse longwave sky radiation. This ne 
longwave radiation has been shown to be emissive (Fig. 2b) and essentially 
constant at 0.10 langleys per minute. During hours of darkness when there ig 
no insolation the energy loss has been termed non-concurrent radiation. Dur 
ing daylight hours this emissive (concurrent) radiation must be deducted from 
the incident radiation to obtain the net radiation. It is evident that the duratio 
and amount of the concurrent and non-concurrent radiations are functions of — 
the total daily insolation. Thus, the insolation effective in producing melt for 
any day is the computed insolation for that day minus the concurrent back 
radiation. 

A variable amount of the incident solar radiation is reflected from the © 
snow surface and is thus unavailable for melting. The percentage reflected, 
called albedo, to be used in reducing the effective insolation is a function of 
the density, structure, moisture content, depth and purity of the snowpack. 
Because sufficient records were not available to permit a proper evaluation 
of albedo variation, a constant albedo value of 60 per cent was assumed for 
the entire snowmelt period. For an albedo of 60 per cent, the computed inso- 
lation for March 30, 1943, of 428 langleys would become 171 langleys. ] 

The net concurrent radiation is a function of the envelope of effective radi 
ation and has been determined from the authors’ emperical equation: 


r = 38.5 + 0.0633 p = 49 langleys 


in which: ris the net concurrent radiation in langleys per day and p is the - 
net insolation as reduced for albedo in langleys per day. Subtracting the net 
concurrent radiation of 49 langleys would leave 122 langleys, which would be 
sufficient to produce 122-+(80 x 2.54) = 0.60 inches of melt. During the sam 
day, the net non-concurrent radiation was (60 x 24 x .1) - 49 = 95 langleys, | 
which is the amount of heat liberated in freezing 0.47 inches of water at 32° 
The insolation nomograph of Fig. 3 provides an easy means of solving for 
the radiation melt potential. To solve for the melt potential using Fig. 3 it ii 
necessary to know the date, the per cent of total insolation and the albedo. / 
straight edge is placed on the nomograph from the date to the per cent of tot 
insolation and the intersection with the net insolation line is noted. A hori- 
zontal line through this net insolation will determine the melt for various 
albedo values and the amount of water frozen by non-concurrent radiation. 
Utilizing the data of March 30, 1943, a straight edge placed from March 30 t 
65 per cent indicates a net insolation of 428 langleys per day. A horizontal 
line through this net insolation value intersects the 60 per cent albedo line < 
a melt of 0.60 inches per day and indicates that 0.47 inches of water would k 
frozen by non-concurrent emissive longwave radiation. The radiation melt 
potential calculations for the snowmelt period of the year 1943 are tabulat 
in Table II. 
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The radiation melt potential values thus determined could not be utilized 
directly in the snowmelt calculations however, because they had to be assimi- 
lated with other melt potential factors which were determined for one-quarte 
day intervals. From a study of daily insolation graphs it was decided arbi- 
trarily to consider that two-thirds of the daily positive radiation melt potent 
occurred in the interval from 9 a.m. to 3 p.m. and one-sixth occurred in each 
of the intervals 3 a.m. to 9 a.m. and 3 p.m. to 9 p.m. It was further conclud 
that the negative melt potential should be equivalent to 36 langleys or 0.18 
inches of melt for the period from 9 p.m. to 3 a.m. and that the remainder of 
the negative melt potential should. be equally divided between the 3 a.m. to 
9 a.m. period and the 3 p.m. to 9 p.m. period. Using this scheme of division, 
the melt potential for March 30, 1943 would be subdivided as indicated in 
Table IV, according to which two-thirds of the .60 inches of melt potential, or 
0.40 inches, is assigned to period 2 and the remaining 0.20 inches is equally 
divided between periods 1 and 3. A total negative melt potential of 0.47 inch 
is divided so that 0.18 inches is applied to period 4 and-the remainder of 0.29 
inches is equally divided between period 1 and period 3. 


Total Melt Potential 


The total melt potential for any time interval is merely the algebraic sum 
of the individual melt potential factors. Table V gives the total melt poten 
as computed for the snowmelt season of the year 1943. 


Snowmelt Runoff Theory 


The streamflow hydrograph resulting from melting snow must reflect the 
variation in total melt potential. The hydrograph is also a function of a snow 
melt unit hydrograph because, like runoff from rainfall, once an increment o 
flow is started on its path it proceeds until it has passed the measuring 
station, unless extremely severe freezing conditions follow melting. 

By preparing a “relative hydrograph” for each day having a positive melt 
potential and making these relative hydrographs reflect the daily variation in 
melt potential as well as the snowmelt unit hydrograph pattern it is possible - 
to divide the daily runoff into the increments which have been contributed 
from snowmelt potential factors existing on previous days. 

To illustrate this procedure in detail numerous examples have been used 
in the following discussion, 


TABLE IV 
RADIATION MELT POTENTIAL FOR MARCH 30, 193 


Period Positive Melt Negative Melt Net Melt . 

Potential Potential Potential 

(1) 3am to 9am el0 ine ~ol ine =e0k ins 
(2) 9am to 3pm Fy Ie) 0 Pyle) 
(3) 3pm to 9pm e10 -olh =.0) 


(4) 9pm to 3am 0 -218 =.18 
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TABLE V 
TOTAL MELT POTENTIAL ~ 19),3 
ite Period Melt Potential in Hundredths of an Inch a Rain in 
Convection Condensation Radiation Positive Inches 
Total 
‘ch 
3 1 6 =-2 -6 
2 28 -10 33 51 
3 22 =) 3 15 
iF 2 -2 -18 
ly i 2 =2 =3 
2 16 ) hh 60 
3 18 ) =a 15 
k 0 ) -18 
5 i -) -1 0 
2 22 5 53 80 
3 30 6 (0) 36 
4 1, 5 -18 al 
6 1 8 -12 0) 
2 18 -18 53 53 
3 8 -18 fe) 
h 22 -16 -18 
7 1 -20 -7 5 
2 “1; -13 38 n 
3 6 -9 5 
h -10 -7 -18 
8 1 =32 =12 -) 
2 8 -7 Tak 26 
3 ~6 -9 a) 
IF =10 -7 -18 
9 a -12 -9 -2 
2 22 5 148 65 
3 18 3 =2 
h 36 10 -18 28 
(0) 1 16 7 =) 19 
2 70 19 ho 128 
3 58 18 =) Te 
Res 58 20 =18 60 
al 1 6 3 -11 
2 66 9 19 9h 
3 30 3 -11 22 
h Uy 6 -18 203 
il 
1 1 2 -1 “i 
2 2 -1 43 hy 
3 8 -7 “4 
h -10 -11 -18 
2 al =26 -1 0 208 
2 =) -18 53 31 
5) -) -11 0 
k -) =2 -18 
3 1 -12 5 - 
; 2 10 -15 39 3h 
3 18 ~10 - 4 
k 10 z -18 
h at 0 fe) 
2 98 9 52 159 
s} 90 =1); (e) 76 
h 28 -12 -18 2 
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Snowmelt Unit Hydrograph 


A trial snowmelt unit hydrograph was determined by isolating the runoff 
from one day’s melt in a manner similar to that used in determining a rainfa 
unit hydrograph. Base flow was removed from consideration by utilizing the 
fall recession curve of Fuhriman (1952). In addition, the effect of the duratic 
of the melt potential producing runoff was considered. The resulting snowme 
unit hydrograph, shown on Fig. 4a, is in close agreement with the rainfall ur 
hydrograph (Beck, 1946), leading to the conclusion that, at least in certain in 
stances, the same unit hydrograph.may be used for rainfall and snowmelt run 
off. The average daily discharge in cfs for melting beginning at four differer 
times during the day are given in Table VI. Fig. 4b compares the recorded 
hydrograph with the reconstructed hydrograph for the snowmelt period of 19 
by means of the procedures outlined herein. 


Ice Index Method of Altering Published Data 


According to the United States Geological Survey Water-Supply Paper No. 
1085 (1947) the runoff records for the Big Eau Pleine River near Stratford, 
Wisconsin are “. . . good except those for periods of ice effect and those bek 
10 second-feet, which are fair.” The notation “good” “. . . indicates that, in 
general, the error in the daily records is believed to be. . . less than 10 per 
cent; ‘fair’ less than 15 percent;. . .” 

Several times during each year actual discharge measurements are made 
at Stratford and compared with the apparent discharge as determined from t 
concurrent gage height and the rating table for the station. Frequently a hea 
correction of a few hundredths of a foot is determined from these actual 

TABLE VI 
AVERAGE DAILY DISCHARGE FOR 
SNOWMELT UNIT HYDROGRAPH 


Discharge in cfs for Melt Beginning at: 


3 am 9 am 3 pm 9 pm 
Day (Period 1) (Period 2) (Period 3) (Period ). 
1 2753 1603 326 5 
2 2218 3068 3937 3667 
3 667 868 1137 15h7 
4 218 282 St 
5 97 19 140 168 
6 6 59 68 80 
7 18 22 28 38 
8 1 hk i 12 
Total 6018 cfs days 6025 6020 602) 
Ratio of 2753/326 1603/326 326/326 5/326 
"First Day" or or or or 
Discharges Boh 492 1.00 0.0153 
(Based on 
Period 
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discharge measurements which correction is then applied to the computed dis: 
charges for a period of time as determined by earlier or subsequent discharg 
measurements. It is quite natural that these corrections should be required 
as a stream bed tends to erode and redeposit its sediment in a slightly differ: 
ent fashion following each passage of high water. The variations however ar 
slight and no alteration of the rating curve has been necessary since 1948, at 
which time only minor changes in the discharges were made for gage heights 
less than 4.5 feet. 

During the winter season ice forms on the surface of the stream and the 
discharge characteristics change from those of open channel flow to those of 
enclosed flow. During periods of prolonged cold weather the apparent dis- 
charge as computed from the recorded gage heights and the rating table will 
tend to increase markedly, while true discharge as determined by field 
measurements is actually decreasing. This is explained by the fact that the 
resistance to flow is markedly increased by the surface contact with the ice, - 
while at the same time, the weight of the ice supported-by the water causes al 
increase in the water level measured in the float well. As an example, an 
apparent discharge of 110 cfs may be determined from the gage reading at 
Stratford, while an actual measurement by current meter may indicate that 
the discharge is 10 cfs. 

Usually two or three discharge measurements are made each year during 
the period when ice effects are noted. The critical period during early snow 
melt runoff, however, is not well represented by actual measurements becaus 
rising water makes the ice unsafe for measurement operations through the 
ice while complete cover exists, and floating ice renders the taking of velocifj 
measurements impractical after the break-up begins. Correct discharge 
values for the early snowmelt runoff season are the most difficult to de- 
termine and are the ones most likely to be in error in the published data. TI 
USGS engineers determine the discharge values reported in the Water-Supply 
Papers from the gage height readings, concurrent air temperatures, reserv 
storage and discharge records, and hydrographs of previous runoff. The 
method of accounting for ice effects presently being used on the Big Eau 
Pleine River by USGS engineers consists of plotting the apparent discharge, 
as determined by recording-gage readings, on semi-log paper, along with tht 
actual discharges determined by the winter discharge measurements which 
are made about two times each winter. A line is then drawn through the 
measured discharges approximately parallel to the apparent discharge and th 
discharge values to be recorded are read from this line. Fig. 5 shows the 
manner in which winter discharge values were determined by this process fo 
the winter 1950-1951. Fig. 6 shows the same discharges plotted on an 
arithmetic plot for the break-up period of 1951. It will be noted in Fig. 6 tha 
at the higher discharges, the use of the USGS method results in extremely 
large reductions in discharge for ice effects. It has been deemed advisable 
therefore, to alter certain critical discharge values from those given in the 
Water-Supply Papers for the purpose of this investigation in accordance with 
the ice index method outlined below. 

All winter discharge measurements on the Big Eau Pleine River were ané 
lyzed. The gage reading which would have been required to produce the actu 
discharge as determined by measurement (without ice effect) was obtained — 
from the rating table. This reading was subtracted from the measured gage 
reading to give the “ice index”. It is reasoned that the ice index is equal to 
the equivalent water depth of the ice cover plus an increment of thickness to 
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FIG. 6 


compensate for the reduction of flow caused by the drag of the water on the ~ 
ice cover. Values of ice index were computed for all days when winter dis- — 
charge was measured, and are plotted against the date of the measurement in 
Fig. 7. The maximum ice index value determined by actual measurement is 
1.65 feet. 

It is reasoned that the ice indices of Fig. 7 may be applied to the recorded 
gage heights to determine the probable discharge for those days when only th 
gage readings are available. For example, on March 17, 1945, the average 
gage reading of 12.95 feet was recorded. This is equivalent to a discharge 0: 
5890 cfs without ice correction as determined from the rating table. Ifa 
maximum ice index of 1.60 feet were considered to apply, the corrected gagé 
reading of 12.95 minus 1.6 or 11.35 feet would provide a flow of 4315 cfs for 
this date. The discharge reported for this date in USGS Water-Supply Paper 
No. 1035 is 1500 cfs, which would occur normally with a gage reading of 7.64 
feet as determined from the rating table. Thus the computed ice index for — 
this date would be 5.31 feet. Such an index could be correct if an ice jam oc: 
curred just below the gaging station and backed up water at the station; or if 
an ice canopy about 5.31 feet in thickness were formed over the stream at th 
gaging station for the full channel width. The absence of any obstruction to 
flow, at the gaging station or downstream, which might cause large ice dam 
to form leads to the conclusion that such reductions in flow as a result of icé 
effect are undesirable for the purposes of this investigation. 

When published discharge data were examined from all the years of reco 
for which good gage readings are available in the initial snowmelt period, it 
was found that certain inconsistencies exist, relative to a study by the ice © 
index method. During the days immediately preceding the occurrence of snd 
melt runoff, the ice index had a reasonable value between 0.7 and 1.7 feet 01 
all years. During the early snowmelt period the maximum ice effect based 
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published discharge data, varied from year to year, at times reading maxi- 
mum values as high as 7.10 feet as in 1938, and as low as 0.74 feet in 1946. 
Furthermore, in some years the ice index was considered to become zero on 
the rising side of the hydrograph, at a discharge of 109 cfs as in 1946, or 
8950 cfs as in 1945; and in other years it became zero on the falling side of © 
the hydrograph, at a discharge of 2240 cfs following a peak apparent discharg 
of 8950 cfs as in 1938. In addition, the ice index for some years was erratic 
alternating between low and high values on successive days. ; 
To eliminate these variations, the following rules for reconstituting the 
average daily discharge for days in the early snowmelt period were adopted. 


1. In the hydrograph summation computations the ice index correction at 
the beginning of snowmelt was assumed to be about equal to that de- 
termined from a previous measurement and was never more than 1.65 { 
feet. 

2. The ice index was assumed to approach zero on the rising limb of the a 
hydrograph at a gage discharge of less than 5000 cfs. 

3. The exact discharge for the last day of ice effect was determined fron 
the hydrograph summation analysis of the total seasonal flow, as out- 
lined later. 


Determination of Daily Melt 


The analytical determination of the amount of snowmelt required each day 
to reproduce the recorded snowmelt hydrograph is a trial and error pro- 
cedure. The complete analysis for the year 1943 is presented to illustrate # 
technique developed. 

Basically, the method consists of applying a runoff for any given day whic 
will make the total runoff for that day exactly equal to the average daily dis 
charge as determined from the USGS Water-Supply Papers, altered, as nece 
sary, for assumed ice conditions. The runoff so applied for any day is then 
followed by a sequence of discharges for succeeding days as determined f 
what has been termed the ‘relative hydrograph”. Thus, the magnitude of fir: 
day’s runoff, resulting from a given period of melting, is of primary im- 
portance because it automatically determines the magnitude of the discharge 
for succeeding days. 


Relative Hydrograph 


The relative hydrograph for any day may be defined as the succession of 
average daily discharges resulting from melt occurring in the proportions — 
given by the variation in melt potential for that day, and so related that the 
runoff corresponding to period 2 is equal to one inch. Since the melt potenti 
variation for each day is usually different from that of every other day, it — 
follows that each day will have a different relative hydrograph depending up¢ 
its melt potential variation. Table VII lists the melt potential data for 1943 
which is necessary for the computation of relative hydrographs, and a summ 
ry of the final runoff data computed by the method of hydrograph summation 
discussed later. The relative hydrograph computations are tabulated in Tak 
VII and are determined in the following manner: 

From Table VII, beginning with March 25, 1943, it is noted that the melt 
potentials for periods 2, 3 and 4 are 0.80, 0.36 and 0.01 inches respectively 
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TABLE VII 
MELT POTENTIAL AND RUNOFF = 19,3 
Melt Potential in Inches Runoff in Inches 
Date Pe To ative Total  Cumilative 
(1) (2) (3) (h). Total Total 
March 
23) 263 063 
) Ripening Period 
075 1.38 eOL 201 
280 236 201 1.17 2055 202 203 
26 201 053 053 3.08 208 ell 
27 ell ell 3219 022 033 
28 026 026 3eh5 220 053 
29 265 205 228 298 eoh3 1.33 1.86 
30 el9 1.26 e72 260 2079 7022 2.2h 4.10 
31 oA 022 1.16 8.38 258 he68 
April 
af Pani ol 8.82 202 he70 
2 e31L o3l 9.13 015 4.85 
3 3h =O 038 9051 007 ho92 
k 1.59 °76 2035 11.86 026 5.18 
TABLE VIII 


RELATIVE HYDROGRAPHS - 19)3 


Relative Average Daily Discharge - cfs Relative 
Date Period First Second Third Fourth Fitth Sf Se Third Fourth Fifth Sieth Seventh Wenth Eighth Runoff 
Day Day Day Day Day ODay Day Day Inches 
Cc. 
25 2 1603 3068 868 282 119 Come h 1.000 
3 BU hyp healt hy 512) 170) 63 Elo! yale) 3} e450 
k (@) 6 19 6 2 al 2013 
Total 1750 886 1399 58 18) 91 35 U 1.463 
29 2 1603 3068 868 282 119 59 = 22 hk 1.000 
’ 3 25 303 87 29 1a 2 1 0977 
h 2 1579 666 218 72 3h = 16 5 0431 
Total 1630 950 1621 529 202 98 Ko) 10 1.508 
30 a o8 329 99 32 WU 7 3 18 
2 1603 3068 868 282 119 59 = 22 h 1.000 
3 183 2215 640 212 19 38 16 h 0563 
hk Dao w® 725%. 238"0 179 38 = «18 6 69 
Total. 2196. 7332 2332 764 .291 a2 59 ly 2.180 
31 2 1603 3068 868 282 119 59 = 22 h 1.000 
3 NG OT, W266 ek 38 16 7 2 023) 
Total 1679 3989 113) 370 152 16 Saale) 6 1.234 
lpril 
3 2 1603 3068 868 282 119 59 = 22 hk 1.000 
3 (ee 17 ame Be TR Fe 8 7 0118 
Total 1641 3531 1002 326 135 67 25 5 1.118 
h 2 1603 3068 868 282 119 59 = 22 4 1.000 
3 1567) £882) 95143» 8180 2867 SBS 5) 3} 0478 
Total!) ..1759 L950" & 462 186 928. 935 7 1.78 
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Referring back to Table VI it will be noted also that one inch of runoff for 

of periods 2, 3 and 4 would produce average “first day” discharges of 1603 
cfs, 326 cfs and 5 cfs respectively. To combine these discharges in the ratio 
of the melt potentials listed above, the “ratio of first day discharges” line of 
Table VI is utilized. The “first day” discharge of period 2 is 4.92 times as 
large as that of period 3, and the “first day” discharge of period 4 is 0.0153 
times that of period 3 for equal runoff, so the “first day” discharge for the 
relative hydrograph will be a combination of discharges from periods 2, 3 

4 in the ratio: 


4.92 x .80: 1x .36: 0.0153 x .01, or 
3.94 : 0.36 : 0.000153. 


For convenience in determining the relative hydrographs the runoff for 
period 2 was assumed equal to one inch in all cases. Thus the relative “first 
day” runoff for period 2 is 1603 cfs. Proportionally: 


1603 period 3 discharge 
3.94 © 0.36 4 


or Period 3 discharge = at (0.36) = 147 cfs. 


and for period 4 the “first day” discharge is: 


1603 

3.94 (0-000153) = 0.062 cfs 
Period 4 computations are carried to at least two significant figures because ~ 
of the extremely large increase in second day discharges for this period as 
compared to the first day discharge. Thus, the relative hydrograph for Mar 
25, 1943 is composed of first day average discharges for periods 2, 3 and 4 
of 1603, 147 and 0.06 cfs respectively, as shown in Table VII. For period 2 
the average discharges for succeeding days are: 3068, 868, 282 cfs for the 
second, third and fourth days respectively as determined directly from Table 
VI, since it was assumed that runoff for period 2 was equal to one inch in all 
cases. For period 3 the discharges for succeeding days are “f” times the 
tabulated discharges for period 3 given in Table VI, where: 


Therefore: 0.450 x 3937 = 1772 cfs for the second day, and 0.450 x 1137 = 51 
cfs for the third day, etc. as shown in Table VIII. For period 4 the same 
process is repeated, with f = 0.0126, and the discharges for succeeding days 
are 46, 19, 6, 2 and 1, as shown in Table VII. When the above discharges aré 
summed up by days for periods 2, 3, and 4, the relative hydrograph for March 
25, 1943 is obtained as indicated by the “Total” line. Thus the relationship 
between daily discharges as a result of melting occurring on March 25, will 
be in the proportions given by the relative hydrograph. 

The “f” values determined above are equal to the runoff contribution of 
each period to the relative hydrograph. The total runoff represented by the 
relative hydrograph, therefore, is: 1.000 + 0.450 + 0.013 = 1.463 inches. It — 
will be noted that the ratio of runoff per period, 1:0.450:0.013 is exactly 
equivalent to the ratio of melt potential per period, 0.80 : 0.36 : 0.01 which 
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as determined originally from the melt potential variation for March 25, 
ables VII and VIII list the complete computations for relative hydrographs 
ir the year 1943, 


ydrograph Summation 


The process of hydrograph summation, by means of which the daily contri- 
itions to the total seasonal snowmelt hydrograph are determined, can also 
= explained best by the presentation of an example. Therefore, the trial and 
rror procedure is presented below as it was computed for the year 1943. 
The average daily discharge for each day of the snowmelt season was tabu- 
ited in sequence as shown in Table IX, Trial No. 1. These discharges are: 
1, 24, 58 and 204 for March 23 through 26, respectively. Base flow dis- 
larges were assumed in accordance with principles not discussed herein. 
ince a base flow of 8 cfs was assumed for March 23 and the total flow for 
lat date was 11 cfs, there remained 3 cfs which was presumed to have been 
ipplied by melting snows. Since this melting occurred in Period 2 (Table 
II), discharges for subsequent days were proportional to the discharges for 
2riod 2 listed in Table VI and were 6, 2 and 1 cfs for the second, third and 
urth days, respectively. These values placed in Table IX left a deficiency 
| the discharge for the following day, March 24, of 12 cfs, which was pre- 
imed to have been contributed from snowmelt. Subsequent discharges, pro- 
rtional to period 2 discharges in Table VI of 28, 8, 3 and 1 cfs were comput- 
i and listed in Table IX. The deficiency in discharge for March 25 was then 
3 cis, which presumably resulted from melting. The relative hydrograph for 
elt occurring on the 25th is given in Table VIII, as discussed previously. 
ince the first day’s flow for the 25th of March was 23 cfs, the second day’s 
charge was: 


23 


1750 (4886) = 64 cfs, 


id the third day’s discharge was: 


23 
1750 


e values 4886 and 1399 being the second and third day’s relative hydrograph 
ilues as given in Table VIII. 

The process outlined above was continued until April 1, at which time the 
immation of discharges was found to exceed the recorded discharge by about 
8 cfs. This excess was presumed to be the result of an error in the as- 
imed base flow distribution, or an error in the average daily discharge data. 
he assumption was made that the discharge value most likely to be in error 
as the discharge corresponding to the last day of ice effect, which was 
arch 28, 1943. The discharge for March 28, used in the first trial, was 1066 
s, which was the discharge computed from the gage reading for an ice index 
1.37 feet. Experience with the method indicated that an increase in the dis- 
arge for March 28, of about 34 cfs might result in a balance and a second 
ial was begun using a discharge of 1100 cfs for this date. 

Results of the second trial are also shown in Table IX. It will be noticed 
at all discharges from melt which occurred before March 28 were unaffected 
‘the change under consideration and could, therefore, be lumped together as 
summation of previous discharge. The second trial resulted in a discharge 


(1399) = 18 cfs, 
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summation of 2394 cfs for April 2, which was 584 cfs more than the recorded 
discharge for this date. It was evident that an over-correction had been ap- 
plied to the discharge of March 28, and a discharge of 1085 cfs was selected 
for this date for the third trial. 

In the third trial the discharges were reasonable until April 6, at which 
time the summation of discharges was 31 cfs too large and a fourth trial was 
undertaken, based on a discharge for March 28 of 1084 cfs. 

The change of one cfs in the assumed March 28 discharge was an over- 
correction and the fourth trial resulted in excessive discharge for April 5 of 
75 cfs. Since it was impractical to apply a further correction to the discharge 
of March 28, the discharge for the second day resulting from melt which oc- } 
curred on March 30 was changed from 7535 to 7545 and the fifth trial was ; 
computed which balanced out satisfactorily. It was assumed that all discharge 
required for balancing after April 7 was the result of base flow and rainfall | 
runoff, 

The flow of 1084 cfs for March 28 is comparable to an ice index of 1.11 : 
feet. Discharges prior to March 28 were computed using an ice index of 1.37, 
which was determined by actual measurement on March 24. The ice index 
was zero after March 28. 

To compute the runoff attributed to any day, the ratio of the “first day” dis- 
charge from Table IX to the “first day” relative discharge from Table VIII 
was multiplied by the total runoff represented by the relative hydrograph. 
Thus, on March 25 the runoff was: 


(1.463) = 0.02 inches, 


ma a 
and on March 30 the runoff was 

2257 

3196 (2. 180) = 2.241 inches | 

After the daily runoff values were computed they were added cumulatively 

from the beginning of snowmelt as recorded in Table VII for the year 1943. ; 
The positive melt potentials were then added cumulatively for a period begin- 
ning two days before the occurrence of snowmelt runoff, as listed in Table b 


Vil. 


Precipitation 


The occurrence of precipitation during the snowmelt season increases the 
difficulty of the melt potential-runoff studies many-fold. This is the results 
of several factors. First, there is no recording precipitation gage within the 
Big Eau Pleine River Basin. Therefore, the precipitation data, used from 
adjacent stations, may be grossly erroneous. Secondly, records of snowfall 
and snow on the ground are available only after October, 1949. Thus it is 
often difficult or impossible to determine whether the recorded precipitation 
was rainfall available for immediate runoff, or snowfall for which a signifi- 
cant melt potential would be required in order to produce runoff. 

Significant amounts of rainfall occurring in the snowmelt period were 
handled in two ways. Whenever possible the rainfall was added directly to the 
melt potential factors for the same period in the analysis. Hourly precipi- 
tation values were used for the periods for which they were available, and 
total precipitation was divided equally for the day for those periods where 
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hourly data were not available. For some years this procedure would not 
permit a solution and it was necessary in these cases to consider that all of 
the rain falling in a given period ran off. Accordingly, the hydrograph of daily 
values was calculated for each increment of precipitation and their sums sub- 
tracted from the total runoff to leave the runoff attributed to melting snow. 
The years 1942 and 1944 were omitted from further consideration because the 
precipitation during the snowmelt season was greater than the amount of snow- 
melt runoff. 


Correlation of Results 


As mentioned earlier, if the techniques developed herein are to be con- 
sidered valid, a reasonable correlation must be evident between the runoff and 
the melt potential. When the cumulative melt potential for the snowmelt 
season is plotted against the cumulative runoff, as tabulated in Table VII for 
1943, a marked similarity among the curves for the various years is noted. 
Fig. 8 shows the curves plotted with the occurrences of significant initial run- 
off at zero on the abscissa scale. It will be noted that the approximate maxi- 
mum initial rate of runoff is 2 inches over the watershed area for the first 
1.5 inches of melt potential, while the minimum initial rate is approximately 
1.35 inches of runoff for the first 1.5 inches of melt potential. In all cases the 
steepest runoff gradient coincides with the earliest portion of runoff, as would 
be expected, because the snow cover is most extensive in the early melt 
season. For most years the melting is complete, or nearly so, at a melt po- 
tential of 7 inches beyond the occurrence of significant initial runoff. If the 
curves of Fig. 8 are reconstructed with a common ordiante at a melt potential 
of 7 inches the curves of Figs. 9 and 10 result. The seven years having the 
most reliable data are plotted in Fig. 9 and a very good correlation is evident. 
In Fig. 10 the years of less reliable data are plotted. The snowmelt-runoff 
equation, represented in Figs. 9 and 10 by the small circles, is: 


Y = 27.4% - 1.88x2 


where Y is the per cent of the runoff at a melt potential of seven inches and 
X is the cumulative melt potential in inches after the beginning of the runoff. 

It was noted from a study of the melt potential runoff data that initial run- 
off did not occur at a fixed accumulation of melt potential, but varied within 
rather wide limits. 


CONCLUSIONS 


1. The relationship between the cumulative melt potential and the total 
snowmelt runoff, taken as the runoff at a melt potential of seven inches, 
may be expressed by the equation: 


Y = 27.4X - 1.88X2 


for the Big Eau Pleine Basin, where Y is the per cent of the runoff at a 
melt potential of seven inches and X is the cumulative melt potential 
after beginning of runoff. 

2. Apparently, factors causing snowmelt may be measured with less than 
ideal precision at locations somewhat removed from the point of appli- 
cation without losing their effectiveness. This may be due largely to 
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the compensating effect of the various errors which are encountered in 
Hydrologic data. 

3. The general problem of ice effect on discharge rates in the spring 
breakup requires further study. The assumption seems valid that actu=— 
al discharge rates for the Big Eau Pleine River at Stratford were higher 
than the rates published by the U. S. G. S. in its Water-Supply Papers; 
however, the ice index method presented must be subjected to further 
study on other rivers before any sweeping claims may be made for its 
validity in other areas. P 

4. Rain falling on the snowpack early in the snowmelt season has not pro- — 
duced rates of runoff as high as those produced by snowmelt alone, for © 
the period of record examined. It must be emphasized, however, that ~ 


and perhaps more intensive instrumentation on small plots. 
6. It appears that the coefficient K in Light’s (1944) equation: 


D = KV(T-32°) 


has a reasonable value of K = 0.002 for the melt conditions of the 
Wisconsin river valleys. 

7. The nomograph presented in Fig. 3 is a useful tool in determining the 
snowmelt caused by insolation, and is as accurate as the cloud cover 
data which is used in it. It is believed applicable to latitudes from 40° — 
to 50°. 

8. The snowmelt unit hydrograph developed from snowmelt runoff may be 
considered equally applicable to the study of snowmelt runoff and rain-— 
fall runoff. 
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CHARACTERISTICS OF A LARGE THROATED SIPHON@ 


Closure by J. C. Stevens 


oC. STEVENS, ! M. ASCE.—The writer is grateful for Mr. Rice’s very 
uminating discussion. It is particularly valuable since the laboratory model 
idies were made by him, at Oregon State College. 

Mr. Rice has resurrected the ghosts of the definitions of discharge coef- 
ients and efficiencies of a siphon. Webster’s International Dictionary de- 
es the latter thus: “1. Efficient quality, power or action. 2. Mech. The 
tio of the energy or work that is got out of a machine, a storage battery, an 
sctrical distributing system, or the like to the energy put in...” 

In the writer’s earlier writings (footnotes 4 and 5) he thought he could 

ttle these moot questions once for all but perhaps that was expecting too 
ich. Mr. Rice concedes that the writer’s definition of efficiency there given 
the “least objectionable, even though it varies with varying heads on the 
ucture.” But how many devices have a constant efficiency? The writer 

a think of none. 

The discharge coefficient will vary with each change of cross-section. For 
S reason it seemed desirable to fix on one, and the outlet was chosen. It at 
st offers means of comparing siphons under comparable conditions—for 
stance, under maximum and minimum heads. 

Testing a model under full head is an excellent method of forecasting the 
iracteristics of a projected prototype. Had the laboratory model been 

ong enough for full-head operation much more valuable data could have 

=n secured in advance of construction. Table III indicates some of the po- 
tialities in this method even crude as the facilities were for: testing and 
ipled with the weakness of the model itself. Such tests render unnecessary 
»use of controlled atmospheric cabinets for sub-atmospheric pressure 
dels—such cabinets have been built. 

The writer cannot agree that the method of correcting the pseudo-energy 

id of Fig. 5 is necessarily erroneous. Mr. Rice indicates that the maximum 
ue of the alpha coefficient for velocity head is 2.0 and that any higher values 
uld imply a reversal of flow. The writer firmly believes that that was ex- 
ly what happened. The highest value in Table VI is 3.5 at the inlet. Here 
velocities were moderate but the sharp right-angle turn the water had to 
ke to enter the siphon (Fig. 1) produced marked eddies with pronounced 
rersals of flow. A value of 3 is indicated at piezometer sets 3 and 4. No. 3 
just below the sharp curve at the end of the crest. Here again marked 
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eddies with flow reversals surely obtained and doubtless these continued on 
Set No. 4. 

The writer will agree that this method of constructing the energy line 
water passing through a tortuous siphon has elements of arbitrariness in it. 
However there are two facts that cannot be ignored. (1) The energy line mi 
slope continuously downstream since energy is being lost continuously, and 
(2) there is hardly any point within the siphon where the true mean velocity 
equals Q/A, hence the coefficient alpha is everywhere greater than unity, 
though Table VI shows it as 1.00 at piezometer set No. 2. This is in the 
middle of the long crest where it would naturally be nearer unity than any- 
where else in the siphon. 
, Mr. Rice’s velocity traverses of flow in the laboratory model cannot poss 
bly represent what actually occurred in the prototype, with its violent eddi 
and direction reversals of velocities. It is possible however that had it 
practicable to have made such traverses in the high head model a greater 
conformity to the prototype might have resulted. 


ERRATA 


Paper 1198 as printed in HY 2 April 1957 contains errors as follows: 


Page 1 under synopsis 3rd line change primary to priming 
2, line 16 change floods to boards 
line 10 from bottom change cross-section to profile 
' 5, 3rd line after the table change Figs. 1 and 2 to Fig. 2 
11, 1st line of the table change 56 in date to 55 
Eq. (6) change final zero to letter O 
12, Table IV heading col. (1) change M to m 
15, Table V 1st col. change 6 to 60.0 
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TRANSITION FROM LAMINAR TO TURBULENT FLOW IN A PIPE@ 


Discussion by J. M. Robertson 


J. M. ROBERTSON,! M. ASCE.—The subject of laminar-turbulent tran- 
tion has received little attention in ASCE publications, although it has been 
\scussed in the general fluid mechanics literature since the time of Osborne 
eynolds and Lord Rayleigh. This lack of concern of the civil engineer stems 
‘om two main causes. Firstly, most theoretical approaches have been in re- 
urd to problems (such as boundary layers and laminar-flow airfoils) of 
rimary concern to the aeronautical engineer. Secondly, the experimental 
formation on transition of primary interest to the civil engineer interested 
| pipe flow was obtained over half a century ago. In fact, Reynolds’ classic 
<periments were conducted seventy-five years ago, while over a hundred 
oars ago G. H. L. Hagen noted a temperature effect on the occurrence of the 
‘ansition between, what may now be termed, laminar-viscous and turbulent 
ow. 2 Although Reynolds did not use our present day terminology, his studies 
d indicate that for “engineering pipes” the critical Reynolds number has the 
(ow) commonly accepted value of 2000 to 2500. This agrees with Hagen’s 
st results (for pipes with square entrances) which have been shown to indi- 
ute a critical value of 2200 to 3000. Since the subject appears to be rather 
dary, this writer will attempt to justify current interest in the transition 
roblem. 

It is usually stated that there is an upper critical Reynolds number above 
hich the flow is always turbulent. Actually there is some uncertainty in this 
ssard. Experimentally, it is observed that as we decrease the disturbance 
itering with the flow, the critical number increases rapidly. Reynolds him- 
If, in what we now call a “Reynolds apparatus” (transparent pipe with bell- 
outh entrance from a rather large tank), had little difficulty achieving lami- 
ir flow at Reynolds numbers R = VD/v of as high as 12,000 to 13,600. Other 
sperimentors have been able to reach even higher values with no great 
ouble. Thus in 1909 using Reynolds’ original apparatus, but taking particu- 
r care with the bellmouth and allowing the water in the tank to settle for up 

half a day, V. W. Eckman(3) raised the upper critical number to a value of 
),000 to 50,000. According to Table II the author has found laminar flow at 
synolds numbers of about 15,000 to 51,000 depending upon the spot source 
id other flow factors. These high values are a tribute to his experimental 
ill, since there is some evidence that a timewise developing flow may be 
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less stable than a steady flow. One wonders how high a Reynolds number he™ 
might have reached without turbulence if the disturbance strengths of his spo 
sources had been further reduced. 
The care required to achieve laminar flow at such high Reynolds numbers 
may appear impractical to the hydraulic designer—however the water ina 
reservoir often has ample time to become quiescent. Furthermore, lamina 
viscous flow may occur even in pumped systems at Reynolds numbers con- ~ 
siderably above the accepted critical values. Thus, we have the results re- 
ported by the University of Iowa 4) with a system in which oil is pumped cor 
tinuously through a pipe; it is impossible to obtain turbulent flow even at the 
upper limit of the pump capacity (R about 7800). The turbulence must be art 
ficially triggered to make it appear in this case. It would appear that ap- 
preciable amounts of laminar-viscous flow may well occur in some of our hy 
draulic structures of clean design. 
The author’s study of transition is therefore of potential value to designe: 
as well as being of intrinsic interest to fluid mechancians. One has only to 
glance at the discussions concerning laminar-flow shapes in aeronautical 
publications to visualize the potential gaines. Fifty to one hundred-fifty per 
cent increases in aircraft ranges from the use of suction-stabilized laminz 
boundary layers are suggested. If comparable achievements are possible in 
some of our hydraulic structures, it might be possible to justify the special 
construction and precautions needed to achieve the results. Only through 
fundamental studies such as reported by the author can one obtain the 
knowledge needed to assess what is possible. Also, we need to be able to tel 
if our pipe flows may occur with appreciable laminar flow regions. In many 
cases one may not design for laminar flow but one should consider its likeli 
hood of occurrence. Differences in terms of pressure drop, heat transfer, 
etc., between a design based on the conventional assumption of turbulent flo : 
(at Reynolds numbers above 2500) and a flow with a considerable laminar- 
viscous region can be significant. 
Study of the laminar-turbulent transition in pipes is of considerable in- 
trinsic interest, since there is some question whether it should occur at all. 
there is no disturbance. As noted above, the experimental evidence shows i 
upper limit to the critical Reynolds number as the external disturbances aré 
reduced. One therefore turns to theory to see what says about this question. 
The author has discussed his observed transition results in terms of the pre 
diction of the small-disturbance theory for the stability of Tollmien- 
Schlichting waves, but for the case of Blasius flow, i.e., for the laminar- 
viscous boundary layer on a plane with zero pressure gradient. The result i 
certainly suggestive, even if not applicable to the present problem. This 
small-disturbance theory has been quite useful in predicting the instability 0 
laminar-viscous flows. The instabilities and resultant wave length—Reynold 
number relations (cf Fig. 9 of the paper) do not immediately indicate tran- 
sition to turbulence, but they have been found to exist. Although such pre- 
dictions were long considered of doubtful physical reality, they were first ve 
fied for flat plate flows in 1940. The phenomenon of selective amplification 
disturbances of the proper wave length and frequency in laminar-viscous flo 
is now accepted as a prerequisite to transition, except where relatively lar 
disturbances are present to override this occurrence. 
Prediction of a critical Reynolds number for pipe flows on the basis of th 
small-disturbance theory has been a rather elusive thing, considering the — 
practical interest in transition and the available experimental information. 
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nongst the early analysis was W. M. F. Orr (1906) who utilized the ap- 
caches of Lord Rayleigh (1879 et seq) and O. Reynolds (1895)—in fact the 
sic relation of the small-disturbance theory is often called the Orr- 
mmerfeld equation. Orr found a critical (for instability of the wave dis- 
‘-bances) Reynolds number of 117 for pipe flow. However, more refined 
alyses have indicated this result to be inaccurate. Thus, C. C. Lin(5) states 
0 theoretical indication of instability has ever been found for the fully de- 
loped parabolic profile in the case of the circular pipe.” 

Only recently has the apparent dichotomy, between our expectation of a 

ite critical Reynolds number for pipe flow and the theoretical prediction, 
en explained. The clarification stems from consideration of the fact that 

ly developed flow (parabolic velocity profile) in the pipe does not appear 

1 blown at the entrance. Instead it develops in an entrance region whose 
igth depends upon the mean flow Reynolds number. About 0.06 R diameters 
required for the entrance length; thus 120, 600 and 3000 diameters are re- 
ired for the flow to develop at Reynolds numbers of 2000, 10,000 and 50,000, 
spectively. In this entrance region the flow is of a boundary-layer type with 
small favorable pressure gradient. Because of this boundary layer type 

w the author has some justification for comparing his results with the sta- 
ity predictions for Blasius flow. For this latter case the most refined 
culations indicate a critical Reynolds number (lowest R value on neutral 
ibility curve of the author’s Fig. 9) of R 6* = 420. The boundary layer of the 
trance-region flow differs from the Blasius boundary layer due to the 
essure gradient. 

The possibility that the entrance-region boundary layer may be unstable 

en though the following Boiseuille flow is stable was investigated by T. 
tsumi(6) in 1952. In terms of his solution for the velocity profiles of the 
rrance-region boundary layer, his calculations yielded neutral stability 
rves, i.e., indicated the existence of flows unstable to certain disturbance 
ves. The minimum critical Reynolds number (R = VD/v) was found to be 
400 at a station rather near the entrance, namely at about x/D = 0.0004 R 
67% of the length of the entrance region) or about 7 D from the start of the 
ye. For flows at higher Reynolds numbers there is a portion of the entrance 
tion which is unstable and in which small disturbances, if of the proper 

ve length or frequency, will be amplified. Thus for a Reynolds number of 
000 the flow is “unstable” (capable of amplifying certain small dis- 
‘bances) for x/D distances between about 1.4 and 20 as may be seen from 
idy of Fig. A of this discussion. If not enough amplification of the dis- 
ances occurs in this distance and they do not become large enough to 

ise transition to turbulence, then Tatsumi’s analysis would say the flow is 
ble. It is apparent that the higher the flow Reynolds number (above 20,000) 
-more opportunity there is for amplification of small disturbances and the 
sater is the possibility that transition to turbulence will occur. A study of 
-author’s transition data, in terms of the parameters of Fig. A, indicates 

t the transition Reynolds numbers were generally above the minimum criti- 
value for stability and tended to increase with x/D. However much larger 
tical x/D values occur. 

The author discusses the fact that the disturbances must be considerably 
plified after they pass the critical Reynolds number (for neutral stability) 
ore turbulent flow sets in. His discussion is in terms of the stability analy- 
of the Blasius flow rather than the actual case in hand. It would appear 

re appropriate to compare the experimental results with the predictions of 


gs 
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the stability analysis of Tatsumi even though this is only for the steady flow ~ 
case. Actually the comparison should be with neutral stability curves for the : 
profiles calculated from Eq. (3) for the temporally and spacially growing . 
boundary layer. However, neutral stability curves are not easily determined. 
Since the stability condition is dependent upon the shape of the velocity pro; — 
file, a suggestion as to which known set of solutions might best apply should 

be apparent through comparison of the velocity profiles. Lacking this infor- 

mation, it is still possible to glean some useful information from comparison | 
with the Blasius case solution. The remainder of the writer’s remarks will 

follow the author in this regard. i 

The author’s presentation of his transition data in Fig. 10 is interesting and 
suggestive. However there are several points about the possible neutral sta- 
bility curves shown that the present writer cannot follow. These points have 
to do with the wave length A chosen as the basis of calculation for the curves. 
Thus 2.1 r, is chosen for spot 2 as “the longest disturbance wave length which 
originates at the pipe junction.” We know very little about the discontinuity 
presumed to exist at the pipe junction and if we did we still could not predict © 
anything about the wave length of the disturbance which it would produce. 
Hence, the present writer is forced to conclude that the author chose a wave 
length which would pass the possible neutral stability curve through the ex-  ~ 
perimental point having the highest transition Reynolds number. The justifi-— 
cation for this assumption is that transition is not expected once the dis- 
turbance has passed the upper branch of the curve. It is also stated that the 
appearance of some of the points (actually seven of the nine points for spot 2) 
below the lower branch of the curve causes no conflict. It is hard to accept 
this statement, since for points in this region the disturbances would have had 
no amplification whatsoever. There seems to be no reason why the wave 
length being amplified should be a unique function of the spot source. Actually 
it is a function of the boundary layer flow as shown by the neutral stability 
curve in Fig. 9. 

The amplification situation for the simpler Blasius case is clarified through 
reference to Fig. B, which is nothing but the author’s Fig. 9 with some ad- 
ditional information. Some of the amplification contour lines (as found by 
S. F. Shen(7)) mentioned by the author have been added, together with straight 
dashed lines indicating the progress of constant wave length disturbances 
across the amplification ridge. The encircled numbers on some of these 
dashed lines indicate approximately the total amounts of anfplification that 
have been calculated for these certain wave lengths. The trends of the lines — 
of constant » are quite suggestive of the run lines on Fig. 10. They are in- 
clined straight lines which would pass through the origin, if extended, since as 
the wave progresses along the flow 6* increases and the ordinate and abscissa 
values increase in constant ratio. Natural transition in the low disturbance 
flows (free stream turbulence intensity less than 0.3 per cent) is found at Rxx 
values above 2000. In the few cases in which natural Tollmien-Schlichting 4 
waves have been measured in terms of their frequency and wave length chara¢ 
teristics, they have been found to occur approximately along the upper branch 
of the neutral stability curve at Reynolds numbers above 2000. It is apparent 
that these have been subjected to amplifications of several hundred times at _ 
least and that they must cause transition rather quickly since damping sets in 
after they pass the upper branch of the neutral stability curve. It would ap- ; 
pear that amplifications of one hundred or less are inadequate to lead to 
transition—the amplification occurring at the critical Reynolds number of 420 
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is small indeed. It is not surprising that we find a considerable gap between ~ 
this critical Reynolds number and that for transition to turbulence. 

Attempts to fit the above concepts and results for small disturbance ampli- 
fication and transition to the results shown in Fig. 10 are not encouraging. 
Two possibilities suggest themselves as explanation for the poor agreement, ' 
The first is that the disturbances causing the transition are not small. The 
second is that the neutral stability curve for Blasius flow is a poor repre- : 
sentation for the actual flow situation. This latter is beyond doubt true, but it 
behooves one to also check on the former possibility. Information on the ef- — 
fect of roughness elements in causing premature transition has been well i 
summarized by H. L. Dryden. (8) Treating the joint discontinuity asa flat 
strip element, the important factor telling us whether it will trigger transition 
is the ratio of the height to the boundary layer disturbance thickness 6* evalu- 
ated at the location of the roughness element. For turbulence spots 2 and 3, 1 
calculations based on the data of Table III indicate that at transition this ratio” 
was of the order of 0.03 to 0.07. Since the change in boundary layer thickness 
could not have been large in the distance from the pipe joint to the transition 4 
location and since Dryden shows that no appreciable effect occurs for ratios — 
below 0.1, it is concluded that indeed the joints were only sources of “small ; 
disturbances.” In the case of the source for spot 1, the situation is quite } 
different although exact calculation is not possible with the available infor- 
mation. The sand roughness initiating this spot was on the entrance bellmouth 
where the displacement thickness is extremely thin. The only question in this” 
case is why the author chose to complicate his study with this avoidable source 
of disturbances and, perhaps, why the sand did not always lead to immediate ; 
transition. ‘ 

The author has been very successful in ferriting out of his rather indirect 
study considerable information on the occurrence of transition. It is unfortu- 
nate that it was not possible to measure the velocity directly. J. Rotta(9) was 
able to use a hot-wire anemomter to record the velocity fluctuations for tran= 
sition studies with air flow at Reynolds numbers up to 18,000. No evidence of 
Tollmien-Schliching waves appear in his results, probably because the hot- 
wire was located only at the end of the pipe. It would be most enlightening to 
see how these appear in the regions where their amplification might be ex- 
pected. Rotta’s measurements do give some interesting indications of an 
intermittent type of laminar-turbulent flow such as noted and discussed by 
Prandtl and Tietj ens, (2 

The writer finds himself at odds with the author in regard to Fig. 7 and 
discussion thereof. The turbulence cannot be “diffused into the laminar 
boundary layer,” since it comes from a breakdown of the laminar-viscous 
flow. Such a breakdown, as opposed to an outside effect, is the whole point 0. 
the small-disturbance theory to which we have been trying to tie the transitio 
occurrences. Fig. 7 furthermore shows the turbulent flow after transition as 
completely filling the pipe cross section where the laminar boundary layer 
flow did not. This is not necessarily so, Assuming that the laminar bounda 
layer profile differs only in detail from that for the Blasius case we may tak 
6 = 2.9 6* and, except for run 20, 6 is therefore less than one-quarter of the 
pipe radius. Although the turbulent boundary-layer may well become four 
times the laminar one in thickness, this does not take place immediately. In 
flat plate flows the length of the transition region is often found to be as grea 
as that of the preceding laminar flow region. Apparently the simplification 
Fig. 7 is satisfactory for the purpose of interpreting the data. 
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In Table II, it appears that the boundary layer thickness and Reynolds 
number occurring for transition of spot 2 is less than that for spots 1 and 3 
(except for run 20). The author’s explanation for this would be appreciated. 
In terms of the indications that the disturbances which are amplified and then 
breakdown to turbulence starts from Rg* values of about 1000 and hence 
2 76*/d values of about 0.17 the author’s discussion of the number of wave 
lengths between the disturbance source and transition would seem to be in 
error. For spot 2 it would appear that the distance is only 3 to 15 wave 
lengths rather than 9 to 40. Also is not the x/D value listed for spot 1 or run 
20 in error? Finally, it seems impossible that for spot 1 of run 21, the 
boundary layer can have a finite thickness at x/D = 0, so this value must also 
be in error? 
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FLOOD FREQUENCIES DERIVED FROM RAINFALL DATA@ 


Discussion by C. O. Clark 
Closure by J. L. H. Paulhus and J. F. Miller 


e. O, CLARK, | A. M. ASCE.—The criteria for designing structures, dams 
ind bridges, subject to rain and floods can be derived from a look at a small 
uudget from an easy chair, or from the arduous, costly, and sometimes un- 
omfortable investigation of all data at and all around the scene in all kinds of 
eather, or from any intermediate combination of observation and non- 
bservation. The consequences of the former, although prevalent, are cata- 
rophic. The mere fact that the subject paper deals with flood frequencies 
erived from something other than floods, and from data outside of the water- 
hed to which it refers, invites a critical examination as to whether it has 
eally produced the flood frequencies. 

The authors have used rainfall data near a site, but outside the subject 
atersheds, to derive a series of calculated flood flows, which calculated 
Ows are not greatly different in their mean value and variation from, and 
ence indicate about the same frequency of flooding as, the actual peaks of that 
me short period. But to reconcile their estimate of normal flood frequency 
ith one observed flood, the authors have to conclude, “data suggest that the 
933 observed flood peak of 11,200 cfs is an unusual event with a return period 
robably much greater than 100 years”. 

The writer finds instead that the 11,000 cfs flood is not at all rare, and that 
is only luck that there were not two or three as large in the sample period. 
his fundamental difference in views is attributed to faults of a common sta- 
stical process by which judgments about frequency are reached by exami- 
tion of a small number of observations, and which in this case was probably 
ed with both sets of data, the set belonging to the stream and the set derived 
om rainfall. 

The statistical process is understood to be applicable to some uses with a 
presentative sample of a large population, and has been thought to yield 
asonable estimates of size of floods of rare flood frequency when the sample 
s fairly large and nothing indicated that it was a poor sample, or non- 
presentative, of the whole. 

In the record of the subject stream, Codorus Creek, there are several 
asons to think that the record is a non-representative sample, some of 
ich reasons involve examination of paths of great storms along the Atlantic 
ast and noting how few hit Codorus Creek as compared with other places. 
lence the data are neither in the rainfall nor the flood experience of a normal 
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sample of storm events; both kinds of data, then, are “biased” or unrepre- 
sentative. Neither does the trend or scatter of the observed data suggest that 
the sample may be successfully fitted with a probability curve of a standard 
type derived from these data alone. This, too, may be correctly said of the 
record of the stream which the writer introduces to demonstrate that the 


normal flood experience of the region, and probably of Codorus Creek, is 
much larger. But neither has this record been used by the writer to select 
the form of frequency curve which he says is regional in application and ap- 
plies to Codorus Creek, too. Of the records, either of Codorus Creek or the 
one introduced, neither is of a sample size sufficient to disprove the form of 
the applied or indicated frequency. All that is shown is how much small 
samples can vary from a curve without proving the curve is wrong or right, 
and without being capable of reliably indicating what curve is right. 

Among the streams of this general area is another stream which from an 
identical quantity of runoff in six hours will produce the same size of peak as" 
the authors assign to Codorus Creek. It is the Passage Creek near Buckton, — 
Virginia, some 90 miles SW of the Codorus Creek, but well within the general 
area of the problems and basins selected by the authors. There is one signifi 
cant difference. This stream had its floods; Codorus Creek did not. At least, 
not in the period 1933-54. 

A comparison of the flood experience of these streams with similar flood 
capacity, but different accidental exposure to opportunity for collecting a 
record during a specific period, is shown on Fig. 1. The circles and crosses 
are the experience, and the only facts. The lines are methods of portraying By 
the normal experience of an infinite period of time in relation to these 
samples. The line attributed to the authors is plotted from Table 4. The lines 
attributed to the writer are from flood frequency studies made in the subject 
area about 1942, by James B. Saylor and the writer, for the purposes of the 
Potomac River Flood Control Survey Report, the essential results of which 
were published by the District Engineer, Washington District, Corps of Engi- 
neers, in his report on that subject to the Chief of Engineers in April 1944. © 
This application is in accord with that study and flood producing capacity of a 
stream of 87 sq. mi. of drainage area and concentration time of 6 hours. } 
(Data to this effect appear on p. 1484 of Transactions ASCE, 1945, Vol. 110.) 

The authors give 74.3 sq. mi. as the watershed of Codorus Creek, and Fig. 
5 indicates a time of concentration of about four hours. So the slightly larger 
drainage area of Passage Creek is to some extent offset by the apparent, but 
surprising, greater “flashiness” of Codorus Creek. In keeping therewith, the 
authors have found Codorus Creek capable of producing a flow of 6700 cfs 
from 1.0 inches of surface runoff uniformly distributed in six hours, and the 
writer had estimated that Passage Creek could produce 6900 cfs from such ai 
event. On this basis, the streams are selected as essentially of the same flo} 
producing capacity, or as having flood frequencies no more different than the: 
two almost coincident lines on Fig. 1. ' 

Now let us return to the evidence that the larger flood of Codorus Creek i 
rare. Or to the evidence that because it lies well above a line fitted to the 
remaining or great majority of points that it should be so regarded. 

Except for the two greater values of the flows for Passage Creek, the ma- 
jority of those data are no greater than the majority of the values for Codoru 
Creek, and there, too, a line fitted to the majority of the data and extended 
would indicate a phenomenal rarity to the larger two of the actually observed 
items. 


ASCE DISCUSSION 1856-15 


20,000 ~ ea 
Flood Experience of Short Periods 
X< Passage Creek nr Buckton, Va., 1932-42 
id Codorus Creek at price. Grove, Fa., eee 
15,000 
10,000 


Passage Creek 
Codorus Creek 


suai a 

Burt meee Fa 

me HEL Geer pile 
HY eae 


Peak Flow.in Oubic Feet Per Second 


1 2 5 10 20 50 100 
Average Return Period in Years 


h 


1856-16 HY 6 November, 1958 


The writer has suggested that there is a pitfall which waits for all of the — 
unwary in frequency calculations from short records, namely that there is no ~ 
record sufficiently normal to be used, nor a sufficient number of observations 
to compensate for errors induced by the abnormalities of a few occurrences, ~ 
What is there to indicate an abnormality in the record of Codorus Creek? 
First of all, a comparison like that here discussed. But what is there to indi-— 
cate that it is Codorus Creek and not Passage Creek, or the writer’s idea of 
normal frequency, that is abnormal? Part of that is in the review of the 
record of this area as a whole. 

Such a review shows that few of the great storms of 1930 to 1950 passed 
very close to Codorus Creek. A further review shows that a storm prior to 
1900, probably 1889, not only passed nearby, but created rumors of great 
floods that persist to this day. Perhaps it was the same storm which sent the 
Potomac so high that backwater in the gorge below the Falls backed right ove 


water mark, there engraved in rock (it was a mark that was discarded as 
falsely high until a lesser flood of September 1944 demonstrated that the 
could, indeed, be drowned out, whether the dam could be or not). 

Continuing in such a review, the 11,200 cfs flood of Codorus Creek occurred 
in August 1933, in general a period of drought, but in a storm that established 
records on a number of nearby Maryland and Pennsylvania streams with shor 
records, but not on any of the streams with moderate length of record at that 
time. The storm is listed as NA 1-24 among the great storms studied as 
“Storm Rainfall in the United States”, published by the War Department, and 
one of the bases for preparing Hydrometeorological Report No. 23, a report © 
of the U. S. Weather Bureau on “Generalized Estimates of Maximum Possible 
Precipitation Over the United States.” 

The storm maps show 6.8 inches over 10 sq. mi. in six hours as the most 
intense part of this storm. Report No. 23 gives 24 inches as possible on 10 
sq. mi. both there and near Passage Creek. The writer does not consider the 
1933 rainfall to be a sufficient approximation of the merely possible as to be 
even unusual. How then, can the flood be unusual ? 

For further comparison, the two larger flood items for Passage Creek are 
for October 1942 and March 1936, Storms SA 1-28 and SA 1-27. In the 
October center, 6.0 inches is the value given for the maximum 10 sq. mi. in 
six hours, but it is shown as occurring after 14 inches of prior rain, as com= 
pared with 3 inches preceding the maximum rates in the NA 1-24 storm. 
March 1936 was such a commonly wet month and year that no rare signifi- 
cance can be assigned either to the 2.7 inches of six-hour rainfall or the 
magnitude of the flow which it produced and which is the second larger on 
Passage Creek. Compared with these not too rare storms, however, the 188) 
storm, SA 1-1, is presumed to have deposited 7.4 inches of water in the samé 
period and area, a wetter season of the year, and somewhat to the northwest 
of Codorus Creek. 

What would Codorus Creek do with 7.4 inches of water in six hours, if one 
inch of runoff produces 6700 cfs? Fig. 3 indicates that from 2.5 to 5.5 inches 
would run off. Wouldn’t that make 17,000 to 36,000 cfs of flow? What frequel 
cy do such flows have? Where would they fit on the suggested normal freque 
cy curves? . 

Out of some 2500 flood items that lie behind the curves which the writer 
presents as normal are over a hundred as large as 2 inches of water in 6 
hours could produce, and a few are over three times that big. For the basins 
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he authors list in Table 1 as Nos. 14-19-20-15-18, values of 5.4-4.1-3.3-2.8- 
2.7 are found in some very short records. In the historical records of the 
same area, but without any guide as to the total number of items among which 
hey would be found, are values as high as 8 inches of runoff could produce. 
All are in the area of study of the subject paper, and in the area for which 
Report No, 23 suggests 24 to 26 inches of rain as being possible on 10 sq. 
niles in six hours (Values for 100 sq. mi. are not given). 

It has been the purpose of this discussion to show that something more than 
he authors have presented is necessary to derive flood frequency either from 
rainfall data or from streamflow data, and that without it either or both kinds 
of data can lead to an erroneous conclusion. There is a wide amplitude of 
yJariation in the magnitude of floods which are possible in the subject region. 
rhe largest item cannot be discarded from the sample only on the basis of the 
remaining experience in the same small sample, nor for the convenience or 
udgment of the calculator. There is nothing needed to explain the wide vari- 
utions which the writer has presented except measures of accidental proximity 
) the storm centers and accidental variations of prior wetness and coinci- 
lence in the degree of storm intensity with that to which a particular stream 
S most sensitive. These alone are sufficient to explain wide variations in any 
small sample, or to explain the extremely long period of time that it would 
ake any one area to experience all of the variations, or even a representative 
sample, of those which are possible. There are just more possible combi- 
ations for hitting a paying combination than there are on a common variety 
f conventional amusement machines, and as many statistical errors to be 
voided in drawing a conclusion on the frequency of outflow as there would be 
rom a Similarly small sample of plays on such a machine. 

Since none would draw conclusions from a small number of personal obser- 
ations of action of such a simple machine, what recommends doing so about 
sither rainfall or streamflow observations, either by the authors’ procedure, 
r by any other? Instead, let this profession be built on the processes by 
hich no single observation alters the conclusion by a determinable amount. 


J. L. H. PAULHUS, ! A. M. ASCE and J. F. MILLER.2—There certainly can 
e no argument against the contention that a good relation can be developed 
etween any two sets of randomly selected unrelated items if the items in each 
et are arrayed in like or opposite order. This fact is so well known that it 
oes not seem to deserve the lengthy treatment given it by Riggs and Benson. 
he effects of the ranking procedure on correlation had already been explained 
nly a short time before by a Dutch climatologist, C. Levert, in a discussion 

“Relation Between Point and Areal Rainfall”, published in the Bulletin of 
he American Meteorological Society, Vol. 38, pp. 618-620, Dec., 1957. How- 
iver, the subject has no bearing on the authors’ paper as correlation was not 
hvolved—the objectives of the study being the same as those of any other 
hethod of extrapolating streamflow records, namely, the approximation of 
ood-frequency data and not of individual annual peak discharges. Other 
omments in the discussion indicated that several important points in the 
riginal paper had been misinterpreted. 


Staff Hydrologist, U. S. Weather Bureau, Washington, D. C. 
Meteorologist, U. S. Weather Bureau, Washington, D. C. 
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The second sentence of the discussion implies that a relation was develop ed 
between synthesized and observed annual peak discharges and then used to ex= 
tend the record. The writers did not develop such a relation at any time. The 
observed and synthesized annual peak discharges were arrayed in like order ~ 
to facilitate comparison for determining whether frequency analyses of the 
two would yield similar results. When it was determined that frequency analy- 
ses of the synthesized peaks would yield results within 20% of those for the 
observed, the record of synthesized peaks was extended by means of rainfall 
data, a rainfall-runoff relation, and a unit hydrograph. At no time wasthe 
record of synthesized peaks extended by a relation between observed and 
synthesized peaks for a specific period. 

The objection to the ranking procedure and the suggestion that “the author i 
method might have been evaluated by applying it to the latter part of along 
flood record and testing the extended results against the known earlier part of 
the period” are very easily answered. Using the only two basins included in _ 
the original study that had records of 30 or more years, the writers have per= 
formed the suggested test while showing that the ranking procedure has no 
bearing on the results. 

In the top half of Table 1 are listed in chronological order the observed and 
synthesized annual peak discharges during the 16-year period, 1938-53, for 
the Whippany River at Morristown, N. J., (29.4 sq. mi.). The synthesized peal 
discharges were obtained from runoff computed by applying a rainfall-runof 
relation for the Anacostia River, about 175 miles SW of the project basin, to 
1938-53 storm-rainfall data from the 2-station network, Plainfield and 
Elizabeth, N. J., roughly 25 miles SE of the project basin. The 6-hour unit 
hydrograph used for converting the runoff to discharge was based on data 4 
subsequent to 1939, as was the testing and adjusting of the Anacostia River 
rainfall-runoff relation. Application of the Gumbel method of frequency analy 
sis to the observed and synthesized peak discharges for the 1938-53 period — 
yielded the results shown at the right in the top half of Table 1. The synthe- 
sized frequency data are well within the 20 per cent limitation, and the 
rainfall-runoff relation and unit hydrograph were therefore considered adapté 
ble to the project basin. The arrangement or order of the observed and 
synthesized peak discharges has absolutely no effect on the flood-frequency 
data. The discharges could be ranked in like or opposite order or scramblet 
in any fashion whatever, and the frequency data would remain the same. ' 

The lower part of Table 1 shows the observed and synthesized peak dis- 
charges for the 1922-37 period and their frequencies. The same rainfall- — 
runoff relation and unit hydrograph used in synthesizing the annual peak dis 
charges for the 1938-53 period were used to synthesize those for the 1922-3 
period, the only difference being that they were applied to the Plainfield- 
Elizabeth storm rainfall data for 1922-37 instead of 1938-53. Again, the 
frequency data for the synthesized annual peak discharges are well within 20 
per cent of those for the observed. 

Table 2 shows the same kind of test applied to synthesized annual peak di 
charges for the North Branch Raritan River near Far Hills, N. J., (26.2 sq. 
mi). In this case, a rainfall-runoff relation for the North Branch Raritan Ri 
near Raritan, N. J., which includes the project basin, was applied to Plainfie 
Elizabeth (about 30 miles ESE) storm-rainfall data for the 1938-53 period tc 
synthesize annual peak discharges for the same period. Data for developing 
the unit hydrograph and testing and adjusting the rainfall-runoff relation wel 
all from the period subsequent to 1939. The synthesized flood frequency dat 
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TABLE 1 


FLOOD FREQUENCIES FOR THE WHIPPANY RIVER AT MORRISTOWN, N. J. 


Flood Frequency 


Period Observed Synthesized 
(cfs) (cfs) 


Annual Peak Discharges 


ear Observed Synthesized 
(cfs) (cfs) 


772 
1259 
1582 
1990 
900 
1947 591 420 | 
1946 1020 860 i 
1945 460 690 
1944 640 1150 
1943 540 470 
1942 1220 1460 
1941 679 510 
1940 1500 830 
1939 529 470 
1938 1170 1990 
8 836 
a aa his S. D. = Standard Deviation 
lode 696 614 
SS SET 
1937 491 560 2 661 617 
me ey ‘wo | oe 
1935 
a 620 600 25 1484 1531 
933 1060 750 
932 485 430 
931 456 400 | 
930 466 420 
929 830 430 
928 1100 1320 
927 394 1140 
926 585 1130 
925 680 380 | 
924 830 650 
923 600 430 
922 700 480 | 
lean 704 665 : 
ed 299 333 
| 
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being within 20 per cent of the observed, the same rainfall-runoff relation and 
unit hydrograph were then used to synthesize annual peak discharges for the © 
1923-37 period from storm rainfall for that period. Again, the synthesized 
flood frequency data are within 20 per cent of the observed. 

Nowhere in the two tests just described would any rearrangement of the 
order of the observed or synthesized annual peak discharges within any one 
period have any effect on the computed flood frequencies. 

The results of the above tests were satisfactory, considering the short 
periods of record used for controls and testing —20-year periods should be the 
minimum and longer periods would be preferable. Generally speaking, this i 
type of test is not a fair one for evaluating the writers’ procedure because the 
off-basin rainfall network does not measure the basin rainfall storm by storm, 
With the rainfall network fairly well removed from the basin, it is impossible 
to distinguish between the errors of the synthesizing procedure and sampling — 
errors in the statistical sense. It would be very interesting to see a regional — 
procedure subjected to this same type of test. ; 

The best test of the writers’ procedure would be to use a basin with a long 
streamflow record, say 40 years or longer, and having a rainfall network for 
5 to 10 of those years. The basin storm-rainfall data and concurrent stream= 
flow record would be used for selecting or developing a suitable rainfall- si 
runoff relation and for constructing a unit hydrograph. The relation and unit 91 
hydrograph would then be applied solely to the storm-rainfall data of an off- ~ 
basin network to synthesize the annual peak discharges for the entire period © 
of the basin’s streamflow record. Comparison of the observed with the syn’ e 
sized flood frequencies would provide the best possible test of the procedure. 

Incidentally, Messrs. Riggs and Benson misinterpreted the purpose of Fig. 
6, which was to avoid the necessity of having several unit hydrographs for the 
same basin. The computed peak discharges were obtained by applying the se= 
lected rainfall-runoff relation and unit hydrograph to the basin rainfall actu- 
ally producing the observed peak discharges. The observed and computed 
peak discharges were then plotted against each other on a concurrent basis 
ea ae after being ranked, as should be obvious from the plotted points of 

ig. 6. 

In the discussion it was stated that an annual flood peak record of about 15 
years is required for applying the writers’ procedure. This is not so. The 
basins used in the study were restricted to those having at least 15 years of 
streamflow record so that the observed annual peak discharges would serve — 
as a sort of “yardstick” for evaluating the synthesized peaks, as in the above 
tests. Actually, concurrent measurements of streamflow and of rainfall from 
an adequate network over the problem basin for a reasonable number of mode 
ate and heavy storms would be the only requirements for application of the . 
technique suggested by the authors. The rainfall and streamflow data would 
determine the choice and adjustment of the rainfall-runoff relation to be use 
and permit the development of a unit hydrograph. One or two years of recor 
would suffice in many cases although longer records would, of course, be 
preferable. Furthermore, the observations need not be continuous but may — 
be made only when major storms or high discharge rates are expected. 

There is absolutely no reason why the synthesized and observed peak dis- 
charges should correspond on an annual basis—no more so than observed pe 
discharges for basins being used in developing a regional procedure. If the 
rainfall network were in the problem basin, agreement between synthesized 
and observed peak discharges might then be expected. The agreement would 
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TABLE 2 


FLOOD FREQUENCIES FOR THE NORTH BRANCH RARITAN RIVER NEAR FAR HILLS, N. J. 


Water Annual Peak Discharges Flood Frequency 
Year Observed Synthesized Return Observed Synthesized 
(cfs) (cfs) Period (cfs) (cfs) 


S, D 906 997 S. 0. = Standard Deviation 

Mode 1111 1127 

ees or ae eee ES RSS tet eee Pe bee Oo) ee 
1 1270 480 i 2 1109 1160 
tose 1860 1960 5 1747 1925 
1935 732 760 10 2170 2432 
1934 2630 2170 25 2703 3073 
1933 1480 1480 

1932 1000 970 

1931 1320 770 

1930 615 870 

1929 1040 710 

1928 1850 1790 

1927 857 1820 

1926 1230 2740 

1925 692 hho 

1924 556 1020 

1923 750 890 

Mean 1192 1258 

Ss. D 574 689 
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not be perfect because of inaccuracies in rainfall-runoff relations and unit 
hydrographs and the probably greater inaccuracies involved in estimating the 
true volume and the true time and geographical distribution of the rainfall 
from the network data. In the study described by the writers the rainfall net 
works were not generally located in the problem basins, because one of the 
primary objectives was to determine if the rainfall data from one network 
could be used to synthesize peak discharges for several basins so that the 
costs of applying the procedure on a per-basin basis would be thus reduced. 
The only limitation on the use of the rainfall data from a given network to 
synthesize the peak discharges for a particular basin is that both the network 
and the basin must have similar precipitation characteristics. Specifically, © 
they must have approximately the same rainfall-intensity frequencies and 
their intense rainfalls must occur at about the same time of year. It certain- 
ly is not required that rainfalls of relative magnitudes over the network and 
over the basin occur on the same dates. 

The reasoning of the writers was that over a long period the basin would 
have a rainfall record similar to that of the network provided both were ina 
region of similar precipitation characteristics. This is certainly a better 
basis for obtaining flood-frequency data than the regional methods mentioned — 
in the discussion. Users of regional methods have to assume that the hydro- 
logic and geologic characteristics, as well as the precipitation characteristics 
of the problem basin are the same as those of the basins used in developing 
the regional method. Of these, the precipitation characteristics are the most 
conservative and may therefore be interpolated with more confidence than th 
hydrologic and geologic characteristics for an ungaged basin. 

The writers’ procedure does make use of causally related data for the 
problem basin, namely, for adjusting or developing a rainfall-runoff relation 
and for constructing a unit hydrograph. If any method can be labeled as 
abandoning the use of hydrologic data causally related, it is the regional 
method of the usual type. Generally, the only independent parameter is the 
size of the drainage area. Even if the problem basin has a 5-year record of 
rainfall and streamflow measurements, the data are usually ignored in regio 
al analysis, 

Another advantage of the procedure described by the writers is that it can 
be adapted to the latest hydrologic characteristics of the basin and might, in 
some cases, yield more useful results than even the record of observed 
streamflow. For example, a small basin with 35 years of streamflow 
measurements has gradually undergone appreciable changes in land use duri 
the first 30 years of record. If the runoff characteristics were seriously af- 
fected, the streamflow for the earlier years would not be comparable to that 
for the last 5 years. A frequency analysis of this record would have little 
significance in design. Neither would the application of the usual regional 
method. However, by adapting a rainfall-runoff relation and constructing a 
unit hydrograph on the basis of the last 5 years of streamflow data, the 
authors’ procedure could reasonably be expected to yield more reliable resul 
than could be obtained from either the streamflow record or the regional 
method. 

The writers admit that there is a great deal of work involved in the appli- 
cation of the procedure, and it would be uneconomical to use it in the design ¢ 
isolated low-cost structures. However, only about two man-months are re- 
quired for an experienced person to make the storm-rainfall survey (50-year 
record), establish the time distribution of the rainfall, adapt a rainfall-runo 
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elation, construct a unit hydrograph, convert the rainfall into peak discharge, 
nd make the frequency analysis. The storm survey and time distribution of 
he storm-rainfall data require about one man-month, so there is a 50 per 
ent reduction in time required if these two items have already been es- 
ablished for an applicable network. There are certainly many structures for 
vhich the expenditure of one or two man-months in obtaining more reliable 
lesign data would be justified. 

In view of the comments in the discussion, it may appear strange that the 
yvriters look upon their procedure as an aid to, as well as a substitute for, 
-egional or generalized methods. The development of regional methods is 
yften hampered by the relatively short streamflow records available, especial- 
y for small basins. The writers’ procedure would provide 50-year or longer 
“ecords as a basis for such methods and thus extend their range of reliable 
pplication to include design values for longer return periods than is usually 
ossible. 
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FLOW CHARACTERISTICS ON THE OGEE SPILLWAY2 


Discussion by Donald P. Thayer 
Closure by Robert B. Jansen 


_ DONALD P. THAYER,! M. ASCE.—The author has done a service to the 
yrofession in bringing up and proposing a solution for one of the most signifi- 
sant problems confronting the hydraulic designer; i.e., the dissipation of 
energy by friction on an overflow spillway or steep chute. 

The writer is in general agreement with the author as he writes, “It is 
reasonable to assume that the energy loss on the spillway face is determined 
y a Simple natural law.” It is felt that the most logical approach would con- 
sist of making a theoretical analysis, taking into account the major factors af- 
ecting the phenomenon. Such an analysis should be of immediate use for pre- 
iminary design and subsequently should be a convenient framework within 
— to correlate experimental measurements. From such correlation there 
hould result emprical coefficients bringing the formula into close agreement 
vith experience. By “experimental measurements”, the implication intended 
S measurement of prototype behavior; the writer’s observation has been that, 
or quantitative evaluation of friction losses, conventional models, based upon 
‘roudian similitude, offer little help. More will be said of this later in this 
liscussion. 


Theoretical Analysis 


A theoretical analysis of the friction loss on an overflow spillway will be 
ttempted. The well-known Darcy formula will be used for evaluation of the 
riction head, 


2 
Nae roee . 
Cito et ee (1) 
where: f = the friction coefficient 


L = the length considered 
D = diameter of the pipe 
v2/2g = velocity head 


The choice of this formula is based upon two considerations: (a) it is di- 
aensionally homogeneous, and (b) it appears to be thoroughly documented by 
xperimental data. It should be noted that this formula is identical, mathe- 
1atically, with the much earlier Chezy formula; the respective constants ap- 
earing in a different form and differing by a constant ratio. For the purpose 


. Proc. Paper 1452, December, 1957, by Robert B. Jansen. 
. Asst. Div. Engr., Div. of Design and Constr., California Dept. of Water 
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of this analysis, however, this formula has two shortcomings: (a) it has been 
derived for and its use has been limited, almost exclusively, to closed circul 
conduits, and (b) its friction coefficient f is a function, not of the roughness 
alone, but of €/D, the relative roughness, where € = the “absolute roughness’ 
ideally a linear measure of the height of the surface irregularities, but pract 
cally an experimentally determined quantity. This formula is often applied to” 
non-circular closed conduits by letting D = 4R, where R is the hydraulic 
radius, and may, with perhaps somewhat less validity be applied to open 
conduits. (1 

The channels considered in this discussion will generally be very wide 
compared to their depth so that without significant error we may say that the 
depth of flow, perpendicular to the bottom, is equal to the hydraulic radius, F 
With this in mind, the Darcy formula may be written 


2 
he = L e bho 
TUR ° 2g 
If the hydraulic slope is defined as s = hg/L, then the formula becomes 
2 
s_= ag = Vv 
SR ge 


In Fig. 1 consider an elemental prism of the flowing water of weight W, 
shown cross-hatched, moving down the spillway slope, A-B, from position M 
to position N. Let the line C-D be the trace of the centroid of the mass as it 
moves. Measured at this centroidal line, the hydrostatic pressure potential 
energy of the elemental prism of water is zero and its total energy may be © 
written as the sum of its potential and kinetic energies, 


= Wve 
2g 
where: y = distance of the centroid above some datum, as yet undefined, 
with positive direction upward. 
As the mass moves down the slope the change in energy may be obtained by 
differentiation with respect to L. 


gE . Wy dv et 
att eee Moree 


E + Wy (4 


For a frictionless condition, the right-hand member of Eq. (5) would be 
equal to zero; with friction it is equal to minus the hydraulic friction slope 
times the weight of the mass, 


dk 


aT = awa 


From geometry of Fig. 1, 


aL = -dy/sina 


where: o = angle of inclination of line C-D, at the point considered, from 
the horizontal. 


From continuity of flow, 
R= a/v 


where: q = specific discharge passing down the spillway. 
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FIG.1- OVERFLOW SPILLWAY 


serting relationships, Eq. (3), (6), (7), and (8) in Eq. (5) and multiplying out 
L, gives 


i mS RS Pelee Ee 


ividing by W and rearranging 


d fv3 
a ed ay (10) 
et 
fv3 
B3 = PORES (11) 
hen 
f 2/3 was 2/3 2 
Be = ae : [ a kA (12) 
ot 
es 
BL [x ] (13) 


yte that (qg)1/3 = vo, the critical velocity, then 
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Note also, that Ve2/g = De, the depth at critical velocity, and substitute this 
relationship, together with the first differential of Eq. (14) for the left-hand 
member of Eq. (10) to obtain 


BdB = d: 
fee aaa (15) 


Integrating 
BAB. sates ' 
5 | Es = Fis FC (16 


The denominator under the integral sign has the factors 


([p - 2].[B+ (2 +4 v3)].[+ G-4 v3) 


These factors become the denominators of partial fractions formed in the 
usual manner. These partial fractions upon integration and the combination 
of conjugate complex quantities give as a value of the left-hand member of 
Eq. (16) 


2 B-2)¢ a . 
foc Bethy - eve H]= 4) Ge 
Then, using this notation, Eq. (16) may be written . 
c 
Eq. (14) may be written 
2a 
roe (19) 


where: hy = velocity head at any supercritical velocity 
he = velocity head at critical velocity 


It should be noted that the function (B), Eq. (17), is inherently negative, vary- 
ing from -2,.4179 when B = 0, to negative infinity when B = 2. 
Eq. (18) and (19) then give dimensionless relationships between y, measure 
in terms of Dc, and hy, measured in terms of he. These equations have been 
evaluated for several different values of F and the results plotted in Fig. 3. — 
In all cases the value of C in Eq. (18) has been so chosen that the curve passe 
through the point (hy/he = 1, y/De = +0.5) thus placing the datum as shown in 
Fig. 1, at the elevation of the crest where critical velocity is attained. It is 
recognized that in many cases critical velocity and depth are not attained at 
the fixed crest. In this case it is a simple matter to compute a critical depth 
at which the given specific discharge would pass over the given crest or, mor 
precisely, to compute a fictitious crest elevation over which the given specifi 
discharge from the given pool elevation would pass at critical velocity. It is 
believed that in many cases the former method is sufficiently accurate. In 


this connection the following relationship may be useful; the equivalent critica 
depth, 


Di = .3U-meGa oe 


SCE DISCUSSION 1856-29 


Values of F 
0.20 


os 


Values of VLolebDrook T 


0.010 0.015 0.020 
Values of Manning n 


FIG.2-DIAGRAM RELATING n, R, f, «, AND F 


where: H = head from the given crest to stagnation pool level 
C = coefficient of crest, as inq = C H3/2 


Fig. 3 has been drawn to give a rapid graphical method of obtaining the re- 
tionship between Manning n, the Colebrook function for f in the Darcy formu- 
, and F as used in this discussion. The use of this diagram should be self- 
rident, except for the angle, a, to be used. Theoretically this is the angle 
e€ median line of the overflowing nappe makes with the horizontal; actually, 
the lower part of the spillway, where the greatest part of the friction loss 
curs, this line is nearly parallel to the spillway and the latter angle should 

satisfactory. The depth, R, varies continuously throughout the spillway so 
at a value approximating the mean in the region where the greater part of 
e friction loss occurs would probably be the best to use. 


Comparisons and Application 


For purpose of comparison, three of the points representing velocity 
easurements at Madden Dam as reported by Mr. Randolph(2) have been plot- 
don Fig..3. The correlation is not particularly good, although the trend is 
the right direction. It must be pointed out that the plotting may be inaccu- 
te since the values were scaled from Fig. 17 in Mr. Randolph’s paper, the 
“iter being unable to find numerical values. Two points from the model 
idy of the Madden Dam, also scaled from this same Fig. 17, are also plot- 
1. These points fall on a much smaller value of F than any of those for the 
ototype, although the roughness of the model and prototype, as reported by 
r. Randolph were in approximately the correct Froudian relationship 
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FIG.3— FRICTION HEAD LOSS ON SPILLWAYS 
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issuming, of course, that the Manning formula is applicable to both). This 
mphasizes the statement made earlier regarding the inapplicability of model 
esults where friction is concerned. 

Curves drawn according to the author’s Eq. (2) are also plotted on Fig. 3; 
ne with K = 0.10 and the other with K = 0.125. These curves cut sharply 
cross the F curves indicating considerably discrepancy. It might be pointed 
ut that the author’s curves are not asymptotic to a vertical line, indicating 
o terminal velocity. For practical purposes, with small relative values of y, 
1is should not be considered a serious defect; the point is more theoretical 
lan practical. 

The analysis made in this discussion disregards two factors which might 
ave a considerable affect on the spillway velocities: 


- (a) The friction of the atmospheric air on the overflowing sheet of water. 

_ (b) The effect of entrained air on the friction coefficient. 

s the author points out, engineers are not agreed on this point and there are 
nflicting data. Viparelli(3 states that for calculation of velocity the value 
n appropriate for the channel surface should be used. He cites earlicr ex- 
riments by Mongiardini with flow of an air-water mixture through a closed 
pe indicating that for conditions of complete turbulence (Reynolds’ number 
eater than 50,000) the friction coefficient is the same as for water alone. 
should be noted that these experiments used a rather small pipe, 20.8 mm 

iameter. On the other hand, Straub and Lamb(4) show that velocities in an 

xperimental chute computed from the “normal” value of n are considerably 

- than those measured, 

It is believed that more measurements of prototype behavior are needed to 
larify this point. : 


Future Investigations 


No other prototype measurements are available to the writer to compare 
ith his analysis. It would be very helpful if anyone having such data would 
ot it on Fig. 3 or forward the data to the writer for such purpose. 

The lack of reported prototype measurements applicable to this subject is 
ther unfortunate, and it is hoped that more data comparable to that collected 
Mr. Randolph will be collected and reported. It is encouraging to note that 

e north spillway bay of the Detroit Dam has been arranged to permit such 
easurements, and it is hoped that data from this source will be made availa- 
e in the near future. 
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ROBERT B. JANSEN,! A. M. ASCE.—The subject of head loss on the spill 
way face is deserving of more attention. Further experimentation and analy- 
sis should be directed toward determination of the mechanics of this process” 
of energy dissipation. The rate of energy variation surely must be determine 
by the relationship of a relatively small number of factors. In the case of 4 n 
ogee spillway, or any straight spillway with constant slope and roughness, the 
pattern of energy change obviously is established primarily by conditions at 
the crest. 

For a given spillway slope and roughness, the head loss should be a functii 
of the distance of fall and the rate of discharge. It was on this basis that Eq. 
(3) was suggested. Application of this formula, with K equal to 0.125, gives ~ 
results fairly close to the prototype observations on the spillway of the Madde 
Dam. Since the face of this spillway was reported to be rougher than usual, © 
it was estimated earlier that the value of K for a typical concrete spillway — 
might be closer to 0.10. Recent developments have shown this assumption to 
be questionable. 

The newly published work of Bradley and Peterka(1) contributes to an 
understanding of this phenomenon. Fig. 15 of their series of valuable papers 
shows the results of prototype observations on the spillways of Shasta and 
Grand Coulee Dams. Also presented are curves, drawn from the considerabl 
experience of the authors, which enable the determination of velocity and dept 
of flow “with a fair degree of accuracy, for chutes having slopes of 0.8:1 or 
steeper, where computation is a difficult and arduous procedure.” The proto 
type tests cover a range of crest heads from 1.6 feet to 13.6 feet with height 
of fall as great as 220 feet. These important measurements provide oppor- 
tunity for comparison with the Madden Dam prototype data, and for further 
analysis of basic relationships. 

The Shasta and Grand Coulee spillways have a face slope of 0.8:1, as com- 
pared with the Madden Dam inclination of 0.75:1. All three have drum-gated 
crests; Shasta with three gates 110 feet by 28 feet, Grand Coulée with eleven 
gates 135 feet by 28 feet, and Madden with four gates 100 feet by 18 feet. 

Since the reported prototype measurements on these three spillways were 
limited to crest heads not greater than 13.6 feet, there is need for more info! 
mation on flow at higher heads. No other prototype tests are known to the 
writer. However, flow at high crest heads is ina range where, at least for — 
relatively low heights of fall, model and prototype results should be in fair 
agreement. That is, discrepancies attributable to roughness and air en- | 
trainment should not enter seriously. It is considered acceptable, therefore, 
to use the results of model experimentation for comparison in this range. 
this purpose, reference is made to two works of the Waterways Experiment : 
Station, Corps of Engineers, Vicksburg, Mississippi. These are Bulletin No. 
37(2) and the model studies for Center Hill pam. Energy losses for large 


crest heads can be calculated by using the nappe profiles presented in these 
papers. 


1. Supervising Hydr. Engr., Dept. of Water Resources, State of California, 
Sacramento, Calif. 
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The data from the above various prototype and model tests have been ana- 
zed by the writer in an attempt to establish a generally reliable formula for 
»mputing the head loss on the face of an ogee spillway. Logarithmic graphs 
ere used to establish exponents. It was discovered that the head loss varied 
rectly as about the 4/3 power of the height of fall, and inversely as the 1/3 
wer of the critical depth at the crest. (These are representative values, 
ith emphasis placed on the apparently more consistent Shasta Dam test data. 
he exponents derived from the Maddem Dam prototype results were 3/2 and 
/2, as reflected in Eq. (3)). These findings lead to the following relation: 


=JZ 4/3 
Del/3 


herein J is a coefficient, Z is the height of fall from the reservoir level and 
; equal ss sum of P and H, and Dg is the critical depth at the crest and is 
C4/3H 
jual to ei/3 * 
_A value of J = 0.15 yields results which agree well with the prototype and 
odel records now available. As a matter of interest, a value of K between 
125 and 0.150 in Eq. (3) would give computed losses in good agreement with 
ese same test data. It appears therefore that the suggested reduction of K 
om 0.125 to 0.100 because of the unusual roughness of the Madden Dam spill- 
y is not justified. A more proper conclusion from the Madden Dam analy- 
Ss only, in view of present knowledge, would involve a coefficient of 0.125 in- 
ead of 0.100 in Eqs. (2) and (3). 

Although Eq. (10) is supported by only a meager experimental foundation, 
has merits which should be noted. It is simple, easy to apply, and gives 
parently dependable results. Use of the term D, includes the effects of both 
and H, which were treated individually in Eq. (3). On the basis of the spill- 
y observations reported to date, Eq. (10) appears to be in a more acceptable 
rm than Eq. (3). It is probably as consistent a criterion as can be offered 
Ww. 

Mention should be made of the influences of the roughness and the slope of 
ie spillway face. The former should not be cause for great concern. In the 
nstruction of a modern concrete spillway, enough attention is usually given 
concrete quality and formwork to ensure a uniform roughness. The effect 
the slope of the spillway face cannot be determined on the basis of the three 
ototypes which have been studied. There is not enough difference in their 
e angles. However, most ogee spillways will fall into a range not much 
eater than the one in reference (0.75:1 to 0.80:1). To extrapolate the known 
ta into the range of flatter slopes, it might be assumed that the energy loss 
ries directly with the length of travel along the spillway face, or inversely 
the sine of the angle between the face and the horizontal. If the term J in 


aL (10) 


. (10) is assumed equivalent to nee then the head loss is given by the 


Owing expression: 


h, =_J z4/3 
ms sina ‘ 1/3 (11) 
lc 


which j is a coefficient and « is the angle of inclination of the spillway face 
‘h the horizontal. 
Eq. (11) is only speculative. The subject of loss on slopes flatter than the 


al gravity-dam slopes is beyond the scope of this paper. 


1856-34 HY 6 November, 195 


It should be noted too that the prototype measurements on which this studj 
is based are limited to values of Z generally less than 200 feet. However, 
this does not necessarily place a limitation on the applicability of the resul 
In using Eq. (10), the calculations of energy loss should be performed pro- 
gressively at intervals down the spillway face until the terminal velocity is” 
reached. This is the point at which the increment in energy loss becomes 
equal to the vertical increment of fall. Since the velocity head is assumed 
constant below this level, the incremental head loss beyond the arrival at 
terminal velocity is taken as equivalent to the increment of fall. This means 
that Eq. (10) need be applied only in the range where the rate of energy loss 
is varying; i.e., only down to the point where the terminal velocity has been 
attained. 1 

Mr. Thayer has made a remarkable contribution to the literature in his 
theoretical analysis based on the Darcy formula. Essentially, it appears thi 
he has independently accomplished the same result as achieved by Bresse n 
his method.(4) Bresse performed his integration assuming a very wide rec- 
tangular channel and adopting the Chezy formula. Mr. Thayer has used the 
Same assumption, and has employed a formula which, as he states, is mathe 
matically identical with the Chezy relation. 

Much has been written about the limitations of the Darcy formula as appli 
in the way that Mr. Thayer suggests. Admittedly arbitrary premises are in 
volved. Open channels are usually rougher and larger than pipes, with high 
values of the Reynolds number, corresponding practically to complete turbu: 
lence. In this case, the Darcy coefficient “{” would be practically independel 
of the Reynolds criterion by which it is commonly determined. The belief hi 
prevailed that some other index, such as the Froude number, should be em= 
ployed for flow with a free water surface. 

It would be hard to imagine more dissimilar conduits than an ogee spilly 2 
and a small pipe such as that for which the Darcy formula was developed. 4 
attempt to use the Darcy relation for analysis of a steep spillway would be 
accompanied understandably by serious doubts. However, comparison of he 
Thayer results with data from other sources reveals a relatively close 
agains oe for flow with low crest heads. Discrepancies appear at higher — 

eads. ; 
Although Mr. Thayer has performed an analysis closely paralleling the — 
work of Bresse and has suggested a criterion basically the same as the mel 
of Bresse, the new approach has advantages over the earlier work. The 
Thayer graphs are presented in such a convenient form as to make compu- 
tation by this means probably more attractive than by the Bresse tables. I 
Thayer, in his integration of incremental energy losses defined by the Dare 
formula, makes unnecessary much of the labor that would be involved in the 
common step method using the same formula. This is especially true if the 
energy relation at only a single point on the spillway face is concerned. Ift 
cumulative loss is to be determined at several points, there is less advanta 
For practical purposes, the result is the same whether integration is ac- 
complished by the calculus or by small finite steps. 

The two graphs of Mr. Thayer provide a rapid means of estimating the — 
energy loss on the spillway face. In using his curves, it is necessary first 
estimate the dimension R and to enter the first graph with this trial value t 
find the factor F. Mr. Thayer states that, “the depth, R, varies continuoush 
throughout the spillway so that a value approximating the mean in the regio 
where the greater part of the friction loss occurs would probably be the bes 
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ise”. It would first appear that this assumption would inject a considerable 
asure of approximation into the Thayer method. However, the writer has 
nd that the suggestion has merit. Also, the extent of trial involved in se- 
ting the value F is not an important handicap. The convenience of the first 
iph makes the trial a simple matter. 
When the sheet of flow is thin, the graphical selection of the proper value 
the Thayer factor F becomes rather difficult. In this range, F varies ap- 
-ciably with small changes in depth. The term F can be expected to vary 
much as 100 per cent on a typical ogee spillway. This is well demonstrat- 
by the first graph of Mr. Thayer. 
For large heads on the spillway crest, in general the second Thayer graph 
(00 Small in scale to yield sufficient accuracy. This is not an unusual 
advantage of a graph prepared for publication. In practice, the plot could 
made to a larger scale, or the basic data could be placed in tabular form. 
wever, it is not certain that the method of Mr. Thayer will give acceptable 
sults when the ratio y/D, is small. More experimental verification is need- 
in this region. 
able 1 is a compilation of data from the several sources which have been 
tioned. The tabulation shows that the energy losses computed from the 
observations on the Shasta spillway prototype are generally illogical 
higher magnitudes of the dimension Z. For example, with a crest head of 
feet the listed energy losses from actual measurements on the prototype 
40, 59, 77, 92, 103, and 108 feet corresponding to heights of 70, 90, 110, 
, 150, and 170 feet, respectively. The increments of energy loss are 19, 
15, 11, and 5 feet, in the same order. Logic dictates that, since the ve- 
ty will continuously increase until the terminal velocity is attained, there 
uld be no decrease in the increments of energy loss in proceeding down 
spillway face. The prototype records presented by Bradley and Peterka 
Shasta and Grand Coulee spillways manifest a general defect in this re- 
ct. In order to eliminate this inconsistency, the writer had adjusted the 
puted losses for the Shasta spillway by accepting the value of hy, corre- 
nding to the lowest recorded value of Z for each crest head, on the as- 
ption that the accuracy of measurement of flows on the prototype would 
ease as the height of fall increased. In making the adjustment, the rate 
nergy loss was not allowed to decrease. In the case of the Grand Coulee 
lway observations, since the initial value of hy, for a crest head of five 
appears obviously too low, the first listed energy loss has been made to 
orm to the Shasta prototype data. The other adjusted losses shown for the 
ot crest head fit the loss differentials of the original Grand Coulee 
surements, except that again no decrease in the loss rate has been permit- 
For a crest head of 13.6 feet, the measured energy loss for the lowest 
e of Z has been accepted, although it may be slightly high, and an incre- 
tal energy loss of ten feet has been adopted from examination of the 
ern of the Shasta columns. 
he tabulated energy losses on the Madden Dam spillway were taken direct- 
om Fig. 17 of the paper by Mr. Randolph. A general nappe profile (Fig. 
Bulletin No. 37 of the Waterways Experiment Station has provided an 
mate of losses at crest heads of 20 and 30 feet. Similar data from the 
el studies of Center Hill Dam, performed at the same laboratory, give 
‘oximate losses for flow with 44 feet of head on the crest. These measure- 
s were taken from the original model, which liad a face slope of 0.61:1. 
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Presented also in Table 1 are energy losses determined on the basis of 
curves drawn in Fig. 15 of the Bradley-Peterka series. This graph shows — 
that the authors themselves have adjusted for the obvious defects in the pra 
type observations on the Shasta and Grand Coulee spillways. 

The last two columns in Table 1 are based on the Thayer method with — 
Darcy coefficients corresponding to Manning n = 0.014, and on the Jansen EK 
(10) with the factor J equal to 0.15. It is believed that Table 1 demonstrate: 
that Eq. (10) offered by the writer gives a fairly reliable estimate of the ene 
gy loss on the face of an ogee spillway, on the basis of comparison with the” 
prototype and model observations which have been discussed. As has been 
stated, the Thayer method when applied to large crest heads is at variance ~ 
with the other listed references. However, information presently on record 
does not permit a definite conclusion in this regard. More observations on 
spillway prototypes are sorely needed. 

One avenue of investigation which can be explored, at least as an appre $ 
mate means of correlation, is application of the law of hydraulic model si- 
militude. For example, disregarding lack of similarity in roughness and if 
air entrainment, flow on a spillway 200 feet high with a crest head of 20 fee 
should involve approximately twice the energy loss as on a spillway 100 fee 
high with a ten-foot head on the crest. 

The superposition of the data of the writer on the second Thayer graph 
would have resulted in a more logical curve if attainment of terminal veloci 
had been taken into account. As has been stated earlier, the equation defini 
the energy loss should be applied only down to the point where the terminal 
velocity has been reached. Eq. (10) is recommended for this purpose, until 
more comprehensive formula can be written. It is recognized that this 
equation is the result of a mere scratching of the surface of this import 
subject. Certainly a more satisfactory criterion ultimately will be discove 
In the meantime, the procedures herein described may serve as a guide. 

The profession is fortunate to have received the outstanding works of 
Messrs. Randolph, Bradley and Peterka, Their important contributions to | 
literature have laid the foundation for a broader understanding of spillway 
hydraulics. In closing this discussion the writer wishes to record appreci- 
ation of the notable treatise by Mr. Thayer. The time which he has given t 
the clarification of the subject matter herein has greatly enriched the papél 
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TABLE 1. - ENERGY LOSS ON SPILLWAY FACE 


Loss in Energy Head, hy, in Feet 
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DISCHARGE CHARACTERISTICS OF RECTANGULAR THIN-PLATE WEIRS2 
Discussion by M. R. Carstens 


M. R. CARSTENS,! A. M. ASCE.—The authors are to be commended for 

1e study of rectangular-notch weirs. 

The authors have attempted to eliminate the effects of viscosity and surface 
snsion from experimentally determined values of the coefficient discharge of 
thin-plate weir by the simple expedient of an effective head and width. The 
uccessful elimination of these two effects should result in the coefficient of 
ischarge of potential (irrotational) flow without surface tension. Thus, the 
ae a7) for the coefficient of discharge for the full-width weir, Eqs. (25), 

6), (27), (28), and (29), should be identical for these five series of tests. The 
mall differences in these equations could be attributed to the following 
auses: 


(a) The use of the effective head does not completely eliminate the effect of 
viscosity and surface tension and the differences in the equations repre- 
sent differences between the five series of experiments in this regard; 

(b) The geometry in each of the five experiments was sufficiently different 
to account for the small differences; and 

(c) The differences are merely an indication of the order of accuracy with 
which the coefficient of discharge can be experimentally determined. 


Since the maximum deviation of the data from each experiment (Figs. 3, 4, 
, and 9) is greater than the deviations between the equations, the writer be- 
eves that (c), above, is the explanation for the small differences in Eqs. (25) 
(29), inclusive. Therefore, the writer would have combined all the data (re- 
uced to potential flow) shown in Figs. 3, 4, 5, and 9. 

Despite the tenuous derivation of Eq. (15) the solution for the coefficient of 
ischarge is amazingly good. Eq. (15) for a full-width weir is 


Yn Ve 
q= Co 2%25 { [+a] - |] | (16a) 
9 


his equation can be placed in the usual form of the discharge equation by 
.ctoring the quantity, h3/2, from the parentheses. 


Yo 4 Wp 
q= © 2V2q ie = | | he (15b) 
&} 2gh 29 


- Proc. Paper 1453, December, 1957, by C. E. Kindsvater and R. W. Carter. 
. Prof., Civil Eng., Georgia Institute of Technology, Atlanta, Ga. 
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The quantity, V12/2gh, can be transformed as follows: 


q = V,(h + P) 


and 
4 
q = (4 Cy 129 2 
3 
from which 
2 
VES Ci ( bi ) 
29h 9 h/P +1 


Substituting into Eq. (15b) 


3, 2 3h 
Cee Oe le oF (GREY + Pia ear ane tate | 
9 h/P +i 9 {Wp +1 (1 


If the theoretical values of C, for a two-dimensional orifice are utiliz 
Eq. (15c), the equation can be solved by successive approximation for Cg 
a function of h/P. In the following table, the theoretical values of Cg are ta 
lated for comparison purposes with the best available experimentally de- — 
termined values, that is, Eqs. (25), (26), (27), (28), and (29). In addition, 
experimentally determined values of Kandaswamy and Rouse are tabulated 
order to compare the theoretical solution with an experimental solution ob 
tained at larger values of h/P. 

The theoretical solution for the coefficient of discharge is within a few > 
cent of the various experimental solutions throughout the entire range of the 
weir, that is, 0< h/P < 10. The writer cannot agree with the authors’ 
statement, 


“The result is a complex equation which must be modified to make it 
practicable and adjusted to make it agree with a particular set of ex- 
perimental data. The fact that a large number of different empirical 
formulas have been evolved by this procedure is evidence that the 
“theoretical” equation is inadequate. The obvious reason for the 
inadequacy is the fact that the basic assumptions are inaccurate.” 


Rather it appears to the writer, that too many people write empirical 
equations. 
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AIR BINDING IN LARGE PIPELINES FLOWING UNDER VACUUM2 


Closure by R. T. Richards 


R. T. RICHARDS,! A. M. ASCE.—The author agrees with Messrs. 
Templeton and Stelson that air evacuation equipment may be required at verti- 
al drops. This is particularly true of low velocity flow. The equations they 
lave submitted give an interesting comparison of the relative head losses at 
vertical and sloping drop. However, the author’s observations on large di- 
eter pipelines indicate that at vertical drops preceeded by relatively short 
adius bends and at water velocities of 6 to 8 feet per second, turbulence ap- 
ears to reduce air collection to a point where it causes no serious head loss. 
ch bends may be found, for example, at the majority of steam power plant 

n condenser discharges. While a few of these condensers become air 
und unless continuously evacuated, the majority are self priming. Hy- 
raulic efficiency here is of great importance. A sloping pipe from the same 
ondenser discharge would probably be air bound. 

Messrs. Templeton and Stelson close their discussion by appropriately 
mmenting that air binding can be a serious problem in conduits flowing 

er pressure as well as under the vacuum conditions discussed in the paper. 
he author’s firm has recently investigated head loss problems in two large 
iameter sloping power tunnels flowing several miles under pressure. The 
liection of air along the tunnel roof appears to have been a major contribut- 
factor in the observed excess head losses. Air relief taps similar to those 
scribed in this paper have been provided on several large diameter sloping 
essure pipes designed by the author; however, to date the test data from 
ese pipes is not sufficient to clearly establish the effect of air binding on 
ese particular installations. If continuous air removal is required it can 
ry easily be accomplished by automatic air relief valves. Removal of air 
om a pipe flowing under vacuum is far more of a problem, especially at 
peline locations where electric power is not readily available for air 
acuation equipment. 


Proc. Paper 1454, December, 1957, by R. T. Richards. 
Civ. Engr., Ebasco Services Inc., New York, N. Y. 
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TURBULENCE CHARACTERISTICS OF THE HYDRAULIC JUMPa 


Discussions by A. J. Peterka and J. N. Bradley, James M. Robertson, 
Donald R. F. Harleman, Philip G. Hubbard 


A. J. PETERKA,! M. ASCE and J. N. BRADLEY,2 M. ASCE.—The authors 
to be complimented for accomplishing an investigation which is difficult 
h in concept and in performance. Their study represents the ultimate ap- 
ach; to explain the workings of the hydraulic jump by means of its internal 
hanics. Although the study is academic in nature, there are certain as- 
ts of the results which help the practicing hydraulic engineer to rationalize 
use of stilling basins shorter than the hydraulic jump iength. 
t is common practice to make the stilling basin floor (usually concrete) 
rter than the true length of the hydraulic jump. Hydraulic model tests and 
otype verification tests have shown that, within limits, this is a desirable 
edure since little in the way of performance is lost and the cost of the 
ing basin is reduced appreciably. In analyzing the action in the jump, the 
ors have shown why it is possible to reduce the apron length to less than 
jump length. The drawings in Fig. 13 indicate that most of the losses oc- 
upstream from the end of the roller and that the downstream end of the 
is relatively ineffective in accomplishing hydraulic losses. However, 
velocity distribution curves in Fig. 6 show that a redistribution of velocity 
rs after the major part of the losses have occurred. For this reason, it 
i possible to reduce the apron length as much as appears possible in Fig. 
For example, in Fig. 13, for Fj = 2, it might be concluded that the apron 
extend only about 3dg downstream from the start of the jump. Fig. 6 for 
2, however, shows that at jump length 3d9 there is a major redistribution 
locity yet to be accomplished. If the apron were ended at this point, the 
bottom velocity might scour an unprotected riverbed. An end sill on the 
bn would help to redistribute the bottom velocity. Even at the end of the 
, 4.5d9, as indicated on Fig. 4, a slight redistribution of velocity is still 
expected according to Fig. 6. 
tests reported by the writers? the jump was found to be imperfectly 
1ed for a Froude number of 2, (dg = 2.1dj). The total energy loss was 
about 7 per cent in a jump measured to be 4.5dg long. Thus, the length 
ie jump was useful more to redistribute the velocities than to accomplish 


roc. Paper 1528, February, 1958, by Hunter Rouse, T. T. Siao, and 

. Nagaratnam. 

ydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 

ydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly of 
. S. Bureau of Reclamation, Denver, Colo. 

radley, J. N. and Peterka, A. J., “Hydraulic Design of Stilling Basins”, 
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losses. Nevertheless, the apron length could not be reduced appreciably fre 
the 4.5d9 given as the length of jump. 

For a Froude number of 4, the writers found the hydraulic jump to have 
better form since dg = 4.5dy. The energy loss was about 40 per cent. The 
incoming jet showed some instability, tending to flutter from top to bottom ¢ 
the jump. Surface waves downstream from the jump were more of a proble 
than scour tendencies. It was found that the apron length could be reduced 
considerably, however, without changing the jump action. With the addition 
chute blocks and an end sill, the apron length was reduced to about 3.6d9. F 
a Froude number of 4, the authors show in Figs. 13 and 6 that most of the 
losses have already occurred at 3.6d9, and only a slight redistribution of 
bottom velocity will occur downstream from 3.6dp9. 

For a Froude number of 6, the writers found the hydraulic jump to have 
ideal form, dg = 8d, and the energy loss in the jump was about 55 per cent. 
The apron length was reduced to 4d9 when cute blocks and an end sill were” 
added to the apron; Figs. 6 and 13 help to justify this length reduction. 

It is unfortunate that the authors’ tests for F, = 8 could not be complete 
since these data would probably have shown that the roller occupied more 
than half the jump length as stated in the conclusions. It was found in the 
writers’ tests that as the Froude number increased, the roller size increas 
until the roller length became the controlling factor in determining the jum 
length. The roller length approached the jump length for Fy, = 8, dg = 10.7% 
and the energy loss in the jump was about 65 per cent. A discharge of 28 
second-feet in a flume 4 feet wide produced a jump which was very turbule! 
and rough and the jump itself surged because the bottom jet did not carry — 
through at all times for the full length of the jump. Bottom jet action was 
tremely erratic, sometimes penetrating the jump for the full length and at 
other times, for only a short distance. At times, the jet appeared to jerk - 
horizontally and flutter vertically. 

Also, the large difference in conjugate depths resulted in the jump havir 
a steep profile. Slugs of water intermittently rolled down the steep face ar 
fell into the high velocity jet. A surface surge was set up each time a slug” 
made contact with the jet. This was strictly a gravitational effect which an 
air model may not duplicate. Actually, flow down the face is common to all 
jumps, but it is usually steady. These findings probably help to explain the 
difficulty found by the authors in adjusting their air model to obtain good di 
for Fj = 8. The writers, therefore, suggest that perhaps the agreement be 
tween the air model and water model at F; = 8 was better than suspected DB 
the authors. In the writers’ opinion, it would be impossible with presently 
known equipment to make accurate measurements in a water jump of Fy =! 
having appreciable physical size. In addition, if the air model boundary we 
set to the jump length, curve of Fig. 4, it is possible that the do test length 
was too short. In the writers’.tests, the jump length was found to be consi 
ably longer for the higher Froude numbers than shown in Fig. 4, It is not 
clear in the description of the water model tests whether the length ratio 
Lj/dg was actually measured in the water model or whether the Bakhmetel 
Matzke curve was accepted as the true length of jump. In the air model, 
shape of the boundary (the jump profile in a water model) could greatly inf 
ence the action within the jump. The jump length given in Fig. 4 for Fy = 
is 4,2d9; the writers’ tests showed 6.1d9. The reason for the difference is 
believed to be due primarily to the greater thickness of incoming jet and V 
of test flume used in the writers’ tests. Excessive friction in too small 
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iodels can lead to erroneous results. The authors found it necessary to 
iden their air model to decrease the wall effects noted in the 1-foot duct. 
his should serve as a warning to other experimenters who are tempted to 
3€ very narrow test facilities when wider sections could be used. 


JAMES M. ROBERTSON, ! M. ASCE.—The study reported in this paper is 
at only of great interest because of the light which it sheds on that most in- 
‘iguing hydraulic phenomenon—the hydraulic jump—but also because it repre- 
ents a thorough study of a separated flow which reattached. Many flow situ- 
‘ions, besides the jump, occur with a separated zone of flow and a roller- 
ke action of the separated flow on the main stream. 

The suitability of air flow in the carefully formed diverging channel as a 
odel of the hydraulic jump is justified by the authors on the basis of expedi- 
icy and verification of the mean flow pattern. With no intention of deprecat- 
ig the value and significance of the study, the writer will point out certain 
spects of the flow model wherein it cannot have been similar to a jump. Two 
atures of the flow model differed from those of the hydraulic jump; one of 
ese was inherent in the type of model chosen. The presence of a curved 

id boundary rather than the upper free surface of the jump certainly caused 

me modification in the nature of the flow while the initial mean and turbu- 
nce velocity conditions for the model jump differed from those normally ex- 
cted in open channel flow. 

The conditions of flow at the initial station of most hydraulic jumps must 
those of fully developed turbulent open channel flow. The probable vari- 
ons in mean velocity, turbulent shear stress, and turbulence intensity for 

s situation are sketched in Fig. A. In contrast, the initial conditions tested 
the air flow analogy involved relatively thin boundary layers at the upper 

d lower surfaces with an intervening core of essentially shearless, turbu- 
ce—free flow. In the Fj = 2 case for instance, the shear stress curves of 
. 7 indicate that the boundary layer thickness along the bottom was at most 
e-quarter of the initial flow “depth” even at relative distances x/dg of at 

st 2. Since turbulent boundary-layer behavior in adverse pressure gradi- 
ts.is sensitive to initial flow conditions, the actual flow following this con- 
ion must have differed considerably from that in a proper hydraulic jump. 
this respect the condition studied would appear to conform closely to a 

p following a sluice gate—a situation still of considerable interest. 

| Since the pressure increases as the flow proceeds in the x-direction (cf. 

. 5) the writer would expect to note evidence of separation along the “bed” 
jthe channel. That none is in evidence might appear to be due to the ex- 
ence of an abnormally thin boundary layer along this surface. Further con-— 
fleration indicates that this explanation does not appear to hold water. The 
telihood of turbulent boundary-layer separation occurring has been shown 
to depend upon the boundary-layer Reynolds number at the start of the 
ssure gradient and upon the rate of pressure rise. The paper does not 

€ sufficient information to properly determine the initial conditions but it 
id appear that at the initial station the momentum thickness Reynolds 

aber probably lay between 103 and 104. If so, it is expected (1) that the 
imum mean-flow velocity could not be changed by a factor less than 0.76 
.69 for the F, = 2 test and 0.78 to 0.67 for the Fj = 6 test. Actually the 
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apparent reduction in maximum velocity is nearer to 0.50 and 0.12 in the two 
cases. It is not clear why separation did not occur along the “bed” at an x/ 9 
location of about 1.2. 

The writer’s second reservation concerning the “air model” of the hy- 
draulic jump has to do with the upper or “free” surface of the jump and the 
roller. Since the streamline separating the roller from the main flow and tht 
bed pressure distributions are well verified, it appears that the model is ade 
quate. This mean flow similarity is certainly a necessary condition for sim 
larity, but it is essentially geometric in character. From a dynamic stand- 
point two other factors enter. The first of these is that in the air model a — 
finite wall shear must have occurred due to the air motion past the rigid upp 
boundary. In the actual jump no such shear occurs, at least if one neglects © 
air friction, The other divergence of the upper surface from the true hy- — 
draulic jump is less simply presented. In the truly-hydraulic jump the surfa 
and its roller appear to be far from steady. Instantaneous views of the uppeél 
surface of the jump tend to give the appearance of a breaking wave falling 
backwards on the oncoming flow. It is felt that the fluctuating nature of the 
free surface may have a significant effect on the flow at some distance below 
the surface even though none is apparent in the mean flow characteristics. 

The upper boundary, or free surface, differences between the air model 
and water prototype may be likened to the differences between the outer 
(distant from the wall) portions of boundary layer and pipe flows. Early 
boundary layer analyses assumed that the velocity profiles (and by impli- 
cation, the turbulence pattern) in the two cases were the same if the bounda: 
layer thickness was equated to the pipe radius. More recently careful stv li 
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ave indicated significant divergences, both in mean velocity and turbulence 
laracteristics, in the region away from the wall. G. B. Schubauer(2) has 
10wn that this difference in turbulence is completely explained through the 
1enomenon of intermittency which occurs only at the outer or free stream 
ige of turbulent flows. As sketched in Fig. B, the division between the turbu- 
tly flowing fluid and that outside the boundary layer, does not always occur 
; the same location, as in the case laminar boundary layers. Instead the 
terface, between turbulent and non-turbulent flows, fluctuates with time and 

3 conceived of as waves or wrinkles which travel along the outer edge of the 
er. 

The fluctuating nature of the outer edge of a turbulent boundary layer 
iffers vastly from the situation that obtains at the center of a pipe where the 
uid is always turbulent. A velocity indicating instrument with fast time re- 
ponse, when placed in the outer boundary-layer region will see patches of 
irbulent and non-turbulent fluid pass by. This condition is termed inter- 
littency and the fraction of time that the flow is turbulent is defined as the 
itermittency factor ¥. Schubauer(2) showed that if we consider only the time 
hen the flow is turbulent in the boundary layer the turbulence in the outer 

gion agrees well with pipe measurements. It is probable that the differ- 
ces in mean velocity profile are also tied in with intermittency. 

Intermittency cannot be immediately applied to our hydraulic jump. Even 

uniform open-channel flow, the free edge of our flow is subject to a strong 

itation and the air-water interface must inhibit if not eliminate the oc- 
rrence of this phenomenon, In the case of the hydraulic jump it would ap- 
ar that the unsteady roller may result in an intermittency-like behavior. 
rtainly the upper surface of the roller is rather intermittent in shape and 
sition. The occurrence of air bubbles in this region should also lead to 
termittency. As a further bit of suggestive evidence it may be noted that in 
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simple boundary-layer flows the wave or wrinkle length of the free boundary 
is of the order of two boundary-layer thicknesses. This is of the same order 
of magnitude as the relative length of the surface roller in the jump. 

Although the writer questions whether the authors have adequately 
“modeled” the hydraulic jump, he agrees that some very suggestive infor- 
mation has been obtained pertinent to the flow behavior in the jump. Actual — 
turbulence studies of a wet hydraulic jump will be needed to ascertain the ex= 
tent to which the model studied was appropriate. 

Aside from its hydraulic jump connotation, the study reported is an ex- — 
tremely valuable piece of basic research. Very few turbulence measurement 
have been made of separating and reattaching flows. And yet, these are quite 
common occurrences (in valves for instance). The only comparable study of 
which the writer is aware was also conducted at the Iowa Institute of Hydrauli 
Research by A. A. Kalinske(3) and co-workers. The sudden expansion and 
large angle diffusers studied in that project involved separation and subse- 
quent reattachment. The experimental techniques used were much cruder 
than in the present instances, but somewhat similar conclusions were reaché 
Thus Kalinske stated, “the maximum total turbulence energy is a small part 
of the total energy change taking place” and “the principal loss of energy oc 
curs at the outside of the high velocity stream as it passes into the expansion 
due to the very high shear developed.” 

The relatively rapid disappearance of the extra turbulence energy genera 
ed in their expanding flow is noted by the authors. Kalinske’s sudden ex- 
pansion results indicate almost as rapid a decay. Since to the observer of 
the hydraulic jump it still appears quite “turbulent” at its downstream end, — 
where Fig. 10 shows a low level of turbulence, one may question whether to — 
believe his eyes. Another possible explanation may account for this anomaly 
In other flow situations (decay of isotropic turbulence) the turbulence energy 
is dissipated through viscous action in the smaller sized eddies while a rela- 
tively small portion of it associated with the large scale eddies, is little ef- 
fected. In the case of the jump at the downstream station the turbulence level 
may be low but large eddies could still be quite prominent even though con- © 
taining little energy. This possibility could be verified through spectrum 
measurements of the turbulence. 

Study of Figs. 6 and 7 should give some idea as to the manner in which a 
separated and then reattached flow redevelops into a stable pattern again. 
Following the roller the turbulence levels are seen to be high, but the mean 
velocity profiles appear to ignore the upper surface. Analysis of mean- 
velocity data from other studies of reattaching flows have suggested to the 
writer that a new boundary layer forms along the wall more or less inde- 
pendently of what has gone on before. It may be that the measurements in th 
jump case did not proceed far enough downstream to indicate this. One of th 
flow cases analyzed, with the aid of H. R. Fraser, is depicted in Fig. C. Thi 
involved the flow of a jet into a partly throttled pipe.(4) Although the pitot — 
tube used was not calibrated for the reverse flows measured, a significant 
separation zone is apparent. The low values of the momentum thickness 6 
and displacement to momentum thickness ratio (6 */g) found at x/D of about 
8 suggest the start of a new boundary layer. Although not at all accurate the 
wall shear stress variation found in the separated region is of the proper sig 
and approximate magnitude. In contrast the apparent behavior of the wall : 
shear. stress along the upper curved surface in the authors’ investigation 
seems incomprehensible. In the separated zone the flow direction near the 
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surface is reversed and yet Fig. 7 indicates a negative wall shear stress 
within and following the roller. 

To calculate the turbulence dissipation the authors utilized the isotropic 
relation (constant K across and along the flow). They did adjust the value i 
a bulk fashion for the several Froude number models studied. Certainly K 
must vary throughout the flow field since the turbulence is far from isotrop: 
and the kind or degree of isotropy can hardly be uniform. P. S. Klebanoff( 
in his very thorough study of the simple boundary layer, measured a majoril 
of the terms appearing in the turbulence dissipation relation. In places som 
of these differed by a factor of ten from the isotropic values. A transverse 
integral of the dissipation, based on plausible assumptions for the unmeasur 
terms, indicated over twice the dissipation calculated with the isotropic re- 
lation. It is hard to believe that the more involved pressure gradient and 
separation flow case in the jump model would come out as close. Therefore 
the dissipation relations of Fig. 11 must be considered as of suggestive 
than absolute value. * 

Even within the limitations which it is felt should be placed on the results 
presented, the paper reports on an extremely valuable and informative studi 
The writer could only wish for more detailed information on the quantities 
measured in order that they could be studied independently of the hydraulic 
jump framework in which the paper is cast. 
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DONALD R. F. HARLEMAN,! A. M. ASCE.—The authors have clearly 
demonstrated yet another aspect of the universal nature of fluid mechanics | 
their detailed study of the hydraulic jump by means of an air model. 

As a supplement to an otherwise comprehensive treatment the writer wo 
like to set down the explicit equation for the exact depth ratio of the hydraul 
jump. This will permit an investigation of the small but systematic deviatit 


of the experimental depth ratio, dp/d,, from that predicted by the elementai 
momentum analysis given by, 


1, Associate Prof. of Hydrs., Hydrodynamics Lab., Massachusetts Inst. of 
Technology, Cambridge, Mass. 
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This deviation is shown graphically in Fig. 4 where the author’s hydraulic 
xperiments, together with those of Bakhmeteff and Matzke, are compared 
vith Eq. (a). At a Froude number F = 6, for example, the measured depth 
‘atio is five per cent lower than that given by the elementary theory. 

The analytical expressions for the hydraulic jump depth ratio taking ac- 
-ount of (1) boundary shear, (2) the excess of mean momentum flux due to non- 
iniform velocity distributions at the beginning and end of the jump and (3) the 
nomentum flux due to turbulence have also been presented by Van Driest.(1 
[he author’s paper, however, permits for the first time a quantitative evalu- 
ition of these terms from direct measurements. 

The exact equation for the jump depth ratio do/d, can be obtained by evalu- 
iting the indefinite integrals of Eq. (13) at the end of the jump where d = dg 
ind x = Lj; Ea. (13) then becomes, 


ae d, Pea da = Xd? xd Ly * (b) 
Se u.dy - fe u, dy + Se Wrdy = AS - Se fe dit) dx 


29 “y=0 


Two typographical errors in the author’s Eq. (13) have been corrected, i.e. 
he dy appearing in the last term should be dx and the sign of the last term 
should be minus, since the bottom shear acts on the fluid in the same direction 
is the downstream hydrostatic force.) 

Upon introducing the dimensionless coefficients 


Jeary " §.tezdy 
aay ri Biot sare 
nd 
(cikay 
T= TFd, 


here 8, and 9 are the momentum flux correction terms for non-uniform 
elocity distribution and Ig is the downstream turbulence flux correction term 
1 has been assumed equal to zero). After dividing through by 7 Eq. (b) may 
e re-written as follows: 


by 
R2Uld,/g - @.U7d,/g +T,Ud/q = tld di)-AFIG)ax © 


y introducing the continuity equation, Uy dj = U2 dy and defining the Froude 
nd Reynolds numbers as 


— U, pt Uid. 
lore ya Toa seetr 


q. (c) may be written as, 


re Sn [(T+)T- +5] | (a) 
ae Aa - (Ar.+12) 
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where d. 


and 


Eq. (d) should reduce to the form of the elementary momentum solution of 
the hydraulic jump if 8, = Bg = 1 and Ig = 0 and Ry = (frictionless flow, 
S = 0). Under these conditions, 


F*= V2 (T+!) 


which is in agreement with Eq. (a). 

The dimensionless quantities 81, B2, Ig and S can be evaluated from the 
author’s experimental data for each of the Froude numbers tested (Fy = 2, 4, 
6) ae the corresponding values of the depth ratio J can be calculated from 
Eq. (d). 

Table I gives the quantities determined from Figs. 7 and 8 of the author’s 
paper. Each of the downstream quantities 89, Ip and S were determined at 
the end of the jump x/dg = 5. 


Table I. 


Momentum, Turbulence and Shear Parameters 


F I Ss J 
1 PL Fe 4 eq. (a) oi 
2 1.01 1.02 OL 0.16 2eaL 2.37 
1.01 1.05 206 1.05 5.02 5.18 
6 1.01 1,02 «15 3.60 7270 8.00 


Also given in Table I are the values of the depth ratio J = do/d, from the © 
exact analysis, Eq. (d), and the elementary analysis, Eq. (e). The exact depth 
ratio is seen to be systematically less than the elementary. 

A final verification of the author’s measurements and the general validity 
of Eq. (d) may be demonstrated by comparing the three predicted (exact) 
values of the depth ratio with Fig. 4 of the author’s paper. This comparison 
is shown in Fig. A and indicates that the predicted values (shown by the dashe 
line) account for the departure of the experimental depth ratios from the 
elementary momentum equation for both the hydraulic and air model tests. 

The bottom shear S is by far the most important corrective term in the 
exact equation for the jump depth ratio. In fact, it can be shown that the mo- 
mentum flux and the turbulence flux corrections tend to cancel each other. — 

In this respect, the effect of bottom shear is similar to the effect of baffle 
piers in reducing the downstream depth of a hydraulic jump. (2) In the 
presence of baffle piers the distributed boundary shear parameter S may be 
replaced by a term descriptive of the more concentrated force exerted by the 
piers on the fluid. 
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Fig. A. Comparison of Exact and Elementary Momentum Equations 
for Jump Depth Ratio. 
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PHILIP G. HUBBARD, ! A. M. ASCE.—As stated by the authors, the analogy 
ween the air-flow system and an actual hydraulic jump may seem to be 
her artificial at first glance, so that some engineers may want further . 
tification or at least an evaluation of the applicability of the results. Be- 
ise of the thoroughness with which the study was conducted and the great 
portance of the conclusions drawn, it would seem desirable to take ad- 
itage of any technique which might supply additional evidence. The hot- 

nm anemometer mentioned in that paper is one such technique, and some 
oratory measurements with that instrument are presented herein. 


Research Engr., Iowa Institute of Hydr. Research, Iowa City, Iowa. 
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Fortunately, the flume used by Dr. Siao for the measurements in water if 
a permanent item of equipment at the Institute, so that it was possible to 
reproduce the conditions of those earlier experiments. In order to obtain 
representative results in a brief period of time, measurements were made 
each Froude number for which results were reported, but only for one repre 
sentative zone in the flow. The zone chosen was near the downstream end co 
the roller, where the longitudinal component of turbulence was measured ~ 
along a vertical line. Although a complete comparison cannot be made on th 
basis of this one component, it is probably the most representative single 
parameter. As mentioned by the authors, the total intensity of turbulence is 
very near equal to twice the intensity of the longitudinal component: 


=f A. if 
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When the measuring instrument was placed in the zone chosen, two perti- 

nent characteristics were immediately evident: The passage of innumerable 
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Fig. B. Record of Velocities Inside Roller, Showing Tranquil Intervals 


- bubbles past the probe had scarcely any discernible effect on the velocity 
lications, so that some earlier misgivings were found to be without basis; 
d the fluctuations were so erratic that only a rough indication could be ob- 
ned of either the mean velocity or the root-mean-square deviation of the 
stantaneous velocity from that mean. Although this problem had been en- 
untered to some extent in the air jump, it was greatly compounded in water. 
had been expected that the frequency spectrum of the turbulence would be 
ifted downward because of the lower velocities used in water, and this shift 
S expected to cause difficulty because the time constants of the turbulence- 
yzing instruments were marginally small. Even when the probable error 
to inadequate frequency range in the root-mean-square analyzer was 
<en into account, however, the measured intensities of turbulence shown in 
>. A were discouragingly below those shown in Fig. 7 of the subject paper. 
In order to examine more carefully the variation of the fluid velocity, the 
ut voltage of the hot-film control circuit was connected directly to the de- 
ction terminals of a cathode-ray oscillograph so that no error was intro- 
ed by coupling networks. Some typical records obtained in this way have 
n traced from photographs for reproduction in Figs. B and C, and it is be- 
ed that they contain the key to the discrepancy already noted. It was ob- 
ved that bursts of high-intensity turbulence were interspersed with periods 
relatively tranquil flow. Strangely enough, this intermittency existed not 
y near the upper edge of the boundary layer, but was quite pronounced 


-Fig. C. Record of Velocity Immediately Downstream from Roller 
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inside of the roller. If the data of Fig. A are modified by eliminating the 
periods of tranquil flow from the record, then the intensity of the remaining 
fluctuations is slightly higher than that indicated in Fig. 7. 

Any evaluation based upon the results reported herein can only be tentati 
pending a more complete investigation, but two conclusions seem warranted 
(1) The diffusion characteristics of the hydraulic jump are probably repre- 
sented quite well by the results of the subject paper, but those flow charac- 
teristics which depend upon the smaller-scale components are probably — 
slightly distorted by the higher indications of the intensity of fluctuation (2) 
Techniques and possibly definitions will have to be developed especially for 
flow situations where intermittency may be a factor. ~ 
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| THE TOTAL SEDIMENT LOAD OF STREAMS 


Discussions by R. J. Garde and M. L. Albertson, D. C. Bondurant, 
John L. Bogardi 


_ R. J. GARDE, 1 J, M. ASCE and M. L. ALBERTSON, 2 M. ASCE.—This 
aper isa welcome addition to scarce literature on sediment transport in 
lluvial channels. It will be of use to those who are working in the field, irre- 
pective of the approach in which they are interested—empirical or theoreti- 
al. 

_ Broadly speaking, the problem of evaluation of total sediment transport 
ate in alluvial channels can be handled in one or a combination of the follow- 
ng three ways: 


1. Based on the present-day knowledge of fluid mechanics and sediment 
transport, a theoretical relationship for the total load can be developed. 

2. By adding suspended load to the contact (or bed) load, the total sediment 
transport rate can be determined. Einstein’s procedure, and the Modi- 
fied Einstein’s procedure would fall under this category. 

3. Based primarily on dimensional analysis and intuition, some empirical 
relationships can be developed to express the total sediment transport 
rate as a function of parameters relating hydraulic characteristics of 
the flow and characteristics of bed material. 


Because of time and energy involved, and also because of the inadequate 
nowledge of sediment transport phenomena, very little work has been done 
1 developing the first method. 

Einstein’s procedure of finding the total sediment load by adding suspended 
9ad to the bed load is quite often used. It is found that the results obtained 
y Einstein’s procedure do not agree with certain data obtained in the field. (8) 
herefore, Hembree and Colby(8) have developed a “Modified Einstein’s Pro- 
edure” to evaluate the total sediment load. 

In evaluating the total suspended load by integration of the sediment distri- 
ution curves, some important assumptions are made. Following are given 
ie basic equations and the assumptions made in their development: 


Equation Assumptions 


d That the density P of the fluid-sediment 
en) dy ot suspension is independent of y and that 
Cém AV the kinematic viscosity V is relatively 
d insignificant compared with the eddy 
viscosity €m. 


ce 


it 


_ Proc. Paper 1530, February, 1958, by E. M. Laursen. 

_ Graduate student, Colorado State Univ., Fort Collins, Colo. 

Director of Research Foundation and Prof. of Civ. Eng., Colorado State 
Univ., Fort Collins, Colo. 
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Equation Assumptions 
OM Say Le ye That the logarithmic distribution of 
y <7 velocity is applicable and that the 


Karman universal constant k has a 
constant value. 


3g ee Aes op =o That the sediment transfer coefficient 
€, is equal to the momentum transfer 
4, Em= Es = Ky [TYP pay coefficient €,). 


At present, work is being done at Colorado State University on the evalu- 
ation of total sediment load by adding suspended load to contact load. Thus 
far, this method has been found to have serious limitations because in a dune 
regime one does not know how to predict accurately either the concentration 
of suspended load at a given elevation or the exponent z in the actual sedimen 
distribution curve. Furthermore, if one assumes a general law of velocity 
distribution as 


vy. 2 as 
Vad ames 1 Log 10 Ks (1) 


in which v is the velocity at a distance y above the bed 
k is the Karman universal constant 
Kg is some length parameter 
V* is the shear velocity 


it is found that both k and K, depend on the flow regime—such as plane bed or 
dune bed—and other flow parameters. | 

Therefore, for immediate use, at least some empirical relationships are 
necessary for predicting the total load. With this view in mind, the writers 
carried out certain preliminary investigations in 1956. The results of these 
investigations are summarized in the following paragraphs. 

While studying the sediment transport through pipes, (1) the writers found 
by dimensional‘analysis the following relationship to be applicable: 


ReVf =$(q > Ce) (2) 
in which D_ is the diameter of the pipe 


d_ is the diameter of sediment particles 


Cy is the concentration of total sediment load expressed in per cent 
by weight 


f is the Darcy-Weisbach resistance coefficient and 


Re is the Reynolds number = - 


For smooth pipes (as compared to standard:corrugated, or helically corru- 


gated pipes), for a given value of 2 it was possible to obtain plots as shown i 
Fig. 1. 
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FIG. t. VARIATION OF Reyf WITH SEDIMENT CONCENTRATION . 
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Although open-channel flow differs from pipe flow because of the existent 
of a free surface, it is reasonable to assume that a similar relationship 
should exist for the transport of sediment in alluvial channels. Therefore, 


V*D (which is proportional to Ref ) was plotted against Cy. A typical var 


V 
ation of ee with C; and 2 is shown in Fig. 2 using Gilbert’s data. 


d 
For each set of data, a definite relationship existed between Yad z Z 


Ct. Plots were made for flume data from various sources such as Kalinske 
and Hsia, (3) Brooks, (5) Barton and Lin, (6) and Gilbert.(4) The range of s 
of bed material included in the study was from 0.011 mm to 1.71 mm. 


transport rate C; by the following empirical relationship: 


(Nee a.) eke 

6) c,/3 a d a a 7 
in which ¥ is a function only of size of bed material. This functional rele 
ship is shown in Fig. 3. 

Fig. 4 shows Eq. (3) plotted for the various data already mentioned. 
(Laursen’s data for 0.04 mm and 0.10 mm are also included.) There is an 
average scatter of 40 to 50 per cent. Possible causes for this scatter are — 
that, in the development of Eq. (3) neither the influence of flow regime (suc 
as dune bed or plane bed) nor the influence of the Froude number on the tota 
sediment transport are taken into consideration. 

The investigations under progress at Colorado State University by the 
writers have shown that k, Kg, and z depend on the condition of the bed and, 
therefore, the rate of suspended load transport should depend on the flow re 
gime. Similarly, Tsubaki and others(7) have shown that under conditions of 
equal tractive force the contact load is smaller if ripples exist on the bed th 
if the bed is plane. Therefore, flow regime at least should be taken into 
sideration in the development of formulae for transport of sediment in alluv: 
channels. 
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D. C. BONDURANT,! A. M. ASCE.—There are several formulae available 
for the computation of sediment transport. Unfortunately, data available for 
the development of these formulae have seldom extended to the range of larg 
streams, and it has been found that their use in respect to such streams re- 
quires a liberal application of experienced judgment. The writer’s primary 
objective in reviewing this presentation was, thus, to determine its values 4 
limitations in respect to streams such as the Missouri River. It is concludec 
that the form of the author’s equation can be very useful but that some re- 
vision of the curve of 1(7“¢/p /) ) is indicated. 

Data for the purpose of the review were available from measurements 
made in the Missouri River in connection with the sediment investigations of 
the Missouri River Division, Corps of Engineers. Of the more than 400 sets 
of measurements of vertical distribution of velocity and suspended sediment, 
76 were selected to provide the maximum range of data. Within these data, ~ 
stream depths ranged from 7 feet to 22 feet, velocities from 3 fps to 9 fps, 
and discharges from 25,000 cfs to 90,000 cfs. Stream widths approximated 
1,000 feet and slopes were from 0.8 to 0.9 foot per mile. From 6 to 8 verti- 
cals were occupied at each range, and the ranges occupied were in 4 separate 
reaches located at (1) River Mile 640, above Omaha, Nebraska, (2) Mile 614.1 
2-1/2 miles above the mouth of the Platte River, (3) Mile 608.2, 3-1/2 miles 
below the mouth of the Platte River and (4) Mile 365.6, below Kansas City, 
Mo. and 14 miles below the mouth of the Kansas River. Reaches were as 
nearly straight as could be found, and were unobstructed. Measurements wer 
made from a boat held in place on a tag-line. 

Figs. 15, 16, and 17, respectively show typical plots of velocity, suspended 
sediment, ‘ed Peds material. Velocities were measured with a Price-type 
current meter; suspended sediment samples were obtained with a U. S. P- -46 
(point-integrating) sampler; and bed-materials samples were obtained with — 
either a B.M.-48 or B.M-54 sampler or both. The B.M. instruments sample 
only the top one-inch (approximately) of the bed, and are so constructed that 
the sample is effectively protected against washing or contamination. . 

Fig. 18 shows a comparative plot of the average suspended sediment 
concentrations at the three upstream reaches as determined from the sample: 
and as computed by the author’s method. The measured concentrations are 
four times as large as the computed; however, the relationship is reasonably 
consistent. Furthermore, the measured and computed size distributions we re 
in excellent agreement. Additionally, four points provided by Mr. J. C. Pyle 
of the U. S. Army Engineer District, Little Rock, Arkansas, showed almost 
exactly this same relationship of four. These latter points represented dis- 
charges of 25,000 cfs, 73,000 cfs, 135,000 cfs, and 190,000 cfs in the Arkansat 
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er. The writer has no information as to the magnitude of other factors 
ylved. 

Jn the basis of the consistency indicated in Fig. 18, it appeared that 

her computations might be unnecessary. However, the bed-materials in 
Kansas City reach, as may be noted in Fig. 17, are substantially larger 

1 those of the upstream reaches, thus they provided an opportunity to ex- 
i the range of the data used. With these data the consistency of the re- 
onship deteriorated rapidly, the ratio varying from 10 to 100 and indicating 
need for a more basic review. 

inspection of the computations revealed that the values of ep /, associ- 
i with the larger bed-materials of the Kansas City reach caused the corre- 
nding function values for the major size fractions to be much smaller than 
eared reasonable. Accordingly, re-computations were made to derive 

1es of the function corresponding to the measured data. 

Figs. 19 and 20 show plots of f(U,/w) computed for the upstream reaches 
Kansas City reach respectively. Since all the writer’s plotted data use 
symbol “U,” for the more formal Vt/p , that symbol will be used for this 
ussion. It may be seen that the points plot with a reasonably limited 

ter except for those points representing the two smaller size fractions 

a few points representing the largest sizes in the respective samples. 

e deviations point to a systematic error in the data in respect to these 
tions, and investigations made subsequent to the analyses of the samples 
ated that such error did exist. 

he samples were analyzed by sieves, for reasons outside the scope of the 
ent discussion. Fall velocities at various temperatures were later de- 
ined by an extensive investigation involving both individual grain drops 
nalyses of individual sieve fractions by the visual accumulation tube 

od. It was noted that the smaller size fractions periodically contained 
tities of heavy magnetic minerals and, if the fraction represented only a 
percentage of the total sample, the average specific gravity of the 

on might be as much as 2.85 and in a few cases greater than 3. Had 
points been corrected for the higher fall velocity which would have re- 
d from the presence of these minerals, they would have been moved to 
eft and would have been essentially in iine with the other points. 

e less consistent variations for fractions of the larger sizes were found 
due to low density agglomerated particles. They were present only in 
ient quantity to affect a fraction which was less than about one per cent 
total sample. 

might be well to note here also that the percentages pertaining to the 
fractions of the bed-material samples were read from the plotted curves 
r than from the original laboratory records. Inspection of the compu- 

s showed that the possible error thus introduced was negligible insofar 
concentration computations based on Missouri River data were con- 

d; i.e., those fractions for which the error of reading might be appreci- 
ere of themselves such a small portion of the whole that a reading 

of even several hundred per cent would make little difference in the 

all result. Time did not permit a return to the original analyses for the 
on computations. 

g. 21 shows individually the lines drawn as the approximate average of 
otted points for the upstream reaches and the Kansas City reach. The 

d between these lines is not excessive, however, they are roughly paral- 
hdicating a systematic difference. A third line on this figure indicates an 
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‘rage for all points representing fractions which were more than 5 per cent 
the total respective samples, but this line is not presented as a recom- 
nded mean. The difference between the lines for the upper and lower 
ches is believed to result primarily from the limits used to determine the 
‘rage concentrations from the measured data. 

If one assumes that the logarithmic distribution of concentrations extends 
y = 0, the concentration at that point becomes infinite; thus, one must as- 
ne a limit at which y equals some real value. For purposes of prior 
nputations, this limit had been adopted as y = 2d with an arbitrary maxi- 
m of 480 g/1 which was not to be exceeded at that level. The average 
icentrations thus determined were used for the present computations. In- 
ar as the three upper reaches were concerned, these averages were 
sonably consistent with averages shown by depth-integrated samples 

ich could not extend closer than 0.5 foot to the bed). In the Kansas City 
ch, where the material was larger and thus necessarily concentrated more 
Wily near the bed, the difference was substantial. 

Experiments of Einstein and Chien26 indicate that the logarithmic distri- 
ion of suspended sediment begins to break at the level where the total 
centration approximates 100 g/1. Inspection of the computations for aver- 
concentrations at Kansas City indicates that assumption of any limiting 

e below 480 g/1 would reduce most of the computed average concen- 
ons, and assumption of a value below 350 g/1 would reduce all of them. 

iS would move the computed values of £(U,/w) for the Kansas City data to 
right, more nearly approximating the position for the upstream data. 

ime has not permitted recomputation of the data to verify, quantitatively, 
variations above discussed. Even with these variations, however, the 
souri River data provide a consistent plot within their range, without ex- 
ive scatter, and verify the form of the author’s equation. A plot of the 
or’s data (roughly scaled from Fig. 7), included in Fig. 21, indicates that 
verage lines from the Missouri River computations fit the author’s data 
flower values of U,/w as well as or better than the curve necessarily con- 
tcted by the author to fit his data for the higher values of U,/w. For the 
er values of U. ,/W the Missouri River results deviate sharply from those 
e author: 

he Missouri River data are sharply restricted in their range in respect 
th sediment size and slope of the energy gradient. The writer does not, 
, suggest that a function derived from these data be given arbitrary pre- 
nce over a function derived from other data. It is to be hoped that other 
ers will further augment the available information to permit the author 
‘ovide an integrated result. 

considering the fit of the Missouri River data, it is well to remember 
all pertinent factors in these data represent measured values in which 
orrect values of other functions, such as roughness or turbulence, are 
icit. It has been found that equally valid computations can be made with 
eoretical type of equations (for example) if the measured values of k (or 
re utilized. The author’s function has the advantage that velocities can 
ore simply measured than k values; thus, in cases where the velocity 

e adequately established, this function should be the more valuable. In 


. A. Einstein and Ning Chien, Effects of Heavy Sediment Concentration 
ear the Bed on the Velocity and Sediment Distribution. University of 
alifornia Inst. of Engr. Research Series No. 33, Aug. 1955. 
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cases where velocity must be computed from assumed roughness factors, the 
necessary dependence upon judgment and the possible variations in results © 
would be about the same as for other available formulae. 


JOHN L. BOGARDI.1—All credit is due the author for his successful in-— 
vestigations into the subject of sediment transportation. In keeping with the 
latest concepts, bed-load and suspended load are not separated inthe — 
treatment, and the possibility for relationships of a more general nature is © 
opened thereby. Numerous aspects of this complex problem are elucidated ~ 
by clear and logical reasoning. The writer fully agrees with every statem 21 
made by the author and only wishes to contribute a few additional remarks, — 

In the writer’s view, the relationship shown in Fig. 7 applies to the 
sediment-transporting capacity of the stream rather than to the actual rate¢ 
sediment transport. Quite obviously, the actual rate of transport will depend 
on a number of factors besides the transporting capacity, e.g. also on the 
availability of sediment in the stream. It follows that the actual rate of 
transport in natural streams is governed also by the hydrological conditions: 
of both the watercourse and the catchment area drained thereby. The writer 
interprets the excellent agreement between predicted and actually observed 
rates of transport on the Niobrara River as an indication of the fact that 
hydrological conditions of the watercourse, as well as sediment supply, are 
such as to enable the stream to carry sediment in an amount equal to its 
transporting capacity. A rate of transport in excess of the transporting ca- 
pacity and observed on Mountain Creek is, as mentioned also by the author, 
due probably to artificial flash floods. Consequently, the author is fully jusi 
fied in arriving at the conclusion that the prediction of the so-called wash lo; 
is not feasible. : 

Relationships (8) and (9) published by the author gave the writer the oppol 
tunity to extend his investigations, conducted on the basis of results publishe 
by Liu(1) and by Albertson, Simons, and Richardson, (2) to total sediment 
transport. In his paper Liu established a definite relationship between U,/W 
and U,d/v for the state of incipient ripple formation. Albertson and his 
collaborators have further shown the existence of a similar relationship be- 
tween the two above dimensionless invariants also for the beginning of 
movement, the development of dunes, transition, and antidunes, respectivel} 
The writer investigated the significance of the channel stability factor 


was named, following the example set by Albertson and in analogy to the shi 
velocity Reynolds number, the shear-velocity Froude number. Using the © 
channel stability factor b = 1/ Fr, the writer succeeded in defining numeric: 
various regimes of movement by a single parameter 8. Investigations of th 
writer, conducted using data and results published by Albertson and his col 
laborators, showed a relationship between the parameter B and the shear- 
velocity Froude number that can be expressed as 
OE a es 
hs U Fry, rc 
1. Asst. Prof., Technical University, Budapest. Head, Experimental Secti 
of the Scientific Research Inst. of Hydroeconomy. 
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ce for all regimes N was 0.882 = constant, the parameter 


oo n8 
U,-6 “OER * 


ines uniquely the bed Pe bans The Beast, between the particle 
ed and 1/ Fr, is shown in Fig. 15. It is to be seen from the figure that the 
“ameter 8 (of the dimension cm-0.882) in effect defines uniquely the bed 
‘mation. 

A number of important conclusions can be arrived at from Fig. 15. 

ples, for instance, cannot develop unless the particle size is smaller than 
roximately 2 mm, and even for these particle sizes only if 66< 8 < 322. 
Fig. 15 relates to a water temperature of 20° Centigrade and to a sediment 
fring a specific gravity vy = 2.65 gr/cc. For other temperatures, e.g. for 
emperature of t° Centigrade, 


B. = 8 Ye, t Vo0 
20 2, 30-088 ve 


rein 99 and v; denote the kinematic viscosity of water at 20° and t® 
tigrade, respectively. 
or sediments having any afar: 4 specific gravity V1 


B. rv)! i Qev) 
Sieey = 


65y Fr a” 1.657 
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The interrelation between U, and d can also be determined from that es 
tablished for B and Fr, Thus 


1-1 

V2) 
a s) 2 
-(g\"" a 


Since, however, a = 0.0582 = constant for every bed configuration, ‘ 


(i) 
of the dimension cm9-9418 sec7l, results in the relationship 


ies eo 0562 


*% , 
This relationship, derived also on the basis of data published by Albertson, | 
shown in Fig. 16. As can be seen from the figure, the configuration of the 
bed is uniquely defined for any arbitrary particle size by the parameter 


introducing the parameter 


= | 0500 
ds a @ 


It should be noted that all sorts of bed configurations are uniquely define 
similarly to parameters # and € by the parameter 


so ee: Se 
re n ( rata 
Fa 


Uyd 
Ry 


The parameter € is also related to a water temperature of 20° Centigrac 
and to a sediment having a specific gravity Y; = 2.65 gr/cc. For water — 
temperatures and sediments having a specific gravity other than mentioned 
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ove, the parameter can be determined from the following two equations: 


7 ~y -1/2 
ey "8565 1.65 


tails pertinent to the above investigation are to be found in the writer’s 
scussion of Liu’s paper. 

In the writer’s opinion, if Fig. 7 presented by the author actually expresses 
= relationship between the sediment-load parameter corresponding to the 
ansporting capacity of the stream and U ,/W, the variation in the sediment- 
ad parameter 
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can also be entered into Figs. 15 and 16. 
Values of U,/w can be established for all five bed configurations in the — 
case of several particle sizes from Fig. “A” published by Albertson and his 
collaborators; (2 moreover, values of U,d/v can also be determined. The 
latter lends itself readily to the computation of values of U, and subsequent 
of those of the parameter gd/U,2. For each value thus determined, those 0: 


Cc 
7/6 ‘ 
a fo iY 
YO ‘ce 
pertaining to various U, /w values can be correlated using Fig. 7 published — 
by the author. 

The relationship between 


and 


can thus be represented in Fig. 17 for particle sizes in Albertson’s Fig. “A”, 
Thereafter, relying upon Fig. 17, the isometric chart of the sediment load 
parameter ' 


can be traced in Fig. 15, Finally using Fig. 15, the isometric chart of the 
sediment load parameter can be entered into Fig. 16 as well. : 
It should be remembered that Figs. 15, 16 and 17 were prepared on the — 
basis of values obtained from Fig. “A” and Fig. 7 published by Albertson am 
the author, respectively, as no numerical values were available to the write: 
Some discrepancy may easily be introduced thereby. 
It should be noted further that the writer completed Albertson’s Fig. “A” 
by extrapolating Shield’s curve towards smaller particle sizes, and, on the 
basis of Fig. 15, slightly modified at the diameter d = 0.0011 cm the curve 
characterizing the five bed configurations by transferring observation po nts 
of Kalinske and Hsia to higher U, /w values. ; 
The writer considers Figs. 15 and 16 very useful for design purposes, 
since for any given watercourse, having established the hydraulic radius 
(water depth), the slope, and the particle size of the bed-material, paramete 
B and € may be computed directly, and the value of one of the two is alone ~ 
sufficient to define uniquely the bed configuration, and also the value of the 
sediment-load parameter, which in turn permits the computation of the tota 
sediment-transporting capacity of the stream. q 
It should be noted further that no accurate limiting values have been de- 
termined towards large and small particle sizes in deriving the relationshy 
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d in constructing the figures. These are felt to require further study and 
periments. For example, the line representing the formation of antidunes 
Fig. 15 is likely to join that representing the limiting condition of a smooth 
d through a curved transition instead of a definite point of intersection. 
milarly, the determination of relationships for values smaller than 7 and 
zher than about 1000 to 3000 of the sediment-transport parameter is likely 
require further investigation. 

The number of particular and general conclusions that can be arrived at 
means of the relationships described above and by Figs. 15, 16, and 17 con- 
ructed therefrom is great. It is, however, the intention of the writer to 
vote a special paper to these problems. 
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EISENHOWER AND GRASS RIVER LOCK MODELS@ 


Discussion by Marvin J. Webster 


MARVIN J. WEBSTER,! M. ASCE.--The authors state that the degree of 
rging and consequently the magnitude of force imposed upon a craft in the 
k chamber reaches a minimum value when the manifold system of ports or 
erals occupies the approximate middle third of the lock chamber walls. 

sts performed by the Bonneville Hydraulic Laboratory, Bonneville, Oregon, 
a 1:25-scale model of The Dalles navigation lock, Columbia River, and in 
prototype furnish additional data to substantiate this statement. 

The Dalles lock has a culvert system in which seven lateral culverts from 
*h wall culvert are located in the middle third of the lock chamber (see 

; 9). Model and prototype tests showed very low hawser forces on a com- 
rcial tow of 2780 tons(1 

The tests in the prototype were made in December 1957.2 The tow consist- 
of two barges and a tug having a total displacement of 2780 tons. Instru- 
ntation for the tests was developed and constructed at the Bonneville Hy- 
ulic Laboratory. The bridles for holding the tow in position were con- 
ucted of 5/8-in. cable and connected to the floating mooring bits, as shown 
Fig. 10. The tow was centered over the lateral culverts. The longitudinal 
ces were measured with two 20,000-lb-capacity dynamometers, one at 

+h end of the tow. No attempt was made to measure transverse forces in 
prototype because of the difficulty of constructing a suitable bridle to 
asure transverse forces and because previous tests in the model on larger 
e tows showed that transverse forces were less than the longitudinal forces. 
stS were made in the 1:25-scale lock model prior to the prototype tests 

., although the displacement was simulated, the shape of the model tow 
ered from that of the prototype. 

Results of the tests during filling for a lift of 84.6 ft in the prototype and 

0 ft simulated in the model are shown in Fig. 11. The maximum lift of 

> Dalles lock is 90.5 ft. The test results show low hawser forces during 
ing operations on the tow used in the tests and good correlation between 
model and prototype for filling. An emptying test in the prototype with the 
ae tow showed a maximum longitudinal force of 1.1 ton. An emptying test 

; not made in the model. Hawser forces also were measured on the same 
in the McNary navigation lock, which has a culvert system similar to The 
es, for a lift of 87.9 ft (maximum lift is 92.0 ft). The maximum 


Proc. Paper 1582, April, 1958, by Martin E. Nelson and Harvey J. Johnson. 
Chf., Hydr. Design and Lab. Section, U. S. Army Engr. District, Portland, 
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On a tow furnished through the courtesy of the Tidewater-Shaver Barge 


Lines, Portland, Oregon. 
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longitudinal force measured during filling was 1.7 tons. The filling valves 
open in 7.6 min and the lock filling time is 15.0 min at McNary compared ¥ 
4 min and 12.5 min at The Dalles, which probably accounts for the lower 
hawser forces at McNary. 

The writer wishes to acknowledge the use of material for this discussi0 
from the U. S. Army Engineer District, Portland, Oregon. Instrumentation 
for the tests was developed and constructed under the supervision of C. 
Robert Grim, Head of the Electronics Unit under the general supervision 0 
Harry P. Theus, Chief, Bonneville Hydraulic Laboratory and the writer. 
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WATER DISTRIBUTION DESIGN AND THE 
McILROY NETWORK ANALYZER2 


Discussions by Quintin B. Graves and Don Branscome, E. Shaw Cole, 
Claud C. Lomax 


QUINTIN B. GRAVES,! M. ASCE and DON BRANSCOME,?2 J. M. ASCE.— 

e authors have presented an interesting approach to the design of water 
tribution systems with particular emphasis on the advantages of the 

Ilroy Network Analyzer in providing a visualization of the hydraulics of a 
er distribution system. 

This paper and a paper by the writers on “Digital Computers for Pipeline 
work Analysis” were presented at the October, 1957 ASCE Convention, 
before different divisions and were published simultaneously in the Pro- 

ngs. 

t is obvious, as the authors point out, that by visual inspection of the 

tem as set up on a McIlroy Analyzer design changes will be quickly indicat- 

nd may be easily adjusted at “Inflow rates well below the design rates”. 

S is an advantage for the McIlroy Analyzer over the electronic digital 

puter primarily in those cases where changes in the distribution of in- 

S or outflows are indicated. In the case of the electronic digital comput- 

, computations on the basis of the largest rates of demand usually will 

ve to be advantageous. This may not be required in the case of the McIlroy 
yzer. 

io doubt, improvements in the techniques of utilizing the McIlroy Analyzer 

e been made since its first introduction, and more probably will be made 

continued use. Also, improvements in programming for digital computer 

ysis of distribution systems will be made. It has been evident from “long 
” calculations with the Hardy Cross Method that convergence of the 

ess does not always occur. Usually, by changes in the initial flow as- 

jptions for the network it is possible to obtain rapid convergence. The 

ers are in agreement with the authors in regard to the use of “electric 

ice power-system calculating boards with linear resistors”. 

me of the advantages of electronic digital computers over other types of 

yzers is their wide distribution throughout the United States. One of the 

ts to occasional users of McIlroy Analyzers is that of travel and associated 

ense for those who do not have one available for use locally. For most 

s of digital computers, such costs will be substantially reduced due in 

to the greater number of such computers available. As greater use is 
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made of digital computers for water distribution analysis, many improve 
in programming will be made thus enhancing their value. 

The authors mention that 50 to 60 punched cards will be required for b 
instructions, regardless of the size of the network being studied. It is po i- 
ble to reduce the number of basic instructions required by using an inter- 
pretive programming system. For most networks only 30 to 40 punched car 
will be required for the basic instructions. 

The authors, in discussing Fig. 6, state “A revised storage rate requires 
a new storage curve”. The writers are of the opinion that this would be mor 
intelligible if the statement read, “A revised uniform pumping rate requires 
a new storage curve.” 

Values of the Hazen-Williams coefficient “C” are not accurately known i 
many small water distribution systems. In such cases the variation in 
voltages or currents, which has been mentioned to the writers by some 
of the McIlroy Analyzer are of little or no concern. In those few cases wher 
basis information is more accurate, a variation of 5 per cent or greater maj 
be significant. 

The authors’ statement that “(However, a 2.00 power rather than 1.85 
power is involved)” in reference to the electronic digital computer program: 
is no longer valid as programs are available utilizing the 1.85 power. 

The authors hold that “an error in programming or data instructions is ff 
evident until the run is made and the analysis is completely summated and — 
checked, ...” In general after a program is once checked out all errors ~ 
should be eliminated and would not introduce errors when applied to other © 
water distribution systems. Errors in data may occur for any given case, — 
however automatic machine checking of data such as the summation of infle 
and outflows at a junction may be included in a program prior to the analysis 
of the system. 

The authors are to be highly commended for this fine contribution to the 
literature on water distribution system design. : 


E. SHAW COLE,! M. ASCE.—The authors have made an important contr 
bution in water distribution design and have shown how complicated problem 
can be solved readily by use of the McIlroy Analyzer. The writer has been 
privileged to work with them on some of their problems and has used the 
analyzer on a number of design problems for other systems. He believes a 
authors have presented thoroughly the advantages of the analyzer, as well as 
the limitations in its use, but would like to mention several points that may 
not be apparent to the casual reader. 

The importance of accurate field data is mentioned. This should be obv: 
and was recognized by Dr. McIlroy at the outset of his work in developing th 
analyzer. Unless the C-Value is reasonably accurate, the answers given by 
the analyzer may be seriously in error, but accurate field data implies mol 
than determination of C-Values. 

Frequently, water distribution systems do not function as intended, due t 
closed valves or obstructions that reduce carrying capacity far below the 
most conservative estimates. Some water works designers may not apprec 
ate the difficulties of maintaining distribution systems at peak efficiency, DI 


1. Pres., The Pitometer Associates, New York, N. Y. 


ICE DISCUSSION 1856-87 


e discovery of major difficulties is becoming more frequent every year. 
le increasing complexity of our systems and the lack of responsibility, or 
mplete records of the field operations, are major causes for this condition. 
The problem of determining the quantity of equalizing storage necessary 
permit uniform pumping is most interesting and illustrates the advantage 
the analyzer. In many water works systems, uniform pumping is not 
actical and pump units are added or subtracted periodically as the load 
anges. The pump changes may only occur a few times during a twenty-four 
ur period but this can be sufficient to make effective use of the existing 
orage. For medium size systems, the cost of pumps is relatively much less 
an storage, particularly if elevated, and it would be interesting to have a 


udy made of the economics of variable pumping under various conditions 
. uniform pumping. 


CLAUD C. LOMAX,! A. M. ASCE.—This paper is a welcome addition to the 
erature on the analysis of fluid distribution networks. Use of the McIlroy 
alyzer to study the effect of demand ratio and pumping station outages, and 
balance pumpage and storage demonstrates its versatility. As engineers 
come more familiar with the potential of this instrument, increased use for 
er than the routine head loss and flow rate determinations can be anticipat- 
. The R. L. Albrook Hydraulic Laboratory of the Division of Industrial 
search, Washington State College, has had such an analyzer since 1951. 

sS analyzer has been used for high and low pressure gas and for water 
tribution networks. The usual analysis has been determination of head 
ses and/or flow rates for present and future demands. Unusual problems 
lude analyses for Anchorage and Fairbanks, Alaska to insure all pipes had 

minimum flow required to prevent freezing. The distribution of cooling 
ter for Priest Rapids turbines was also analyzed. The City of Spokane in- 
ded the pump characteristics as influenced by the system pressure. Bend, 
egon used this analyzer to assist in selecting a location and water surface - 
vation for a reservoir being added to their system. The City of Los 

eles used the analyzer to study drawdown characteristics of their reser- 
rs. Pressure drops and distribution of bleed air for B-52 jets were 
cked. 

Other possible uses include economic studies related to pump selection, 
ervoir size and location, and main size. Stagnation potential of reservoirs 
be evaluated. Leaks and blocking can be pinpointed when a few field 
ssures have been taken. Hazen-Williams C values can be determined or 
irmed from field data for a badly tuberculated system. In time of dis- 
er this analyzer could be invaluable in re-routing gas and/or water for the 
st efficient use of the remaining system. A booster pump which will raise 
ssure at a selected point in the network has been added to the R. L. Albrook 
ratory’s installation. 
he authors’ comparison of the digital computer and McIlroy analyzer con- 
ns several of the conclusions of this Laboratory. The great value of the 
yzer over a digital computer in presenting an immediate and visual analy- 
of a system should be emphasized. While such a picture is very useful to 
engineer, it is also valuable in helping administrative personnel, fire 
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chiefs and water superintendents understand the distribution systems, The 
analyzer is probably a better design tool than the digital computer since the 
effect of changes are easily and quickly evaluated. 

There are disadvantages to the McIlroy analyzer which are not easily — 
solved. The number of pipe sections which can be included are limited by tf 
available Fluistor sockets. The identical pipe sections which can be placed 
on the analyzer are limited by the Fluistor stock, and it is not practical to — 
maintain a large stock of identical Fluistors for several pipe sizes. Solutio 
of the Darcy equation takes more time, and a trial and error fitting of the 
Fluistors is necessary. 

A program for the IBM 650 digital computer to supplement the McIlroy ~ 
analyzer has been developed by the R. L. Albrook Hydraulic Laboratory an 
the Computing Center at the State College of Washington. It is believed that 
the computer will be better for system analyses when there are few changes 
anticipated and when (1) the Darcy equation is used, (2) the piping system ec 
sists of many identical sections, (3) the number of pipes greatly exceeds th 
analyzer capacity and/or (4) great accuracy is required. The combination 
the analyzer and the computer will provide complementary rather than com 
petitive services. : 

The possibilities of the digital computer and the analyzer have not been 
fully developed. As increased use of each continues it is hoped future com- 
parisons will be made. There seems to be an established need for each. 
Factors influencing the selection of one or the other are not clearly es- 
tablished at this time. 
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NORTHEASTERN FLOODS OF 1955: METEOROLOGY 
OF THE FLOODS2 


Discussion by H. Alden Foster 


H. ALDEN FOSTER, ! M. ASCE.—Hydraulic engineers and hydrologists are 
ticularly interested in the effects of hurricanes or other large cyclonic 
urbances on the production of floods; hence it may be in order to point out 
tain lessons that may be learned from the authors’ analysis of the storms 
955. 

ne interesting fact is the extreme magnitude and expanse of the forces 
lved in such cyclonic storms. Thus, the total volume of water carried 
the continent in one 24-hour day of hurricane Connie is given as 562,000 
re-mile inches, or about nine cubic miles,—a lot of water even if spread 
a large land area. Also, during a single week in the life of such a hurri- 
, a large portion of the atmosphere over the tropical Atlantic Ocean could 
rawn into the storm pattern. Evidently, hurricanes mean not only ex- 

ive winds but also immense volumes of water that may be released to the 
ams and rivers in a short period of time. 

here is an extreme variation in location of the tracks of these tropical 
ms. They originate over the open ocean far south of the northeastern 

s, but during their travels may end up almost anywhere in this area. 

r intensity is maintained largely because atmospheric conditions near the 
Atlantic coast during late summer or early fall may at times be simi- 
o those commonly occurring in tropical regions where these storms 
inate. 

he travel paths of hurricanes are quite unpredictable for more than a day 
0 in advance. These paths are controlled by the general tendency of all 
Dspheric masses in motion to curve to the right in the northern hemi- 

fre, but they are also strongly affected by atmospheric conditions in their 
ty, such as occurrence of high or low pressure centers or troughs. 

he storm may pass over the open ocean throughout its journey to the 
fhern latitudes, in which case it is most likely to cross the coastline of 
Island or southern New England. Or it may pass inland south of Cape 
eras, and wear itself out while passing northward more or less parallel 
coastline. Where the storm track crosses the coastline, exceptionally 
tides may be expected, not necessarily at the location of the path of the 
ror “eye” of the storm but more likely some distance ahead of this lo- 
n,—that is, further north or east. 

e precipitation rates caused by a hurricane are not necessarily a function 
wind velocity. Even after the storm has grown old and is about to 
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expire, it may still carry an immense load of moisture in the circulating — 
atmosphere and is quite capable of dumping this water in places where it c 
cause severe floods. 

Theoretically, the major portion of the precipitation should take place 
ahead of and somewhat to the right of the eye of the storm, although topo- 
graphic or orographic influences may tend to affect this tendency; or the air 
masses carrying excessive quantities of moisture may be deflected from th 
natural paths by atmospheric conditions in the general vicinity. As a result 
the track of the storm does not necessarily indicate the location of the maxi- 
mum total precipitation during the storm. For instance, the total precipitatior 
for Connie was generally to the right of the track, while for Diane most of the 
precipitation in southern New York and New England appears to the left of the 
track. This seems to be due to the fact that, on August 16 when Diane was 
moving northward over Virginia, the precipitation was mostly located far to — 
the northeast, from Pennsylvania all across southern New England. By the © 
next day, Diane had made an abrupt turn to the right and passed south of ar 


area where the rain had mostly fallen the day before. If Diane had not change 
her course so abruptly but had followed a track eastward near the southern 
boundary of Massachusetts, the precipitation in southern Connecticut might 
have been considerably increased by bringing in some of the water which actu 
ally must have fallen over the ocean south of Long Island in the right-front 
quadrant of the storm as it moved along. 

This discussion brings up the important question as to whether the storms 
of 1955 produced the maximum precipitation that may ever be expected in the 
affected areas. Studies made in 1955 by Dr. A. A. Koch under the writer’s — 
supervision show that no storm of record in these Northeastern states has 
yet closely approached the maximum possible value for precipitation. Itis — 
quite possible that greater precipitation and more severe floods may occur in 
these areas in the future, though the probability of occurrence of such floods” 
may be very small. On the other hand, the maximum possible storm tides 
have been approached on several occasions at certain locations on the coast 
of southern New England, although other Atlantic coastal regions have not yet 
experienced the greatest possible storm tides. 

It is interesting to note in a news report of August 23, 1958, that the 
Weather Bureau has started a comprehensive study of storm tides along the 
3,000 miles of the Atlantic coast, from Texas to Maine. It is to be hoped the 
such a study will show whether the probability of occurrence of hurricanes i 
this region may be greater in the future than indicated by past records. 
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NORTHEASTERN FLOODS OF 1955: FLOOD HYDROLOGY2 


Discussions by Gordon R, Williams and H. Alden Foster 


GORDON R. WILLIAMS, ! M. ASCE.—Ample evidence is given in this paper 
hat floods resulting from extraordinarily large concentrations of precipi- 
ation are hydrologically and hydraulically a different phenomenon from floods 
laving recurrence intervals of 10 years or less. No longer should engineers 
und mathematicians compile frequency relations on the assumption that all 
loods comprise a homogeneous statistical population. The writer has previ- 

usly emphasized this point. (1) For river basins with areas of less than about 
,000 square miles in rolling or mountainous country, the major floods often 
efine frequency curves which have steeper slopes than the curves for ordi- 

ry floods. These relations are particularly pronounced in many 
alifornia(2,3) streams and elsewhere.(4) Rowe(3) has used the term bilineal 
requency curve. Curves with this characteristic are also common in New 
ngland streams.(5) Unfortunately the number of points available to define 

€ upper curves are too few so the problem of determining reliable recurrent 
ntervals for rare floods is not solved. 

The most promising approach to frequency determinations is the regional 

r composite frequency curve method that has been adopted by the U. S. 
eological Survey(6,7) and others. (3) Briefly this method uses the mean 
ual flood or the 5-year flood as an index to evaluate the magnitudes of 
oods of equal frequency. Within similar physiographic regions, floods of 
ual frequency have quite similar magnitudes when expressed as ratios to an 
dex flood. Thus it is possible to use the sum total of flood experience in a 
egion to develop a composite dimensionless frequency curve and to obtain 
ore reliable estimates of the frequencies of major floods than’ can be ob- 
ined by a single record. Dalrymple, in another paper in this symposium, 
s used the ratio to the mean annual flood to indicate the relative magnitudes 
f the 1955 floods. 
Correlations between index floods and size of drainage area are improved 
y introducing factors such as the lagtime for the unit-hydrograph 8) ora 
ctor which is a measure of average basin slope.(7) 
At least two conclusions can be drawn from this paper. One is that the 
onservative meteorologic, hydrologic, and hydraulic design criteria(9,10) 
t have been used in the past by the Federal engineering organizations are 
Ow justified. Modifications to the criteria will be inevitable as operating 
perience is accumulated but the basic conservatism must be retained. The 
econd conclusion is that experience and judgment in hydraulic engineering 
annot be replaced by stereotyped design procedures. 


Proc. Paper 1663, June, 1958, by Elliot F. Childs. 
. Prof. of Hydr. Eng., Massachusetts Inst. of Technology, Cambridge, Mass. 
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H. ALDEN FOSTER,! M. ASCE.—Mr. Childs’ paper brings out very clear- 
ly the fact that hydrological knowledge is never complete or final. The basic 
data for design of hydraulic structures need to be continually revised to con: 
form with new experience. Nature always demands the last word in such 
matters. ] 

The writer will confine his comments to the portions of the paper dealing 
with questions of probability or frequency of floods. j 

The author refers to the “probable maximum precipitation” which is used 
as a basis for deriving the “maximum possible flood”. This seems to invol} 
some confusion of ideas, although the latter expression has been commonly 
used in the past by the Corps of Engineers. It does not seem advisable to 
place any finite limit to future flood magnitudes, In any case, the “probable 


maximum precipitation” should be used to estimate the “probable maximum 
flood”. | 


1. Prin. Assoc., Parsons, Brinckerhoff, Hall & Macdonald, New York, N. ¥ 
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The difficulties of combining extraordinary floods with the “run-of-the- 
nill” variety when preparing a probability curve from a particular flood 
record were recognized by Allen Hazen nearly 30 years aga. (1) It is closely 
related to the question of “sampling errors” and “confidence limits”.(2,3) 

Any sample record of 25 years or more is likely to have one or more extreme 
values well outside the general trend of the rest of the sample. The writer 
1as been unable to develop any method by which the floods caused by extra- 
ordinary storms, such as hurricanes, could be separated from the record, 
2ecause some of the smaller floods may have been produced by hurricane 
storms at a considerable distance from the gaging station which would cause 
2xcessive floods in other locations. 

In preparing a flood-probability curve, one or two floods that obviously do 
10t conform with the general trend should be ignored in selecting the curve- 
9f-best-fit. This curve should only be used for such purposes as determining 
he mean annual flood damage where omission of the extreme floods would 
lave relatively small effect on the results because of their obviously small 
srobability. Of course the hurricane floods should not be ignored in spillway 
lesign, but here they should be supplemented with estimates of the “maximum 
srobable flood”. The probability of the latter is necessarily small, but we 
are not concerned with whether it is the 1,000-year flood or the 100,000-year 

ood; as long as we are convinced that it may occur sometime, or even be 
eo oica, we cannot afford to ignore it. 

These remarks do not include the so-called “historical floods”. It should 
e generally possible to include them in the probability study by the methods 
Beebe in the ASCE “Hydrology Handbook”, page 102ff. See also (4). 

The following comments are quoted from a report(5) on the Northeastern 
Floods of 1955 prepared under the supervision of the writer: 


“Climatic Changes. —The great increase in precipitation and stream 
discharge in the 1955 floods as compared with previously recorded 
storms indicates that the storms of August 1955 were of a very ex- 
ceptional nature. From a statistical standpoint, it appears that they 
would have a probability of occurrence much smaller than would be 
indicated by theoretical studies based on previously available records. 
Two alternative inferences may be drawn from this situation: 


(1) The 1955 storms were the result of atmospheric conditions of a 
most unusual character, following the simultaneous occurrence 
of a number of unfavorable controlling factors that normally 
would not be expected to occur together; or 

(2) There have been changes in general climatic or atmospheric con- 
ditions in recent years that tend to increase the frequency and 
intensity of storms of this type along the North Atlantic seaboard. 


“If the first inference is accepted, it would be possible to make a 
reasonable estimate of the probability of occurrence of such storms in 
the future, for use in determining the effect of the 1955 floods on the 
mean annual flood damage to be expected on rivers in the area under 
consideration. If the second inference holds, however, any estimates of 
flood probability based on previous records would have to be revised to 
include the effect of the assumed changes in meteorological conditions. 

_... The recent apparent changes in climate affecting the paths of hurri- 
canes along the North Atlantic seaboard cover too short a time period to 
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justify the assumption of any long-term major climatic changes in this — 
region. Nevertheless, the possibility of such cyclic variations in climate 
should not be neglected in estimating flood probabilities. 4 
“Probability of Occurrence of 1955 Floods.—If the question of cli- 
matic change is eliminated from consideration, the excessive magnitude” 
of the 1955 storms must be treated as due to an unusual combination of 
most unfavorable controlling factors, but excluding any additional 
factors that have not been effective during the period of record. Asa 
result, if we have prepared a suitable probability curve representing the 
previously recorded floods, the 1955 floods should be located on the 
same curve. But since these floods are so much larger than any previ- 
ously observed, such a method of plotting will result in assigning a very 
low value of probability to them. If the present hurricane situation is — 
assumed to be the result of climatic changes, the question of determin-— 
ing the probability of occurrence of the 1955 floods becomes much more 
difficult. Such a change would have to be assumed as commencing not — 
earlier than 1938, the year when the first severe hurricane in more than 
100 years struck the Long Island coastline. The intervening period of — 
17 years (between 1938 and 1955) is much too short to permit any reason 
able estimate of probability for the recent storms. All of this brings 
out the difficulties of estimating mean annual flood losses in any locali 
by statistical methods, even if the problem of errors of sampling could ~ 
be satisfactorily taken care of.” ) 


Although the 1955 floods in many places were exceptionally large in com- 
parison with previous recorded floods, their effect on the determination of 
the Mean Annual Flood Loss for any particular property may be relatively — 
small, Studies referred to above indicated that, for a number of properties 
in the City of Waterbury, Conn., the addition of the 1955 floods to the recore 
available for 1930-1953 would have increased the estimate of mean annual 
flood losses for those properties by less than 3 per cent. 
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